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1 . 0  BACKGROUND 


1.1  OBJECTIVE  OF  THIS  STUDY 


The  work  reported  here  was  performed  as  Phase  I  of  a  study  of 
applications  of  a  new  type  of  wing,  known  as  the  joined  wing.  The 
particular  application  studied  here  is  an  advanced  cruise  missile 
capable  of  being  launched  from  a  standard  21  inch  diameter  submarine 
torpedo  tube.  The  present  report  discusses  alternative  configurations 
for  joined  wing  cruise  missiles.  In  Phase  II  of  the  study  wind  tunnel 
tests  will  be  performed  on  the  most  promising  of  these  configurations. 

The  joined  wing  concept  has  previously  been  described  only  in  patents 
such  as  Ref.  1  (other  patents  are  pending),  and  in  reports  of  limited 
circulation,  e.g..  Refs.  2,  3,  4,  and  5.  Therefore,  Section  2  of  this 
report  presents  a  general  description  of  the  concept,  including  aerodyna¬ 
mic  and  structural  aspects.  The  theoretical  and  experimental  results  pre¬ 
sented  in  Section  2  were  obtained  prior  to  the  present  study.  Section  2 
also  summarizes  previous  work  on  related  configurations. 

Section  3  presents  some  special  requirements  of  cruise  missiles  and 
describes  how  these  relate  to  the  joined  wing  configuration.  Section  4 
discusses  one  key  aspect  of  cruise  missile  design,  i.e.,  the  requirement 
to  fly  close  to  the  surface  of  waves  or  terrain  to  minimize  the  proba¬ 
bility  of  being  detected  by  radar.  A  method  of  optimizing  the  missile 
design  for  closest  terrain-following  is  described. 

Section  5  presents  comparisons  of  cruise  performance  between  a 
series  of  conventional  advanced  cruise  missiles  and  their  counterparts 
employing  joined  wings.  Readers  interested  in  performance  rather  than 
guidance  and  control  should  read  Section  5  before  section  4. 

Conclusions  and  recommendations  are  given  in  Section  6. 

Appendices  A,  B,  and  C  present  detailed  material  relating  to  the 
optimum  terrain-following  described  in  Section  4. 
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2.0  BASIC  PRINCIPLES  OF  THE  JOINED  WING 

2.1  JOINED  WING  CONFIGURATIONS  AND  ADVANTAGES 

The  joined  wing  is  a  general  concept  which  involves  the  combination 
of  two  wings,  a  fuselage,  and  a  fin,  such  that  the  wings  form  a  diamond 
shape  both  in  plan  view  and  in  front  view.  Figure  1  shows  a  design  for 
a  joined  wing  manned  aircraft.  Figure  2  illustrates  a  design  for  an 
advanced  cruise  missile  employing  a  foldable  joined  wing  plus  a  canard. 
Many  arrangements  for  joining  the  wings  are  possible.  Figure  3  shows  a 
design  with  the  joint  located  inboard  of  the  tips,  while  Fig.  4  shows 
that  the  joint  can  be  a  winglet,  thus  returning  some  aerodynamic  benefit 
as  well  as  providing  the  necessary  structural  tie  between  the  wings. 

Advantages  claimed  for  the  joined  wing  include: 

•  Light  weight 

•  High  stiffness 

•  Low  induced  drag 

•  Good  transonic  area  distribution 

•  High  trimmed  Cl  max 

•  Reduced  wetted  area  and  parasite  drag 

•  Direct  lift  control  capability 

•  Direct  sideforce  control  capability 

•  Good  stability  and  control 

These  claims  have  been  supported  by  independent  analyses,  design 
studies  and  wind  tunnel  tests,  as  described  in  Sections  2.2.  through  2.9. 
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FIGURE  2.  ADVANCED  CRUISE  MISSILE  EMPLOYING  JOINED 


FIGURE  3.  INBOARD-JOINT  CONFIGURATION 


FIGURE  4.  JOINED  WING  AIRPLANE  WITH  WINGLETS 


2 . 2  STRUCTURAL  PRINCIPLES 


The  internal  structure  of  the  joined  wing  is  unconventional,  involv¬ 
ing  chordwise  tapering  of  wing  thickness,  as  shown  in  Fig.  5.  This 
feature  must  be  clearly  understood  in  order  to  apply  the  joined  wing  con¬ 
cept  so  that  it  yields  a  decisive  weight  advantage  over  the  conventional 
wing -plus -tail  arrangement.  The  necessity  for  chordwise  taper  of  skin 
thickness  is  explained  below. 

Consider  Fig.  6.  This  shows  that  the  lift  load  acting  upon  any 
chordwise  section  of  a  joined  wing  at  a  distance  y  from  the  aircraft's 
plane  of  symmetry  may  be  resolved  into  two  components, 

(1)  An  "inplane"  component,  AL  sin  0,  inclined  at 
the  local  "tilt  angle,"  6.  This  component  acts 
within  the  plane  of  the  truss  formed  by  the 
joined  wings. 

(2)  An  "out-of -plane"  component  AL  cos  0. 

The  inplane  component  is  well  resisted  by  the  truss  structure,  but  the 
truss  provides  no  special  advantage  in  resisting  the  out-of-plane  compo¬ 
nent.  Therefore  the  joined  wing  structural  material  must  be  disposed  so 
that  each  wing  resists  out-of -plane  loads.  This  is  achieved  by  concen¬ 
trating  the  wing  box  material  at  the  maximum  possible  distance  from  the 
local  plane  of  the  truss  (i.e.,  the  plane  x-x  in  Fig.  6).  This  increases 
the  wing's  effective  beam  depth  about  its  bending  axis,  as  shown  in 
Fig.  7. 

The  weight  advantage  resulting  from  chordwise  tapering  is  consid¬ 
erable.  This  has  been  demonstrated  by  Ref.  3,  which  compares  a  conven¬ 
tional  airliner  aluminum  wing-plus-tail  with  an  aluminum  joined  wing 
having  equal  airfoil  section  and  equal  parasite  and  induced  drag.  Both 
wings  had  taper  ratios  ct^p/croot  -  0.3.  Reference  3  employed  a  finite 
element  structural  analysis  program  to  optimize  both  wing  systems  for 
maximum  strength/weight.  Comparing  complete  weights  (including  control 
surfaces  and  non-structural  components)  it  was  found  that  the  complete 


SHADED  REGION  INDICATES  EQUIVALENT  SKIN  THICKNESS, 
ie:  AVERAGE  OF  LOCAL  SKINS,  STRINGERS  AND  SPARS. 
EQUIVALENT  SKIN  THICKNESS  IS  EXAGGERATED  FOR 
CLARITY,  TO  EMPHASIZE  THE  TAPER. 


FIGURE  5.  CHORDWISE  TAPER  OF  WING  SKIN  MATERIAL  THICKNESS 
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FIGURE  6.  IN-PLANE  AND  OUT-OF-PLANE  LIFT  COMPONENTS 


FIGURE  7.  LARGE  EFFECTIVE  BEAM  DEPTH/CHORD  RATIO  OF  JOINED  WING 


joined  wing  is  24%  lighter  than  the  complete  conventional  wing-plus- 
tail.  It  is  also  much  stiffer;  as  indicated  in  Fig.  8,  the  tip  deflec¬ 
tion  of  the  conventional  wing  is  2.8  times  that  of  the  joined  wing. 

Reference  3  did  not  attempt  to  optimize  the  taper  ratio  selected 
for  the  joined  wing.  Further  impro  'ements  are  possible  if  the  taper  of 
the  joined  wing  is  decreased  to  allow  for  the  loads  transmitted  through 
the  inter -wing  joint.  (In  general,  some  taper  should  be  retained,  unless 
the  cost  of  manufacture  dictates  otherwise.) 

Because  the  effective  beam  depth  is  primarily  determined  by  the 
chord  of  each  of  the  joined  wings  and  their  tilt  angle,  9,  it  is  possible 
to  employ  very  thin  airfoils.  Figure  9  shows  a  3%  thickness/ chord  ratio 
airfoil  which  has  an  effective  beam  depth  of  32%  of  its  chord.  Thus  low 
drag  airfoils  can  be  employed  for  supersonic  flight  without  the  struc¬ 
tural  weight  penalties  of  thin  cantilever  wings. 

A  secondary  benefit  of  the  chordwise  variation  of  skin  thickness 
shown  in  Figs.  7  and  9  is  that  it  thickens  the  skin  (a)  near  the  upper 
leading  edge,  where  airloads  are  high  and  where  a  good  aerodynamic  sur¬ 
face  contour  must  be  maintained,  and  (b)  around  the  flap  attachment 
region,  where  high  local  loads  due  to  flap  airloads  occur. 

For  minimum  wing  weight  a  high  tilt  angle  (6)  must  be  maintained 
over  most  of  the  span.  This  implies  that  a  tall  fin  and  relatively  short 
fuselage  is  a  desirable  combination. 

The  material  distribution  shown  in  Figs.  7  and  9  is  not  optimum  close 
to  the  wing  centerline  and  near  the  inter-wing  joint.  This  is  because 
torsional  and  other  secondary  structural  loadings  become  dominant  in  these 
regions.  Nevertheless,  the  basic  principles  described  above  have  been 
shown  (Refs.  3  and  5)  to  yield  joined  wing  structures  which  are  optimum 
over  most  of  their  spans. 

Buck! ing :  Due  to  the  inplane  load  components  the  rear  wing  is  in 
compression  during  flight  at  positive  "g"  load  factors.  Overall  column 
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VERTICAL  DEFLECTION 


FIGURE  8.  DEFLECTIONS  OF  JOINED  AND  CONVENTIONAL  WINGS  HAVING  EQUAL 
TOTAL  VERTICAL  AIRLOADS  AND  EQUAL  L/D’S,  WITH  JOINED  WING 
WEIGHT  -  0.76  x  CONVENTIONAL  WING-PLUS -TAIL  WEIGHT 
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TILTED  BENDING  AXIS 


AIRFOIL 

THICKNESS/CHORD  =  3% 


JOINED  WING  PROVIDES  EFFECTIVE  BEAM 
DEPTH  EVEN  FOR  VERY  THIN  AIRFOILS 


FIGURE  9.  CONSTRUCTION  OF  JOINED  WING  EMPLOYING  THIN  AIRFOILS 
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buckling  is  therefore  theoretically  possible.  However,  this  has  not 
proved  to  be  a  problem  on  any  configuration  studied  to  date.  If  buckling 
should  become  significant  (e.g.,  for  an  ultra -high-aspect-ratio  vehicle) 
it  could  be  alleviated  by  moving  the  joint  inboard,  as  in  Fig.  3. 

Torsion:  Torsion  of  the  front  joined  wing  is  resisted  by  its  own 
torsional  stiffness,  and  also  by  flexural  stiffness  of  the  rear  wing, 
as  shown  in  Fig.  10.  The  net  result  is  that  the  joined  wing  pair  is 
extremely  stiff  torsionally  and  flexurally. 

Wing  Box  Chord :  A  conventional  cantilever  box  spar  structure  extend¬ 
ing  from  15%  to  65%  chord  is  definitely  not  optimum  for  the  joined  wing. 
The  box  should  extend  from  the  leading  edge  (or  as  close  to  it  as  the 
de-icing  system  will  permit)  to  just  ahead  of  the  leading  edge  of  the 
flaps  and  ailerons. 

Fuselage  Structure  Weight:  The  joined  wing  is  a  special  form  of  the 
tandem  wing.  DeYoung  (Ref.  6)  and  Bottomley  (Ref.  7)  compare  fuselage 
weights  of  tandem  versus  conventional  configurations.  Both  these  refer¬ 
ences  conclude  that  for  a  given  total  lift,  and  with  a  fixed  fuselage 
length,  the  tandem  configuration's  fuselage  skin  is  lighter.  This  is 
because,  as  illustrated  in  Fig.  11,  the  two  uploads  of  the  tandem  wings 
impose  less  severe  bending  moments  than  do  opposing  wing  and  tail  loads. 
Figure  11  illustrates  the  resulting  saving  in  the  weight  of  fuselage 
structural  material. 

Since  minimum  fuselage  skin  gauge  is  determined  by  additional 
factors  such  as  pressure  loads  and  landing  loads,  only  part  of  the 
weight  saving  shown  in  Fig.  11  is  achieved.  However  DeYoung  (Ref.  6) 
shows  that  substantial  weight  reductions  can  be  achieved  despite  this 
limitation.  Compared  to  the  tandem  wing,  the  joined  wing  structure  has 
the  further  advantage  of  relieving  fuselage  torsional  and  lateral 
bending  moments. 

It  is  desirable  that  the  fuselage  should  be  short,  for  minimum 
weight.  As  shown  in  Section  2.7  this  is  compatible  with  the  aerodynamic 
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FIGURE  10.  BENEFICIAL  TORSIONAL /FLEXURAL  INTERACTION  OF  JOINED  WINGS 
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FUSELAGE  BENDING  MOMENT 


requirements  of  joined  wing  aircraft,  primarily  because  such  aircraft 
do  not  require  a  long  fuselage  to  provide  a  large  tail  moment  arm. 


2.3  INDUCED  DRAG  OF  JOINED  WINGS 


The  joined  wing  yields  substantially  less  induced  drag  than  a  con¬ 
ventional  monoplane  wing  of  equal  lift  and  span  operating  at  the  same 
dynamic  pressure.  This  has  been  demonstrated  by  Prandtl-Munk  biplane 
theory  and  supporting  wind  tunnel  tests  as  described  below. 

For  an  isolated  monoplane  wing  induced  drag  is  generally  expressed 
in  coefficient  form  as: 


cDi 


ireA 


(1) 


or  equivalently,  in  dimensional  form  as: 


Di 


L2 

irqb2e 


(2) 


2 

where  L  ■  lift;  q  =  dynamic  pressure,  (1/2) pV;  b  »  span;  and  the  coeffi¬ 
cients  Cj)^,  Cl  and  aspect  ratio  (*  b^/S)  are  all  referred  to  a  common 
area  S.  The  span-efficiency  factor  e  has  an  upper  limit  of  1  for  an 
isolated  undlhedralled  monoplane  wing,  but  as  shown  by  Prandtl  and  Munk 
(see,  e.g..  Ref.  8),  for  biplane  type  configurations  substantially  higher 
values  of  e  are  obtainable. 

Equations  1  and  2  are  valid  for  monoplane  wings  which  do  not  have 
spanwise  variations  of  incidence  or  camber  (i.e.,  "untwisted"  wings). 
However,  for  twisted  wings.  Ref.  9  shows  that  an  alternative  form  of 
induced  drag  equation,  Eq.  3,  is  more  accurate. 


CDi 


twisted  monoplane  wing 


(Cl-C^)2 

TTAex  +  (CDi)min 


(3) 


17 


where  Cl*  and  ex  are  constants.  Note  that  Eq.  3  yields  an  offset  drag 
polar  even  though  no  viscous  drag  is  included.  (Offset  drag  polars 
also  occur  for  cambered  airfoils  tested  in  two-dimensional  conditions, 
but  that  offset  is  due  to  viscous  effects  and  is  additional  to  the 
inviscid  induced  drag  effects  discussed  in  this  section.) 

Equation  3  is  generally  more  accurate  than  Eq.  1.  However  Prandtl- 
Munk  biplane  theory  has  been  developed  on  the  assumption  that  the  induced 
drag  of  each  wing  of  a  biplane  pair  can  be  described  by  Eq.  1,  when  that 
wing  is  isolated.  Therefore,  to  facilitate  reference  to  the  theory,  we 
shall  employ  Eq.  1  for  most  of  this  discussion,  explaining  at  the  end 
how  twist  and  camber  corrections  can  be  included. 

Prandtl-Munk  Biplane  Theory:  Since  the  fundamentals  of  biplane 
induced  drag  are  well  presented  in  Refs.  8,  10,  and  11,  only  a  brief 
summary  will  be  given  here  to  set  the  stage  for  a  specialization  of  the 
theory  for  the  joined  wing. 

The  induced  drag  of  biplane  wings  having  optimal  span  loadings  is: 


Di 


ML2  l| 

blb2  +  b2 


(4) 


where  Lj,  L2  are  the  lifts  on  wings  of  respective  spans  bj  and  b2, 
p  -  air  density,  V  -  speed,  and  o  is  a  constant,  determined  by  the 
Tref f tz-plane  configuration,  i.e.,  the  front  view  of  the  wings  as  seen 
by  an  observer  looking  parallel  to  free-stream  direction.  Tables  of  0 
for  biplane  wings  are  given  by  Laitone  in  Ref.  12,  which  corrects  errors 
in  similar  tables  given  in  Ref.  10. 

For  wings  of  equal  spans  bi  •  b2  “  b  with  wing  No.  1  carrying  a 
fraction  fj  of  the  total  lift,  W; 
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This  expression  is  minimized  by  putting  fj  =  0.5.  However,  the 
minimum  is  much  flatter  than  the  corresponding  minimum  for  a  wing-plus- 
tail  combination,  as  illustrated  in  Fig.  12.  The  excess  drag  above  the 
minimum  is  trim  drag.  Hence,  trim  drag  is  relatively  less  important  for 
an  equal  span  staggered  biplane  (i.e.,  a  tandem  wing  or  joined  wing  air¬ 
plane)  than  for  a  staggered  biplane  of  very  unequal  spans  (i.e.,  a  con¬ 
ventional  aft  tail  airplane) .  For  preliminary  design  purposes  it  will 
often  suffice  to  neglect  trim  drag  for  joined  wing  aircraft. 

Hunk's  stagger  theorem  (Refs.  8  and  12)  states  that  the  induced  drag 
of  a  biplane  pair  of  wings  carrying  lifts  Li  and  L2  maintaining  a  constant 
Tref f tz-plane  configuration  is  unaffected  by  displacements  of  the  surfaces 
parallel  to  the  freestream.  The  theorem  implies  that  sweep  need  not 
affect  induced  drag,  provided  span  is  constant.  Hunk's  stagger  theorem 
has  been  verified  by  Lockheed  boxplane  tests  (Ref.  13)  and  also  by  Vought 
tandem  wing  tests  (Ref.  14)  both  of  which  are  discussed  in  Section  2.9. 
Hoerner  (Ref.  15,  page  7-8)  presents  some  data  which  appear  to  contradict 
the  stagger  theorem,  showing  that  sweep  influences  induced  drag.  However, 
Hoerner  states  that  "swept  wings  (or  certain  other  wing  shapes)  can  have 
the  same  induced  drag  as  straight  wings  of  the  same  aspect  ratio  —  pro¬ 
vided  the  lift  distribution  is  the  same."  To  accomplish  this,  proper 
twist  and  camber  are  required,  as  was  incorporated  in  the  swept  tandem 
wing  wind  tunnel  model  described  in  Ref.  16.  The  induced  drag  of  this 
model  was  found  to  be  in  excellent  agreement  with  the  minimum  achievable 
induced  drag  predicted  by  theory. 

Letcher  (Ref.  17)  has  calculated  the  minimum  achievable  induced 
drag  of  joined  wings  having  optimum  twist  and  camber,  and  also  optimum 
relative  load-sharing  factor,  fj.  Kuhlman  (Ref.  18),  by  an  independent 
method,  verified  Letcher's  results  and  extended  them  to  include  joined 
wings  with  winglets  as  shown  in  Fig.  13.  Figure  13  indicates  that  a 
joined  wing  fitted  with  winglets  can  achieve  high  values  of  e.  For 
example,  with  d/h  -0.5  and  h/b  ■  0.321  (corresponding  to  ±20°  dihedral) 
e  -  1.328. 
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FIGURE  12.  INDUCED  DRAG  OF  CONVENTIONAL  AND  MULTI -WING 
CONFIGURATIONS  WITH  OPTIMAL  SPAN-LOADINGS 


FIGURE  13.  CALCULATED  SPAN-EFFICIENCY  FACTOR  FOR  JOINED  WING 
WITH  SYMMETRIC  INCLINED  WINGLETS  ("CARPET"  PLOT) 
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Corrections  to  Prandtl-Munk  Theory:  Even  for  untwisted  wings,  two 
corrections  must  be  applied  to  values  of  e  calculated  by  the  Prandtl- 
Munk  biplane  theory,  such  as  those  given  in  Fig.  13: 

(1)  A  correction  for  the  downward  drift  of  the  vortex 
sheet  shed  by  the  front  wing  (the  theory  does  not 
allow  for  this  drift) . 

(2)  A  correction  for  wing-body  interference  and  for 
loss  of  leading  edge  suction. 

The  first  correction  is  significant  for  configurations  having  stag¬ 
gered  wings  of  large  (vertical)  gap.  Tandem  wing  test  data  analysed  in 
Ref.  14  indicate  that  the  e  given  by  Fig.  13  should  be  increased  by  102 
to  14%  for  typical  joined  wing  configurations,  in  which  the  rear  wing  root 
is  higher  than  the  front  wing  root. 

The  second  correction  is  standard  for  conventional  configurations. 
Typically  e  is  reduced  by  102  to  202  at  low  and  moderate  Cl's.  For  the 
joined  wing  the  second  correction  should  be  less  than  for  a  conventional 
configuration  since  most  of  the  induced  drag  (typically  2/3  of  the  total) 
is  developed  on  the  rear  wing,  which  is  not  intersected  by  the  fuselage. 

At  high  Cl's  there  is  an  appreciable  loss  of  leading  edge  suction,  and 
the  second  correction  becomes  larger.  However,  at  Cl's  within  the  normal 
cruise  range  for  most  joined  wing  configurations  Corrections  (1)  and  (2) 
approximately  cancel,  so  Fig.  13  can  be  employed  directly  to  compute  Cd^. 

Twisted  Wings:  As  will  be  explained  later,  properly  designed  joined 
wings  have  spanwise  variations  of  incidence  and  camber.  Therefore  to 
compute  the  drag  polar  accurately  over  a  range  of  Cl  the  parameters  ex,  Clx, 
must  be  known.  These  parameters  can  be  found  via  vortex -lattice  methods 
(e.g.,  Refs.  19,  20)  by  computing  Cj^  for  a  range  of  Cl  and  graphing  the 
results.  Often  Cj^  is  required  only  over  a  small  range  of  Cl-  If  the 
wing  is  twisted  and  cambered  to  yield  minimum  Cj>^  at  some  specified 
design  Cl  then  in  the  region  of  the  design  Cl  it  should  be  sufficiently 
accurate  to  compute  Cj)^  from  Eq.  1  with  e  given  by  Fig.  13,  plus  the 
corrections  described  above. 
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Experimental  Validation:  The  predictions  of  the  Prandtl-Munk  theory 
relating  to  induced  drag  of  multi-wing  configurations  have  been  compared 
with  wind  tunnel  test  data  in  Refs.  13,  14,  and  16  and  have  been  shown 
to  be  in  good  agreement. 

2.4  PARASITE  DRAG 


Compared  to  the  conventional  aircraft  the  joined  wing  displays  some 
advantages  and  some  disadvantages  with  respect  to  parasite  drag. 

Adverse  effects  include: 


(1) 

Lower  wing  Reynolds  numbers. 

(2) 

Increased  ratio  of  wetted  area  to 
due  to  dihedral  cosine  effect. 

lifting  area 

Beneficial  effects  include: 

(1) 

Reduced  wing  area  due  to  high  trimmed  CTmav 
(see  Section  2.6). 

(2) 

Reduced  fuselage  wetted  area  due 
length. 

to  reduced 

(3) 

Improved  wing-fuselage  interference. 

(4) 

Smoother  fairing  of  landing  gear. 

(5) 

Fin-nacelle  integration. 

Wing  Reynolds  Numbers:  Comparing  a  conventional  configuration 
having  wing  area  S y  and  tail  area  Sx  with  a  joined  wing  configuration 
of  the  same  total  horizontal  projected  area  and  the  same  span  yields 
Reynolds  numbers  in  the  ratio  0.5[1  +  (Sj/Sy)].  Typically  this  yields 
joined  wing  Reynolds  numbers  which  are  60%  to  65%  those  of  the  wing  of 
the  conventional  airplane.  This  effect  is  adverse  to  the  joined  wing. 
A  possible  benefit  would  arise  if  laminar  flow  control  is  considered, 
since  it  is  easier  to  establish  laminar  flow  at  low  Reynolds  number. 
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Dihedral  Cosine  Effect:  Considerations  of  structure  weight  and 
induced  drag  favor  the  use  of  steep  dihedral  angles  for  joined  wings. 
Such  dihedral  angles  increase  wing-plus-tail  wetted  area  by  typically 
3.5%  compared  to  a  conventional  configuration. 

To  overcome  these  adverse  effects  joined  wing  configurations  must 
be  designed  to  take  full  advantage  of  the  benficial  effects  discussed 
below. 

High  Trimmed  Ct. may  :  This  is  discussed  at  length  in  Section  2.6. 
The  aerodynamic  design  principles  given  in  Section  2.6  permit  substan¬ 
tial  reductions  in  joined  wing  wetted  area  compared  to  wetted  area  of  a 
conventional  wing -plus -tail  of  similar  lift  capability. 

Reduced  Fuselage  Length:  The  joined  wing  airplane  can  obtain  ade¬ 
quate  static  stability  without  a  long  tail  moment  arm,  as  shown  in  Sec¬ 
tion  2.7.  Hence  a  short  fuselage  of  lower-than-usual  fineness  ratio  (FR) 
can  be  employed.  Such  a  fuselage  may  be  advantageous.  Figure  14  illus¬ 
trates  this  by  comparing  the  drag  of  fuselages  having  untapered  mid¬ 
sections.  At  subcritical  Mach  numbers,  for  equal  volumes,  the  drag  of 
the  FR  *  6  fuselage  is  only  91%  of  that  of  the  FR  -  9.7  fuselage.  With 
proper  nose  and  tail  shaping  the  short  fuselage  can  avoid  a  significant 
decrease  in  drag  divergence  Mach  number. 

If  fuselage  maximum  cross-sectional  area  is  fixed  (instead  of 
fuselage  volume)  the  drag  benefits  of  a  short  fuselage  are  even  more 
striking.  The  FR  *  6  fuselage  then  has  only  67%  of  the  drag  of  the 
FR  ■  9.7  fuselage. 

Improved  Wing-Fuselage  Interference;  Attaching  the  front  wing  to 
the  fuselage  ahead  of  the  maximum  fuselage  cross-section,  as  shown  in 
Fig.  15,  places  the  wing  root  in  a  favorable  pressure  gradient.  This 
reduces  the  boundary  layer  thickness  at  the  wing-fuselage  junction  and 
increases  wing  root  local  Cl  max  . 

Smoother  Landing  Gear  Fairing  :  As  explained  in  Section  2.8,  the 
malnwheels  of  a  joined  wing  airplane  can  be  located  well  aft  of  the  c.g. 
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FIGURE  15.  WING-FUSELAGE  INTERFERENCE  COMPARISON 
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Indeed  the  main  wheels  can  be  retraced  into  the  tail  cone,  thus  avoiding 
the  unsightly  and  drag-producing  "landing  gear  bump"  that  characterizes 
many  current  transports. 

Fin-Nacelle  Integration:  If  properly  designed,  integration  of 
engine  nacelles  with  the  fin  can  save  wetted  area,  by  elimination  of 
pylons  and  masking  of  fin  area.  Fin-nacelle  integration  is  well  suited 
for  joined  wing  aircraft  because  the  fin  is  not  a  cantilever  structure. 
Thus  more  than  one  engine  can  be  integrated  with  the  fin. 

2 . 5  WAVE  DRAG 


For  transonic  and  supersonic  flight  special  emphasis  should  be  given 
to  those  aspects  of  the  joined  wing  which  yield  low  wave  drag.  Appro¬ 
priate  guidelines  for  design  are  as  follows. 

Wave  Drag  at  Zero  Lift:  As  indicated  by  Fig.  16  the  joined  wing 
airplane  has  an  exceptionally  smooth  area  distribution  around  M  *  1.0. 

It  is  therefore  predicted  to  have  a  lower  transonic  drag  rise  than  a 
conventional  configuration.  Test  data  supporting  this  prediction  are 
analyzed  in  Ref.  2,  p.  13.  The  data  were  obtained  from  transonic  tests 
on  a  boxplane  configuration  (Ref.  13)  described  in  Section  2.9. 

A  most  important  feature  of  the  joined  wing  is  its  compatibility 
with  very  thin  airfoils  (see  Fig.  9).  As  shown  in  Ref.  15  (Sec.  17, 
p.  12)  and  Ref.  21,  p.  277,  airfoil  transonic  wave  drag  coefficient  is 
approximately  proportional  to  (thickness/chord)  ,  so  dramatic  drag 
reductions  can  be  achieved  by  exploiting  this  aspect  of  the  joined  wing. 

Wave  Drag  Due  to  Lift :  As  explained  in  Ref.  22,  wave  drag  due  to 
lift  is  minimized  by  distributing  the  lifting  surfaces  over  a  large 
fraction  of  the  total  airframe  length.  This  design  objective  is  readily 
achieved  using  the  joined  wing. 

Wing  Interference  Drag:  The  configurations  of  Figs.  1,  2,  3,  and  4 
all  employ  a  wing  joint  arrangement  with  considerable  stagger,  so  that 

!  j 
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Che  wings  do  not  overlap  in  plan  view.  This  stagger  is  important:  at 
supercritical  speeds  it  eliminates  the  possibility  of  shock  waves  from 
one  wing  impinging  on  the  other,  with  a  consequent  shock- induced  sepa¬ 
ration  and  increase  of  drag.  Staggering  the  wings  at  the  joint  is  also 
desirable  at  subcritical  conditions  to  avoid  the  increase  of  boundary 
layer  thickness  that  otherwise  would  occur  at  the  junction  of  two  air¬ 
foils. 

Figures  1,  2,  3,  and  4  show  only  a  few  of  many  possible  wing  joint 
arrangements.  Others,  including  "wraparound"  continuous  curved  tips, 
and  "slat"  type  arrangements  are  discussed  in  Ref.  2. 

2.6  DESIGN  FOR  HIGH  LIFT 


Certain  guidelines  must  be  followed  to  enable  the  joined  wing  to 
obtain  high  lift  without  excessive  wing  area.  These  guidelines  include 
considerations  of  trim,  as  well  as  flow  separation.  This  subsection 
first  explains  the  requirements  for  obtaining  high  trimmed  lift  by 
selecting  appropriate  overall  airframe  geometry.  Airfoil  and  flap 
design  considerations  are  discussed  subsequently. 


Trim  Equations  and  Their  Consequences  for  Airframe  Geometry: 

Figure  17  shows  the  lift  forces  Lp,  Lr,  acting  at  the  aerodynamic  cen¬ 
ters  of  a  joined  wing  pair.  When  each  wing  (individually)  is  producing 


zero  net  lift  it  produces  a  pitching  couple  Mg,  as  indicated.  For  1  g 
steady  flight  the  lift  and  moment  equations  about  the  c.g.  become  (with 
notation  as  in  Fig.  17) 


LpS-p  +  MgF  +  MgR  -  Lr£r 


(6) 


Lp  +  Lr  -  W 


(7) 


It  will  be  found  for  almost  all  statically  stable  joined  wing  configu¬ 
rations  that  when  the  front  wing  attains  its  the  rear  wing  lift 
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required  to  maintain  trim  is  less  than  the  rear  wing  could  generate  were 
trim  not  required.  In  other  words:  "The  front  wing  stalls  first." 

This  provides  good  stall  recovery  but  nevertheless  it  is  undesirable 
for  the  front  wing  to  stall  while  the  rear  wing  is  still  far  from  stall¬ 
ing.  If  this  occurs  it  implies  that  the  rear  wing  area  is  oversized 
from  the  viewpoint  of  generating  the  total  maximum  lift  required  to  meet 
the  requirements  of  1  g  and  maneuvering  flight.  One  would  prefer  the 
front  wing  to  stall  when  the  rear  wing  is  almost  at  its  Cy. maY  .  To 
achieve  this,  certain  principles  must  be  followed  in  designing  the  over¬ 
all  airframe  geometry,  especially  the  planform.  These  principles  can  be 
derived  by  manipulating  Eqs.  6  and  7  as  follows. 

From  Eq.  6, 


whence: 


lR 


MoF  +  ^R 


(7) 


cLr  _  sp  .  £f  +  (mqf  +  mqr) 
cLf  SR  %  CLF  '  sr£r 


The  first  term  in  the  right  side  of  Eq.  8  is  always  less  than  unity  for 
a  statically  stable  aircraft.  Even  for  identical  front  and  rear  wing 
areas  and  equal  (positive  and  negative)  sweep  angles  £p/£R  <  1,  because 
the  overall  vehicle  aerodynamic  center  is  nearer  to  the  a.c.  of  the 
(isolated)  front  wing  than  the  a.c.  of  the  (isolated)  rear  wing.  This 
is  because  the  CLa  of  the  rear  wing  is  (1  -  d e/da)  times  its  isolated 

The  following  principles  should  be  followed  to  obtain  the  highest 
possible  value  of  the  first  term  in  Eq.  8. 

(1)  Decrease  the  ratio:  (frontwing  Clq) / (rearwing  CLq) 
by  employing  less  sweep  on  the  rear  wing. 

(2)  Select  front  and  rear  wing  airfoils  such  that  the 
airfoil  with  the  highest  two-dimensional  CLa  is  on 
the  rear  wing. 
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(3)  Increase  Che  ratio  of  rear  wing/front  wing  chord  to 
move  the  a.c.  aft. 

(4)  Reduce  the  effective  de/da  at  the  rear  wing  by  using 
a  large  centerline  gap  (i.e.,  make  the  fin  tall).  It 
is  also  beneficial  to  employ  winglets,  since  these 
reduce  the  average  de/da  at  the  rear  wing. 

(5)  Take  advantage  of  the  large  pitch  damping  of  the  joined 
wing  configuration  to  move  the  c.g.  aft  without  sacri¬ 
ficing  maneuver  margin.  (Generally  maneuver  margin  is 
the  key  stability  parameter,  not  static  margin.)  If 
control  configured  vehicle  (CCV)  technology  is  allow¬ 
able,  it  should  certainly  be  employed,  so  that  the 
maneuver  margin  can  be  made  negative.  This  increases 
trimmed 

Turning  now  to  the  second  term  in  Eq.  8,  the  key  consideration  is 
to  make  the  Mq  terms  large  and  positive.  This  leads  to  the  following 
guidelines : 

(1)  Do  not  employ  equal  flap  deflections  across  the  span  of 
each  wing.  Instead  employ  maximum  flap  deflections  on 
the  inboard  section  of  the  front  wing  and  near  the  out¬ 
board  section  of  the  rear  wing,  as  shown  in  the  lower 
half  of  Fig.  17. 

(2)  Leading  edge  devices  near  the  root  of  the  front  wing  are 
beneficial. 

(3)  A  low  aspect  ratio  canard  can  provide  appreciable  nose-up 
moment  with  little  forward  a.c.  shift.  (This  is  discussed 
in  more  detail  in  Section  2.7.) 

(4)  Where  practicable,  the  fuselage  should  be  shaped  (i.e., 
cambered)  and  the  wing  incidence  selected  to  provide  posi¬ 
tive  Chq. 

The  foregoing  recommendations  have  emphasized  overall  airframe 
geometry  rather  than  local  aerodynamic  design  (e.g.,  of  airfoils  and 
flaps) .  It  is  however  particularly  important  that  the  front  wing  root 
airfoil  should  be  selected  to  give  a  high  local  as  discussed 

below. 
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Airfoil  Design  for  High  Lift.  Figure  18  (reproduced  from  Ref.  23, 
p.  16-4)  shows  how  the  outward  flow  in  the  boundary  layer  of  a  swept- 
back  wing  increases  near  the  wing  root.  This  effect  is  doubly 

beneficial  for  the  joined  wing  because: 

(1)  It  increases  Cl  max  directly  on  the  front 
wing. 

(2)  It  provides  a  nose-up  moment  which  allows 
the  rear  wing  to  be  trimmed  at  a  higher  Cl. 

It  is  desirable  to  design  joined  wings  such  that  this  effect  is 
maximized.  As  shown  by  experimental  data  graphed  in  page  3  of  Sec.  16 
of  Ref.  23  and  by  Fig.  19,  the  magnitude  of  this  effect  is  correlated 
with  the  sharpness  of  the  wing  leading  edge.  A  sharp-nosed  airfoil  is 
beneficial,  yielding  a  22%  increase  in  Cl  max  with  40°  sweep. 

Design  guidelines  for  airfoil  selection  thus  include: 

(1)  Near  the  front  wing  root  use  a  sharp-nosed  air¬ 
foil,  or  a  sharp-nosed  segmented  Kruger  flap  to 
generate  LE  vortex-shedding  which  will  assist 
sweeping  the  boundary  layer  outboard. 

(2)  Alternatively,  employ  a  leading  edge  root  exten¬ 
sion  and/or  a  close-coupled  canard.  The  canard 
should  be  of  low  aspect  ratio  to  avoid  excessive 
forward  shift  of  the  a.c.  which  (for  a  statically 
stable  aircraft)  will  limit  the  lift  that  can  be 
trimmed  on  to  the  rear  wing.  As  discussed  in 
Section  2.7  a  highly  swept  canard  with  a  sharp 
leading  edge  is  particularly  desirable.  Use 
should  be  made  of  existing  data  on  close-coupled 
canard  configurations  (e.g..  Ref.  25). 

Induced  Camber:  Addoms  (Ref.  26)  has  investigated  airfoil  selection 
for  biplanes.  He  shows  that  due  to  the  curvature  of  the  flow  induced  by 
each  wing  of  a  biplane  pair,  biplane  airfoils  require  increased  camber  to 
achieve  equal  to  that  attained  by  monoplane  airfoils.  For  joined 


33 


FIGURE  19.  INCREASE  OF  Cl  max  WITH  SWEEP  ANGLE 
(REPRODUCED  FROM  REF.  14) 


wings  the  induced  flow  curvature  varies  because  of  the  spanwise  varia¬ 
tion  in  gap  and  stagger  of  the  airfoils.  Thus,  as  shown  in  Fig.  20  the 
sign  of  the  flow  curvature  may  change  depending  on  the  spanwise  ordinate 
of  the  joined  wing  airfoils. 

The  induced  camber  effect  makes  it  impossible  to  perform  fair  com¬ 
parisons  between  monoplanes  and  multi-wing  configurations  if  identical 
airfoils  are  employed  for  both.  For  example,  in  the  comparative  mono¬ 
plane  vs.  boxplane  tests  described  in  Ref.  13,  identical  airfoils  were 
employed  for  both  configurations.  As  a  result,  the  boxplane  model 
suffered  from  premature  flow  separation  on  the  front  wing. 

Since  the  induced  camber  varies  with  Mach  number  some  compromise 
may  be  necessary  to  obtain  high  at  both  high  and  low  Mach  numbers. 
Multi-element  and  slat  airfoil  design  techniques  are  appropriate,  and 
the  vortex-lattice  program  of  Ref.  20  which  yields  optimum  joined  wing 
twist  and  camber  lines  for  specified  pressure  distributions  has  proven 
to  be  useful  in  airfoil  selection. 

Apart  from  induced  camber,  multi-plane  airfoil  systems  experience 
mutually  induced  changes  in  the  magnitude  of  local  airspeeds.  As  ex¬ 
plained  in  Ref.  27  these  can  be  employed  with  advantage  in  the  design  of 
two-dimensional  airfoil-plus-slat  combinations,  to  yield  high  Cl's.  If 
desired,  it  should  be  possible  to  obtain  corresponding  advantages  by 
overlapping  the  tips  of  joined  wings  partially  so  that  the  front  wing 
tip  acts  as  a  slat  for  the  rear  wing  tip. 

2.7  STABILITY  AND  CONTROL 

Longitudinal  Stability:  Figure  21  shows  the  location  of  the  aero¬ 
dynamic  center  for  typical  joined  wing  configurations  as  predicted  by 
the  vortex-lattice  program  of  Ref.  19.  Figure  21  also  shows  C,^,  Cl^, 
CLa  and  C^  for  one  configuration,  as  computed  by  the  above-mentioned 
vortex-lattice  program.  The  pitch  damping  1/4  pSUc^  is  more  than 
double  that  of  a  conventional  wing-plus-tail  configuration.  This  large 
damping  causes  the  maneuver  margin  to  exceed  the  static  margin  by  a 


FIGURE  20.  INDUCED  CAMBER  EFFECTS  OF  REAR  WING 
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FIGURE  21(a).  JOINED  WING  AERODYNAMIC  CENTER:  (a)  GEOMETRY 
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FIGURE  21(b).  JOINED  WING  AERODYNAMIC  CENTER:  (b)  RESULTS 
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greater  distance  than  is  typical  of  conventional  aircraft.  Hence  the 
joined  wing  can  achieve  satisfactory  longitudinal  stability  and  control 
with  very  small  (or  even  zero)  static  margin. 


Wind  tunnel  tests  on  tandem  and  joined  wing  configurations  (Refs.  4, 
14)  show  reasonable  agreement  with  the  a.c.  location  predicted  by  vortex- 
lattice  methods  at  low  and  moderate  Cl's.  However  at  high  Cl's  the  a.c. 
moves  aft  and  a  moderate  "pitch  down"  tendency  occurs.  (No  "pitch-up" 
trend  has  been  observed.)  The  pitch-down  is  of  value  as  a  stall  preven¬ 
tative,  but  has  the  disadvantage  that  by  moving  the  a.c.  aft  the  static 
margin  may  become  excessive.  This  is  undesirable  because  it  implies 
large  control  deflections  for  maneuvering  and  also  because  it  limits 
trimmed  Cr ,  as  discussed  previously.  Therefore,  if  a  canard  is 
employed  it  should  be  of  low  aspect  ratio  and  high  sweep  with  a  sharp 
leading  edge,  to  give  a  large  nonlinear  (vortex  lift)  contribution  to 
CLq.  This  nonlinear  CLa  component  can  be  tailored  to  correct  the  rear¬ 
ward  shift  of  the  a.c.,  thus  giving  essentially  constant  static  margin 
up  to  the  stall. 

Stall  Characteristics:  Flight  tests  of  joined  wing  models  and  of  a 
manned  joined  wing  glider,  plus  wind  tunnel  model  tests,  indicate  good 
stall  characteristics  with  no  tendency  for  wing-dropping . 

Lateral  Stability:  By  tailoring  the  positive  and  negative  dihedral 
angles  of  the  front  and  rear  wings  almost  any  desired  level  of  dihedral 
effect  (C^g)  can  be  obtained.  Joined  wing  models  show  good  spiral  sta¬ 
bility,  and  display  no  tendency  to  Dutch  roll.  Directional  stability 
also  appears  to  be  good,  probably  because  of  the  endplate  effect  of  the 
rear  wing  on  the  fin. 

Control :  Versatile  control  capabilities  can  be  obtained  if  each 
wing  half  is  equipped  with  a  control  surface.  As  shown  on  Fig.  22,  this 
can  provide  pure  pitching  moment,  direct  lift,  rolling  moment,  and  direct 
side-force.  Any  desired  combination  of  these  can  be  obtained  via  suit¬ 
able  coupling  of  the  individual  control  surfaces. 


PITCH  CONTROL 


\ 


DIRECT  LIFT  CONTROL 


DIRECT  SIDE  FORCE  CONTROL 


FIGURE  22.  VERSATILE  CONTROL  CAPABILITIES  OF  JOINED  WING  AIRPLANE 
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To  demonstrate  one  advantage  of  this  versatile  control  arrangement, 
consider  the  control  required  for  take-off  rotation.  By  suitable  choice 
of  front/rear  elevon  gearing  the  longitudinal  control  can  be  adjusted  to 
yield  a  combination  of  direct  lift  and  pure  pitching  moment  control.  Hence 
the  airplane  can  be  rotated  for  take-off  without  any  net  download  on  the 
rear  wing.  This  permits  the  main  wheels  to  be  located  further  aft  of  the 
c.g.  than  on  a  conventional  airplane,  at  which  location  it  can  be  conveni¬ 
ently  retracted  into  a  part  of  the  fuselage  where  it  does  not  impinge  upon 
passenger  or  cargo  volume.  For  example  the  main  wheels  could  be  retracted 
into  the  "dead -space"  of  the  fuselage  tail  cone. 

This  form  of  control  also  yields  faster  response  and  more  accurate 
control  of  the  flight  path  because  it  avoids  the  initial  wrong-way  height 
response  of  aft  tail  airplanes. 

The  direct  side  force  control  shown  on  Fig.  22  will,  in  general,  pro¬ 
duce  a  side  force  which  does  not  act  exactly  through  the  c.g.  Hence  some 
rudder  deflection  will  be  required  to  eliminate  the  residual  yawing  moment. 
Nevertheless  this  type  of  control  is  valuable  as  a  vernier  adjustment  for 
precise  alignment  of  the  flight  path,  e.g.,  on  the  approach.  It  may  also 
be  of  use  for  remotely  piloted  vehicles  where  the  remote  operator  may  be 
confused  by  the  requirement  to  bank  in  order  to  turn.  However,  we  do  not 
advocate  this  type  of  control  for  manned  aircraft,  where  all  four  elevons 
should  be  used  as  shown  in  Fig.  22  to  produce  a  rolling  moment  in  order  to 
initiate  a  turn  in  the  usual  manner. 

A  very  effective  airbrake  can  be  obtained  via  the  use  of  a  split 
rudder  located  near  the  trailing  edge  of  the  rear  wing  (Fig.  23).  The 
split  rudder  increases  parasite  drag  and  also  increases  induced  drag  by 
degrading  the  lift  distribution  on  the  rear  wing.  Any  resulting  nose-up 
pitching  moment  can  be  counteracted  by  elevons  or  spoilers  on  the  front 
wing,  or  outboard  elevons  on  the  rear  wing. 

Wing  Warping :  For  some  applications  (e.g.,  remotely  piloted  vehicles) 
wing  warping  control  may  be  of  interest,  because  it  avoids  the  expense  of 


separate  control  surfaces,  and  also  avoids  the  increase  in  radar  signature 
associated  with  control  hinge  lines.  For  the  joined  wing  such  control  can 
be  simply  obtained  via  actuators  inside  the  fuselage  which  twist  (for 
example)  the  front  wings  against  the  tip  constraint  supplied  by  the  rear 
wings . 

2.8  SPECIAL  DESIGN  ASPECTS 

Under  this  heading  we  discuss  some  specialized  topics  such  as  joined 
wing  landing  gear  design,  crashworthiness,  folding,  and  suitability  for 
VTOL  aircraft. 

Landing  Gear:  The  joined  wing  airplane  having  elevons  on  both  wings 
can  rotate  for  take-off  without  incurring  a  download  on  the  rear  wing. 
Hence  the  main  wheels  can  be  located  far-*ift  of  the  c.g.,  possibly  even 
retracting  into  the  tail  cone.  Thus,  as  shown  in  Fig.  24  the  legs  of  the 
main  landing  gear  can  be  shorter  than  those  of  a  conventional  airplane, 
while  still  providing  the  necessary  ground  clearance  during  rotation.  A 
short  landing  gear  is  usually  a  light  landing  gear,  and  retracting  it  into 
the  tail  cone  obviates  the  familiar  drag-producing  "landing  gear  bulge." 

The  feasibility  of  this  type  of  landing  gear  depends  upon: 

(1)  Overturning  criteria,  which  demand  a  wider  track 
for  main  wheels  mounted  far  aft  of  the  c.g. 

(2)  Directional  stability  during  take-off  and  landing, 
when  the  airplane  may  be  pivoting  in  yaw  about  its 
main  wheels. 

Crashworthiness :  The  joined  wing  airplane  offers  a  potential  advan¬ 
tage  in  this  regard.  Placing  the  front  wing  root  towards  the  nose  of  the 
fuselage  enables  the  joined  wings  to  brace  the  passenger  compartment  so 
that  it  is  less  likely  to  buckle  in  a  crash.  The  front  wing  becomes  the 
major  energy  absorbing  structural  element. 

Ditching  characteristics  of  joined  wing  aircraft  are  expected  to  be 
satisfactory,  since  the  airplane  floats  supported  by  the  front  wing,  wing- 
lets,  and  tip  tanks. 
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CONVENTIONAL  AIRPLANE 


JOINED  WING  AIRPLANE 


CLEARANCE 


FIGURE  24. 


EFFECT  OF  MAINWHEEL/C.G.  LOCATION 
ON  LANDING  GEAR  LENGTH 
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Folding :  Figure  25  illustrates  a  folding  system  suitable  for  manned 
aircraft.  Other  systems  have  been  devised  for  missile  applications  where 
it  is  desired  to  accommodate  the  folded  wing  within  the  missile  fuselage.* 

Fuel  Volume:  For  a  given  thickness/chord  ratio  the  volume  of  a 
joined  wing  is  less  than  that  of  a  conventional  wing  of  similar  perfor¬ 
mance.  Thus  to  ensure  adequate  fuel  volume  within  the  wing  structural 
box  it  should  extend  over  as  large  a  fraction  of  the  chord  as  possible. 
(This  is  desirable  for  structural  reasons  also,  as  shown  in  Section  2.2.) 
The  tip  joining  fairing  can  also  serve  as  a  tip  tank.  Joined  wing  fuel 
volume  has  not  proved  to  be  a  significant  constraint  on  any  configuration 
studied  to  date. 

Suitability  for  VTOL;  The  joined  wing  arrangement  is  well  suited  to 
a  VTOL  aircraft  employing  two  main  lifting  units  displaced  laterally  from 
the  c.g.  This  minimizes  interference  between  the  wings  and  the  lifting 
units,  which  may  be  ducted  fans  or  propellers. 

2 . 9  PUBLISHED  WORK  ON  RELATED  CONFIGURATIONS 


The  joined  wing  concept  involves  complete  integration  of  structural 
and  aerodynamic  aspects.  In  this  it  differs  from  the  related  concepts 
that  have  been  published.  For  example,  the  Lockheed  boxplane  (Fig.  26) 
described  in  Refs.  13  and  28  provides  low  induced  drag,  but  has  less 
structural  rigidity  than  the  joined  wing,  as  indicated  by  the  comparative 
front  views  of  Fig.  27.  Reference  28  shows  that  the  low  rigidity  of  the 
boxplane  yields  unacceptable  flutter  characteristics.  The  rhomboidal 
wing  of  Zimmer  (Ref.  29),  shown  in  Fig.  28,  achieves  a  more  rigid  struc¬ 
ture  than  the  boxplane  but  lacks  the  structural  and  aerodynamic  benefits 
obtainable  through  integrating  the  fin  with  the  wing,  and  through  the  wing 
structural  thickness  distribution  shown  in  Fig.  5.  Thus  the  rhomboidal 
wing  configuration  has  area  rule  and  induced  drag  characteristics  which 
are  only  slightly  better  than  those  of  a  conventional  configuration. 


* 


A  cruise  missile  folding  system  is  described  in  Section  3. 
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:  442  58  m*  (4764  ft*) 


BOTH  UPPER  AND  LOWER 
WINGS  ARE  SWEPTBACK 


VIA  PYLONS 


FIGURE  28.  RHOMBOIDAL  WING  CONFIGURATION 


Reference  16  describes  wind  tunnel  tests  on  a  boxplane  with  and 
without  endplates  joining  the  wings.  Some  additional  analysis  of  the 
data  is  given  in  Ref.  14. 

The  wind  tunnel  data  presented  in  Refs.  13,  14,  28  and  29  show  good 
stability  and  control  characteristics  and  lift  and  drag  characteristics 
which  generally  agree  with  theoretical  predictions.  The  boxplane  tran¬ 
sonic  test  results  of  Ref.  13  show  premature  separation  on  the  front  wing 
at  high  Cl's:  as  explained  in  Section  2.6  this  can  be  ascribed  to  the 
airfoil  selected  being  the  same  as  that  of  a  "baseline"  comparable  mono¬ 
plane,  thus  having  insufficient  camber  to  compensate  for  the  "induced 
camber"  effect  of  the  boxplane.  (Further  analysis  of  boxplane  data  is 
given  in  Ref.  2.) 

Dexter  and  Cahill  (Ref.  3 CO  present  tests  at  Mach  numbers  between 
0.2  and  0.9  on  two  very  similar  configurations,  one  of  which  is  shown  in 
Fig.  29.  These  configurations  had  negative  tilt  angles  and  contained 
numerous  crudities  yielding  high  drag  (e.g.,  too  sharp  wraparound  radius, 
flat  plate  airfoil).  As  a  result  the  measured  span-efficiency  factor  for 
the  configuration  of  Fig.  28  was  only  0.502.  For  these  reasons,  plus  the 
configurational  differences,  the  tests  of  Ref.  30  have  very  limited  rele¬ 
vance  to  the  joined  wing  concept. 

Apart  from  the  data  referenced  above,  the  only  other  significant  data 
that  we  have  been  able  to  locate  on  related  configurations  is  an  unpub¬ 
lished  wind  tunnel  test  report  on  the  Warren-Young  wing  (Fig.  30).  This 
figure  shows  one  of  the  wind  tunnel  models  tested  at  low  Mach  number  by 
the  University  of  Bath,  England.  A  manuscript  copy  of  the  unpublished 
test  report  (Ref.  31)  was  kindly  supplied  by  Dr.  J.  S.  Henderson  of  that 
institution.  The  Warren-Young  project  differs  from  the  joined  wing  air¬ 
plane  because  it  does  not  take  advantage  of  the  height  of  the  vertical 
tail  to  obtain  a  large  tilt  angle.  A  further  difference  is  that  the 
Warren-Young  arrangement  connects  the  wings  by  means  of  a  large  horizon¬ 
tal  tip  surface.  The  wind  tunnel  tests  of  Ref.  31  confirmed  previous 
radio-controlled  model  tests  which  indicated  good  flying  qualities  and 
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gentle  stall,  but  showed  a  nonlinear  Cj*  -  variation  below  the  stall. 
Between  Cl's  of  0  and  1.0  (based  on  total  wing  area),  the  aerodynamic 
center  moved  aft  by  approximately  35%  of  the  geometric  mean  chord  (i.e., 
total  area/ span) .  By  contrast,  the  corresponding  low  Mach  number  boxplane 
tests  of  Refs.  13  and  16  showed  very  little  aerodynamic  center  shift. 

Sketches,  photos,  and  drawings  of  other  configurations  related  to 
the  joined  wing  have  been  published,  but  the  accompanying  data  have  been 
fragmentary  and  incomplete.  The  references  cited  above  are  therefore 
believed  to  represent  the  principal  sources  of  technical  information  on 
related  configurations. 

Currently  several  research  programs  on  the  joined  wing  are  being 
planned.  Additional  data  including  further  wind  tunnel  results  are 
explected  to  be  available  soon.  It  is  therefore  suggested  that  the 
author  of  this  report  should  be  contacted  to  obtain  the  most  complete 
and  up-to-date  information. 


3.0  JOINED  WING  CRUISE  MISSILE  CONFIGURATIONS 

3.1  INTRODUCTION:  CRUISE  MISSILE  DESIGN  REQUIREMENTS 

Having  described  the  basic  principles  of  the  joined  wing  in  Sec¬ 
tion  2,  this  section  applies  these  principles  to  the  design  of  cruise 
missiles.  Problems  and  solutions  relating  to  overall  configuration 
selection  and  mechanical  design  are  discussed:  detailed  performance 
and  control  analyses  are  given  later,  in  Sections  4  and  5  of  this 
report . 

For  the  purposes  of  this  report,  a  "cruise  missile"  is  defined  as 
a  missile  which  spends  the  major  portion  of  its  flight  in  a  "cruise"  mode, 
i.e.,  flying  at  nearly  constant  altitude  and  nearly  constant  speed,  using 
aerodynamic  lift  to  support  its  weight.  The  advanced  cruise  missile  con¬ 
figurations  discussed  in  this  report  have  launch  weights,  payloads,  and 
folded  dimensions  similar  to  those  of  Tomahawk  class  vehicles.  The  speci¬ 
fied  low-altitude  cruise  flight  condition  (M  *  0.7  at  sea  level)  is  also 
similar  to  that  of  the  Tomahawk.  However,  the  present  study  is  not 
restricted  to  vehicles  comparable  with  the  Tomahawk;  for  example,  we 
consider  configurations  having  Tomahawk-like  fuselages  but  wings  which 
are  two  or  three  times  larger  in  area  than  those  of  the  current  Tomahawk, 
to  improve  maneuverability. 

General  design  requirements  for  cruise  missiles  have  been  well  sum¬ 
marized  in  Refs.  32  and  33.  From  these  references,  plus  discussions  with 
personnel  of  the  Joint  Cruise  Missiles  Project  Office,  we  formulated  the 
following  design  goals  for  an  advanced  low-altitude  cruise  missile  capable 
of  being  launched  from  a  standard  21  inch  submarine  torpedo  tube. 
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Principal  Design  Goals 


•  Low  wave  (or  terrain)  clearance  to  avoid 
detection 

•  Low  radar  signature  to  avoid  detection 

•  Robustness  of  the  folded  vehicle 

•  Reliable,  rapid,  and  symmetric  unfolding 

•  Simplicity 

•  Low  cost 

•  High  sustained  "g"  capability  for  terminal 
maneuvers 

•  Long  range 

•  Cruise  at  M  ■  0.7,  little  emphasis  on  tran¬ 
sonic  and  supersonic  capability 

The  design  goals  listed  above  serve  as  criteria  by  which  the  rela¬ 
tive  merits  of  joined  wing  and  conventional  cruise  missiles  can  be 
compared . 

3.2  COMPARATIVE  AIRFRAMES 

In  order  to  compare  the  joined  wing  versus  the  conventional  wing- 
plus-tail  arrangement,  3  configurations  of  each  type  were  selected.  The 
conventional  configurations  (denoted  3GA,  6GA,  9GA)  had  gross  wing  areas 
of  respectively  11.5,  23.0,  and  34.5  ft^,  as  shown  in  Fig.  31.  The  joined 
wing  configurations  (denoted  3F,  6F,  9F)  had  gross  horizontal  projected 
areas  of  10.0,  20.0,  and  30.0  ft^.  Figure  32  shows  the  planforms  of  these 
configurations.  Figure  33  shows  three-views  of  Configuration  3GA;  this 
configuration  is  dimensionally  similar  to  the  Vought  YGBM-110  Tomahawk, 
but  incorporates  a  pop-down  inlet  similar  to  that  of  the  General  Dynamics 
YGBM-109  Tomahawk.  Configuration  3GA  is  a  baseline  configuration,  resem¬ 
bling  both,  but  not  identical  to  either  of  the  existing  Tomahawk  missiles. 
Configurations  3GA,  6GA,  9GA  have  aspect  ratios  of  9.6:  Configurations  3F, 
6F,  9F  have  aspect  ratios  of  8.35. 


GROSS  WING  AREAS 


CONFIG.  9GA 

SWING  =  34.5  FT  2 

/**  - 

_ 


CONFIG.  6GA 


SWING  =  23.0  FT2 


CONFIG.  3GA 

s  o 

WING  =  11.5  FT* 


FIGURE  32.  PLANFORMS  OF  JOINED  WING  CONFIGURATIONS  WITH 
VARIOUS  WING  AREAS 


FIGURE  33.  CONVENTIONAL  CRUISE  MISSILE  CONFIGURATION  3GA 
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An  additional  configuration  was  also  studied.  This  configuration, 
denoted  as  6H,  is  shown  in  Fig.  34.  It  is  a  hybrid  configuration  com¬ 
bining  a  joined  wing  with  a  canard. 

The  configuration  numbers  denote  very  approximately  the  non-sustained 
maneuver  capability  of  each  configuration  at  M  =  0.7,  S.L.,  maximum  weight 
conditions.  Thus,  Configurations  3GA,  3F  can  pull  approximately  3  g's, 
Configurations  6F,  6GA,  6H  can  pull  approximately  6  g's,  etc.  These 
values  are  not  necessarily  sustained  g's.  The  engine  thrust  may  not  be 
sufficient  to  overcome  the  airframe  drag  during  A  steady  maneuver  at  the 
designated  number  of  g's. 

The  initial  cruise  weight  of  all  configurations  is  2864  lb.  Thus 
the  more  maneuverable  configurations  sacrifice  some  fuel  because  of  the 
larger  weight  of  their  lifting  surfaces.  (The  increased  lifting  surface 
volume  also  reduces  fuel  tankage,  as  discussed  in  Section  5.) 

Performance  comparisons  between  the  configurations  are  given  in 
subsequent  sections  of  this  report.  The  remainder  of  this  section  dis¬ 
cusses  some  mechanical  and  configuration  design  possibilities  of  joined 
wing  cruise  missiles.  These  include  folding,  wing-warping  type  control, 
variable  camber,  and  staging. 

3.3  JOINED  WING  FOLDING  SYSTEM  FOR  CRUISE  MISSILES 


For  a  submarine-launched  cruise  missile  the  dimensions  of  the  folded 
configuration  are  constrained  so  that  it  can  fit  inside  a  standard  21" 
diameter  torpedo  tube.  For  the  present  study  it  was  assumed  that  the 
missile  is  stored  within  a  21"  diameter  casing  prior  to  launch,  as  is 
the  General  Dynamics  Tomahawk.  The  fuselage  diameter  was  therefore  set 
equal  to  that  of  the  G.D.  Tomahawk,  i.e.,  20.375  inches.  Thus  the  wings 
must  fold  into  a  tube  of  20.375  inches  diameter. 

The  design  goals  listed  in  Section  3.1  profoundly  influenced  the 
derivation  of  the  folding  system.  To  attain  the  design  goal  of  relia¬ 
ble  unfolding,  the  front  and  rear  wings  should  be  connected  throughout 
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deployment.  The  design  goal  of  symmetric  unfolding  requires  cross-links 
between  port  and  starboard  wings.  A  further  design  goal  is  that  the 
folded  configuration  should  be  robust.  For  example  it  should  be  capable 
of  being  lashed  down  for  storage  in  naval  vessels  which  may  experience 
violent  accelerations  due  to  attack  by  bombs  or  depth  charges.  This  rules 
out  folding  systems  in  which  the  folded  wings  rest  on  top  of  the  fuselage. 
Since  the  wings  must  therefore  fold  within  the  fuselage  they  must  do  so 
in  a  way  which  does  not  encroach  into  the  space  required  for  the  warhead, 
the  guidance  system,  or  other  equipment.  The  joined  wing  folding  system 
described  below  satisfies  these  requirements,  and  occupies  only  50%  of  the 
fuselage  volume  required  by  a  folded  monoplane  wing  of  the  same  span  and 
area,  employing  the  same  airfoil. 

The  proposed  folding  system  involves  three  key  concepts,  illustrated 
in  Figs.  35  through  38,  and  described  below. 

(1)  Port  and  starboard  fin  halves,  hinged  at  their  fuselage 
junctions  as  shown  in  Fig.  35,  so  that  they  can  be  rotated 
in  a  rolling  sense.  Figure  35  shows  three  alternative 
arrangements:  the  curved  fin  arrangement  occupies  the 
least  fuselage  volume,  but  is  less  stiff  and  has  higher 
drag  than  the  twin  fin  or  split  single  fin  arrangement. 

The  latter  has  the  least  drag  and  is  also  the  stiffest 
configuration. 

(2)  A  two-axis  hinge  at  the  root  of  the  rear  wing:  as  shown 

by  photos  showing  deployment  of  a  half-model  (Fig.  36)  this 
two-axis  hinge  permits  the  rear  wing  halves  to  pivot  out¬ 
wards  while  the  fin  halves  rotate  upwards. 

(3)  A  slot  in  the  lower  surface  of  the  rear  wing,  which,  for 
aerodynamic  cleanness,  is  covered  in  cruise  by  a  flush-fitting 
flap.  This  slot  guides  a  spigot  located  on  the  upper  sur¬ 
face  of  the  front  wing  tip,  to  keep  the  front  and  rear  wings 
connected  during  deployment. 

Figures  37  and  38  illustrate  other  details  of  the  folding  system. 

Both  these  figures  show  configurations  with  joined  wings  of  31.3  ft^  total 
gross  projected  horizontal  area.  This  is  larger  than  the  largest  of  the 
configurations  shown  in  Figs.  32  and  34.  We  selected  such  a  large 
wing  area  to  demonstrate  that  the  folding  system  is  suitable  even  for 
highly  maneuverable  cruise  missiles. 
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FIGURE  35.  ALTERNATIVE  FOLDING  FIN  ARRANGEMENTS 


FIGURE  36(a).  PHOTOS  OF  HALF-MODEL  OF  JOINED  WING  CRUISE  MISSILE 
SHOWING  WING  DEPLOYMENT  SEQUENCE  (CONTINUED  OVERLEAF) 


FIGURE  36(b).  PHOTOS  OF  HALF-MODEL  OF  JOINED  WING  CRUISE  MISSILE 
SHOWING  WING  DEPLOYMENT  SEQUENCE  (CONCLUDED) 
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FIGURE  37.  JOINED  WING  FOLDING  SYSTEM  FOR  HIGHLY  MANEUVERABLE  CRUISE  MISSILE 

(TWIN  FIN  CONFIGURATION) 


« 


The  fuselage  slot,  through  which  the  wings  deploy,  is  covered  in 
cruise  by  a  flush-fitting  panel  shaped  like  an  elongated  letter  L.  The 
upright  portion  of  the  L  seals  the  cut-out  in  the  fuselage  wall  through 
which  the  fin  deploys;  the  horizontal  portion  extends  as  far  forward  as 
the  trailing  edge  of  the  front  wing.  Sealing  the  junction  of  the  front 
wing  and  fuselage  can  be  done  by  similar  methods  to  those  used  on  exist¬ 
ing  conventional  missiles  with  folding  wings,  e.g.,  sliding  plates  or 
elastomeric  seals.  Note  that,  immediately  prior  to  the  wings  deploying, 
the  L-shaped  panel  rotates  downwards  through  approximately  180  degrees. 

The  vertical  arm  of  the  L  thus  provides  directional  stability  during 
deployment  until  the  fins  reach  their  cruise  position,  after  which  the 
L-shaped  panel  swings  back  up  and  closes  to  provide  a  smooth  fuselage 
contour. 

The  model  shown  in  Fig.  36  has  demonstrated  that  deployment  and 
locking  can  be  done  rapidly  and  smoothly,  with  no  tendency  for  the 
mechanism  to  stick  at  any  point  before  the  fully  deployed  position. 

Mode  of  Operation:  In  the  folded  position  the  fin  halves  are  stored 
within  the  fuselage  and  the  rear  wing  is  nested  above  the  front  wing. 

The  nested  wings  do  not  extend  across  the  full  width  of  the  fuselage. 

The  upper  and  lower  portions  of  the  fuselage  are  connected  by  a  center- 
line  structural  connection  approximately  1.5"  wide  which  acts  as  a  shear 
web  and  can  accommodate  fuel  lines  and  electrical  wiring.  By  locating 
the  front  wing  below  the  fuselage  centerline,  the  conventional  warhead 
can  be  accommodated.  The  low  front  wing  location  also  benefits  the  joined 
wing  induced  drag  and  structure  weight . 

Deployment  commences  with  simultaneous  opening  of  the  fin  and  the 
front  wing,  causing  the  folded  rear  wing  to  swing  out  approximately 
parallel  to  the  fuselage  (see  Fig.  36).  Deployment  continues  with  the 
overlapping  wings  rotating  in  a  plane  containing  the  top  of  the  fins  and 
the  front  wing  root.  This  plane  is  initially  parallel  to  the  fuselage 
waterline,  but  as  the  fins  rotate  upwards  the  tilt  of  this  plane  in¬ 
creases.  The  wings  are  linked  by  a  sliding  single  axis  bearing.  This 
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bearing  slides  in  a  slot  which  is  located  in  the  lower  surface  of  the 
rear  wing.  Spring-loaded  seals  (discussed  below)  close  these  slots 
after  passage  of  the  bearing.  Figure  37  specifies  that  the  inter-wing 
bearing  shall  be  mounted  in  an  elastomeric  block.  This  block  allows 
the  bearing  axis  to  tilt  2  or  3  degrees  in  the  spanwise  and  fore-and- 
aft  directions  to  avoid  binding  resulting  from  hinge  misalignments  or 
from  wing  flexing. 

At  the  end  of  deployment  the  two  fins  and  the  front  and  rear  wings 
are  secured  by  spring-loaded  locking  pins.  The  fuselage  wall  cut-outs 
are  also  filled  by  spring-loaded  plates  of  single  curvature,  so  that  the 
fuselage  exterior  surface  is  smooth. 

No  analyses  have  been  performed  of  the  dynamics  of  deployment,  but 
the  joined  wing  cruise  missile  is  expected  to  be  as  satisfactory  as  its 
orthodox  counterpart  provided  the  following  recommendations  are  followed. 

(1)  The  port  and  starboard  halves  of  the  front  wings  should 
be  joined  by  a  cross-linkage  to  ensure  symmetry  during 
deployment .  This  will  also  yield  approximately  symme¬ 
trical  deployment  of  the  rear  wings,  even  if  the  two 
halves  of  the  fin  unfold  at  different  rates. 

(2)  Deployment  should  be  rapid,  to  provide  centrifugal 
stiffening  of  the  outboard  section  of  the  rear  wing. 

(During  deployment  this  section  is,  in  effect,  a  swept- 
forward  cantilever  wing,  and  as  such  will  be  prone  to 
divergence.  A  helicopter  blade  encounters  a  similar 
dynamic  situation  during  one  quadrant  of  each  cycle, 
but  does  not  diverge  due  to  centrifugal  stiffening. 

Considering  the  lack  of  root  bending  constraint  and 
the  high  aspect  ratio  of  the  helicopter  blade,  we  do 
not  anticipate  that  divergence  will  occur  for  the  much 
stiffer  joined  wing.) 

(3)  The  booster  should  remain  attached  during  deployment  so 
that  thrust  vector  control  can  be  retained  until  the 
aerodynamic  controls  become  effective. 

(4)  The  flap  covering  the  undersurface  slot  on  the  outboard 
portion  of  the  rear  wing  should  be  capable  of  closing 
even  if  the  rear  wing  is  bent  by  airloads.  This  implies 
either  a  hinged  flap  with  the  hinge  and  flap  broken  up 
into  a  number  of  short  segments,  or  alternatively  a 
sliding  cover  as  sketched  in  Fig.  39. 
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ROLLER  SUPPORT 


FIGURE  39.  SLIDING  COVER  FOR  UNDERSURFACE  SLOT  IN  REAR  WING 


Comparisons  with  Conventional  Wing  Folding  Systems;  Figure  40  shows 
schematically  two  wing  folding  systems  employed  on  conventional  cruise 
missiles.  The  Vought  Tomahawk  employs  the  "scissors"  arrangement,  the 
General  Dynamics  Tomahawk  employs  "penknife"  type  folding  with  a  single 
post.  (Two-post  penknife  folding  has  been  used  on  some  other  airframes, 
e.g.,  the  Brunswick  Samson  decoy  vehicle,  but  gives  too  narrow  a  root 
chord  for  the  relatively  slender  cruise  missile.)  While  the  joined  wing 
folding  system  described  here  is  more  complex  than  the  General  Dynamics 
system,  it  requires  only  the  same  total  number  of  hinges. 


G .  D .  Tomahawk 


Joined  Win 


2  Wing  hinges 
2  Horizontal  tail  hinges 
1  Upper  fin  hinge 
1  Under  fin  hinge 


6  *  Total 


2  Front  wing  hinges 
2  (Two-axis)  rear  wing  hinges 
1  Port  fin  hinge 
1  Starboard  fin  hinge 


6  *  Total 


The  front  joined  wing  cannot  deploy  while  leaving  the  rear  wing 
folded.  To  avoid  the  corresponding  situation  (wings,  but  not  tail, 
deploying)  on  the  conventional  cruise  missile  a  special  interconnecting 
wing-tail  linkage  would  have  to  be  incorporated. 

Folded  joined  wings  occupy  less  fuselage  volume  than  comparable 
conventional  wings.  This  is  important  for  advanced  cruise  missiles  with 
high  maneuver  capability  (and  therefore  large  wings) .  Consider  a  mono¬ 
plane  wing  of  span  b,  mean  chord  c,  and  thickness/chord  ratio  t/c.  The 
wing  area  is  S  *  be  and  the  volume  is  kSc(t/c),  where  k  is  a  constant 
which  depends  on  the  airfoil  section.  Now  suppose  we  replace  this  wing 
by  a  joined  wing  pair  also  of  total  area  S,  span  b,  and  employing  the 
same  airfoil.  The  chord  of  each  joined  wing  is  therefore  half  the  orig¬ 
inal  chord.  The  ratio  of  wing  volumes  is: 


FIGURE  40.  CANTILEVER  WING  FOLDING  SYSTEMS 


i 
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Total  joined  wing  volume 
Monoplane  wing  volume 


2  x  k  x  (s/2)  x  (c/2)  x  (t/c) 
k  x  S  x  c  x  (t/c) 


Thus  the  joined  wing  occupies  only  half  the  fuselage  volume  of  the  mono¬ 
plane  wing.  More  volume  is  therefore  available  for  fuel. 

A  more  precise  calculation  than  that  given  above  would  take  into 
account  horizontal  tail  volume,  the  relative  wing  areas  required  for 
similar  "g”  capabilities,  and  the  ullage  (unused  volume*)  associated  with 
each  type  of  wing.  As  shown  in  Section  5,  these  factors  do  not  substan¬ 
tially  change  the  above  conclusion. 

Wing-folding  for  highly  maneuverable  vehicles  of  conventional  lay¬ 
out  poses  severe  mechanical  problems.  Consider,  for  example.  Configura¬ 
tion  9GA  of  Fig.  31.  The  root  wing  chord  at  the  fuselage  sides  is 
2.0735  ft  *  24.88  in.  This  exceeds  the  maximum  body  diameter  of  20.375  in 
The  tip  chord  (1.98  ft  ■  23.76  in.)  also  exceeds  the  fuselage  diameter. 

To  compound  the  problem,  the  wing  span  (18.18  ft)  exceeds  the  fuselage 
length  (17.84  ft).  Accommodating  this  wing  inside  the  fuselage,  while 
also  finding  room  for  the  tail,  may  be  exceedingly  difficult.  By  contrast 
the  joined  wing  shown  in  Figs.  36-38  has  a  gross  projected  area  (31.3  ft^) 
comparable  with  the  9GA  configuration  wing  (34.5  ft^),  yet  it  folds  leav¬ 
ing  ample  room  for  the  warhead,  guidance  system,  fuel,  etc.  If  desired, 
a  still  larger  joined  wing  could  be  employed,  by  overlapping  the  folded 
wings.  Wing  areas  of  approximately  60  ft^  could  be  accommodated  in  this 
way. 

Inlet  Flap:  The  flush  inlet  shown  on  Figs.  37  and  38  is  similar  to 
that  employed  in  the  Vought  YBGM-110  Tomahawk.  The  folding  system  does 
not  encroach  upon  this  inlet,  and  also  is  compatible  with  pop-down  inlets 
such  as  employed  in  the  G.D.  YGBM-109  Tomahawk.  Advanced  cruise  missiles 

*Airfoil  area  is  typically  0.60  to  0.65  times  the  product  of  chord  and 
maximum  thickness.  Thus  the  ullage  due  to  folding  an  airfoil  into  a  box¬ 
shaped  compartment  of  length  £  is  0.35  to  0.40  times  £c^ (t/c) . 
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may  require  larger  inlets.  This  does  not  present  an  obstacle  to  folding 
the  joined  wings.  All  that  is  required  is  to  incorporate  a  spring-loaded 
hinged  flap  in  each  side  of  the  inlet,  as  shown  in  Fig.  38.  The  flap 
hinges  inwards  when  the  wings  are  folded,  permitting  the  wings  to  occupy 
the  empty  space  inside  the  inlet.  For  inlets  with  double  curvature 
several  short  hinged  flap  segments  should  be  employed. 

Wing  Taper:  Untapered  wings  are  more  convenient  to  fold  than  tapered 
wings  and  are  cheaper  to  manufacture.  These  advantages  are  offset  by 
structural  weight  penalties.  Hence,  for  joined  wing  configurations  having 
relatively  small  wing  area,  untapered  wings  may  be  acceptable.  (A  con¬ 
figuration  with  untapered  wings  is  described  in  Section  3.6.  This  con¬ 
figuration  reduces  wing  area  through  employing  a  canard  to  generate  part 
of  the  lift.)  For  initial  studies  a  moderate  taper  ratio  (~  0.5)  is  recom 
mended . 

3.4  WING  WARPING  CONTROL 


Since  the  front  wing  of  a  joined  wing  pair  has  a  constraint  at  its 
tip  due  to  the  tip  joint,  it  can  be  warped  for  pitch  and  roll  control  via 
actuators  mounted  near  its  root.  Such  wing-warping  control  eliminates 
hinged  control  surfaces  thus  simplifying  the  wing  structure  and  the  con¬ 
trol  system,  and  also  reduces  radar  cross-section  by  eliminating  radar 
reflections  from  control  surface  hinge  discontinuities.  The  feasibility 
of  wing-warping  has  been  studied  in  a  preliminary  fashion:  the  results 
appear  promising,  as  described  below. 


A  structural  analysis  of  the  torsional/flexural  interactions  of  the 
front  and  aft  wings  was  performed.  This  was  aimed  at  predicting  the  span- 
wise  incidence  variation  induced  by  torques  applied  at  the  front  wing  root 
A  number  of  simplifying  assumptions  were  made,  e.g.,  untapered  wings  with 
equal  chords  and  equal  (positive  and  negative)  sweep  angles,  isotropic 
material  properties,  and  infinitely  rigid  tip  joints  were  assumed.  With 
these  assumptions  the  spanwise  incidence  variation  was  found  to  be  para¬ 
bolic  across  both  wings.  For  example,  if  the  front  wing  root  is  twisted 
such  that  its  incidence  is  increased  by  4  degrees,  the  corresponding 
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induced  incidence  at  the  tip  is  1  degree,  while  at  the  mid  semi-spans 
the  incidences  are  2.25  degrees  and  0.25  degrees  on  the  front  and  rear 
wings  respectively.  This  parabolic  variation  is  typical  of  untapered 
wings  made  of  isotropic  material.  The  torsion  applied  to  the  front  wing 
is  opposed  by  the  flexural  stiffness  of  the  aft  wing,  so  by  judicious 
tailoring  of  fiber  orientation  other  (nonparabolic)  variations  can  be 
obtained,  if  desired. 

Two  prime  aspects  have  been  studied:  (a)  the  aerodynamic  forces 
and  moments  generated  by  warping;  (b)  mechanical  design  of  the  warping 
mechanism  and  integration  with  the  wing  folding  system.  Initial  indica¬ 
tions  regarding  (a)  are  that  the  steady  g's  obtained  per  degree  of  root 
twist  are  typically  1/3  to  2/3  of  the  g's  per  degree  of  aft  tail  deflec¬ 
tion  for  the  aft  tail  cruise  missile  having  the  same  total  lifting  sur¬ 
face  area.  For  the  Vought  Tomahawk  the  aft  tail  pitch  control  deflection 
was  limited  to  ±5  degrees.  It  appears  reasonable  to  allow  a  larger  limit 
for  wing  root  twist  (e.g.,  15  degrees)  hence  adequate  longitudinal  control 
is  expected  to  be  available  from  wing  warping  provided  premature  flow 
separation  does  not  occur  near  the  root  of  the  front  wing.  To  ensure  this 
the  root  airfoil  must  be  properly  designed  and  a  good  wing  root  seal  must 
be  maintained.  The  importance  of  proper  sealing  of  the  wing  root  for  main¬ 
taining  good  fuselage  "carry-over"  of  wing  lift  has  been  demonstrated  by 
wind  tunnel  test  data  on  a  folding  missile  given  in  Ref.  34. 

The  wing  warping  mechanism  must  satisfy  some  difficult  spatial  con¬ 
straints  since  it  should  be  located  within  the  fuselage  below  the  front 
wing  root  in  order  to  leave  the  space  above  the  root  free  for  warhead  or 
fuel.  In  addition,  it  is  desirable  that  the  actuators  and  gears  remain 
fixed  during  deployment.  In  the  scheme  shown  in  Fig.  41,  the  front  wing 
is  mounted  in  a  yoke  which  has  a  transverse  bearing  to  permit  wing  deploy¬ 
ment.  The  yoke  can  be  rotated  about  an  axis  parallel  to  the  wing  quarter 
chord  line  by  means  of  a  spur  gear  segment  driven  by  a  reversible  servo 
motor.  The  conceptual  design  shown  in  Fig.  41  complies  with  the  spatial 
constraints  listed  above.  It  is  also  mechanically  and  structurally 


FIGURE  41.  WING-WARPING  CONTROL  SYSTEM 


simple  and  should  have  less  backlash  than  conventional  control  systems, 
which  use  push-rods  to  actuate  hinged  flaps,  elevators,  and  ailerons. 

3.5  VARIABLE  CAMBER 


Variable  camber  has  long  been  recognized  as  yielding  significant 
aerodynamic  advantages  for  aircraft  that  must  develop  high  L/D  over  a 
wide  range  of  lift  coefficients.  Unfortunately,  a  conventional  canti¬ 
lever  wing  presents  considerable  mechanical  difficulties  in  mechanizing 
variable  camber.  This  is  because  numerous  actuators  must  be  mounted 
within  the  wings.  For  a  joined  wing  the  possibility  exists  of  varying 
camber  via  two  sets  of  actuators,  one  set  located  within  the  fuselage  at 
the  wing  root,  the  other  set  located  inside  the  tip  joint  fairing.  Such 
a  system  would  not  induce  constant  camber  across  the  entire  span,  but 
conceivably  the  spanwise  variation  of  camber  might  be  small,  and  accept¬ 
able  in  view  of  the  saving  in  the  numbers  of  actuators  required.  Addi¬ 
tionally,  since  the  actuators  are  outside  the  airfoils,  thin  airfoils 
could  be  employed  (e.g.,  for  supersonic  cruise). 

For  a  cruise  missile  the  prime  benefit  of  variable  camber  is  that, 
as  fuel  is  consumed,  and  wing  loading  reduced,  camber  can  be  reduced  to 
maintain  the  highest  possible  L/D,  thus  maximizing  range.  This  possi¬ 
bility  was  considered  in  the  present  study  but  was  not  pursued  because: 

(1)  To  demonstrate  the  feasibility  of  variable  camber  would 
have  required  more  extensive  mechanical  design  studies 
than  could  be  justified  for  the  funds  available. 

(2)  An  attractive  alternative  form  of  variable  geometry 
exists,  which  is  well  suited  to  the  joined  wing;  this 
is  "staging,"  which  is  described  below. 

3.6  STAGED  JOINED  WING  CONFIGURATION 


Figure  42  shows  a  configuration  employing  joined  wings  which  deploy 
similarly  to  those  of  Fig.  36.  The  configuration  also  includes  a  canard 
surface  which  is  attached  to  a  nose  tank  which  is  ejected  when  empty. 
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FIGURE  42.  HYBRID  TWO-STAGE  CONFIGURATION  WITH  UNTAPERED  JOINED  WINGS 
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Such  a  configuration  is  feasible  with  the  downward -looking  guidance  and 
terrain-following  system  (TERCOM)  used  in  current  Tomahawk  cruise  missiles 
and  projected  cruise  missiles.*  The  objectives  of  this  staged  configura¬ 
tion  are: 

(1)  To  improve  range,  speed  and  maneuverability  by  reducing 
missile  drag  and  weight  during  the  latter  portion  of 
the  mission. 

(2)  To  reduce  radar  cross-section. 

(3)  To  reduce  the  required  operating  range  of  on  the 

joined  wings,  permitting  the  use  of  airfoils  of  lower 
drag  coefficient. 

(4)  To  reduce  the  fuselage  volume  taken  up  by  the  folded 
lifting  surfaces  (for  a  given  thickness/chord  ratio 
and  given  total  area  three  lifting  surfaces  occupy  less 
volume  than  two,  or  one). 

(5)  To  provide  a  modular  missile  in  which  the  lifting  sur¬ 
faces  do  not  intrude  on  the  fuselage  sections  occupied 
by  the  guidance  equipment  and  warhead.  Figure  43  shows 
an  inboard  profile  which  illustrates  this  advantage. 

Such  a  modular  arrangement  facilitates  the  use  of  alter¬ 
native  warheads  and  guidance  systems. 

(6)  To  provide  a  decoy  capability.  If  desired,  the  nose 
section  can  be  powered  (e.g.,  by  a  booster)  so  that,  as 
the  missile  closes  with  its  target,  the  nose  section 
flies  off,  thus  providing  an  additional  radar  signal 

to  confuse  the  enemy  during  a  critical  phase  of  the 
attack.  By  using  the  canard  plus  pop-out  fins,  the 
decoy  can  maneuver . 

Figure  44  illustrates  the  proposed  method  of  fairing  the  nose  tank. 
This  leaves  a  smooth  fuselage  contour  after  the  nose  has  separated.  Ex¬ 
perimental  data  given  on  p.  5-7  of  Ref.  15  indicates  that  the  drag  penalty 
for  the  1/32"  step  surface  irregularity  that  exists  before  separation 
should  be  small,  representing  only  a  1.2%  increase  over  the  drag  of  a 
smooth  fuselage.  For  conservatism,  the  performance  calculations  of 
Section  5  assume  a  2%  drag  increase,  corresponding  to  a  1/20"  step. 

*Down ward-loo king  radar  altimeter  systems  are  preferred  to  forward- 
looking  systems  because  they  give  a  smaller  radar  signature  to  the  enemy. 


FIGURE  A3.  INBOARD  PROFILE  OF  HYBRID  CONFIGURATION  OF  FIG.  A2 
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FIGURE  44.  DETAIL  OF  NOSE  TANK  ATTACHflENT  GEOMETRY  FOR  TWO-STAGE  CONFIGURATION 


For  the  configuration  shown  in  Fig.  42,  ejecting  the  nose  tank 
reduces  the  missile's  total  wetted  area  by  19%  and  empty  weight  by  5%. 
Detailed  estimates  of  the  performance  benefits  of  staging  are  given  in 
Section  5. 

A  key  factor  determining  the  practicality  of  staging  is  the  longi¬ 
tudinal  stability  of  the  configuration.  This  is  best  measured  by  the 
maneuver  margin,  but  for  initial  studies  the  static  margin  provides  an 
acceptable  index  of  stability.  Figures  45  and  46  illustrate  the  effect 
of  fuel  consumption  and  staging  on  static  margin  for  the  configuration 
of  Fig.  42.  Figure  45  applies  for  the  canard  surface  having  a  fixed  sweep 
angle,  Ac/4  "  35°  as  drawn  in  Fig.  42.  Figure  46  applies  for  a  canard 
having  a  sweep  angle  which  progressively  increases  from  35°  as  the  nose 
tank  is  emptied. 

With  constant  canard  sweep  the  static  margin  becomes  unstable  as  the 
nose  tank  is  emptied,  as  shown  by  the  upper  lines  in  Fig.  45.  Immediately 
after  ejecting  the  nose  tank  the  static  margin  becomes  slightly  unstable, 
but  it  becomes  much  more  stable  as  the  aft  fuel  is  consumed. 

Although  a  slightly  unstable  static  margin  is  tolerable,  it  would 
be  preferable  to  have  a  small  positive  static  margin  throughout  the  flight. 
Figure  46  shows  how  this  can  be  obtained  by  varying  the  canard  sweep  as 
nose  fuel  is  burned.  The  lower  part  of  Fig.  46  illustrates  how  exces¬ 
sively  large  static  margin  can  be  delayed  via  fuel  management,  i.e. , 
splitting  the  aft  fuel  into  two  tanks  and  emptying  the  more  forward  tank 
first.  A  similar  benefit  can  be  obtained  with  a  single  aft  tank  by  trim¬ 
ming  the  missile  to  fly  more  nose  up  as  fuel  is  consumed.  (This  can  be 
done  via  the  direct  lift  capability  of  the  joined  wing.)  Neither  fuel 
management  nor  attitude  adjustment  can  change  the  static  margin  when  the 
fuel  tanks  are  nearly  empty.  However  this  limitation  is  not  expected  to 
be  serious  because: 
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GROSS  WEIGHT 


FIGURE  45.  STATIC  MARGIN  VARIATION  WITH  WEIGHT: 
FIXED  CANARD  SWEEP  ANGLE 


GROSS  WEIGHT 


« 


FIGURE  46.  STATIC  MARGIN  VARIATION  WITH  WEIGHT: 
VARIABLE  CANARD  SWEEP  ANGLE 
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(1)  Only  a  small  percentage  of  missions  will  require 
maximum  range. 


(2)  Near  maximum  range  the  low  gross  weight  will 
improve  maneuverability. 

(3)  With  direct  lift  control,  excessive  static  margin 
does  not  impair  maneuverability  as  much  as  it 
does  for  aft  tail  control. 

From  the  above  discussion,  staging  appears  to  be  practical  and  promis¬ 
ing.  The  joined  wing  configuration  is  particularly  well  suited  to  staging 
because  of  the  inherently  wide  c.g.  range  over  which  it  can  be  trimmed. 

(If  desired,  100%  of  the  lift  can  be  carried  on  the  front  wing.) 

Staging  is  not  suitable  for  aft  tail  configurations  because  the  c.g. 
range  shown  on  Figs.  45  and  46  (approximately  25")  is  large  compared  to 
the  normal  trim  range  (with  conventional  tail/wing  area  ratios).  It  does 
not  appear  feasible  to  alleviate  this  situation  by  variable  wing  sweep, 
e.g.,  going  from  sweepforward  to  sweepback  as  the  c.g.  moves  aft.  This 
demands  excessive  sweep  angle  variations  and/or  spans. 

The  next  section  of  this  report  discusses  terrain-following.  Readers 
who  are  interested  in  performance  rather  than  guidance  and  control  may 
prefer  to  read  Section  5  before  Section  4. 
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4.0  AIRFRAME  DESIGN  FOR  TERRAIN-FOLLOWING 


4.1  APPLICATION  OF  OPTIMAL  CONTROL  THEORY 

A  cruise  missile,  like  a  submarine,  depends  critically  on  stealth. 

This  is  achieved  by: 

(1)  Designing  the  vehicle  to  have  low  observables. 

(2)  Flying  close  to  the  surface  of  the  sea  or  terrain. 

In  this  section  we  investigate  the  consequences  of  Requirement  (2)  on 
the  best  selection  of  missile  airframes. 

If  very  low  altitude  flight  is  to  be  performed  successfully  the 
flight  path  of  the  missile  must  reflect  the  irregularities  of  the  terrain 
or  waves.  Terrain-following  with  low  clearance  requires  an  effective 
guidance  system.  However  as  will  be  shown,  even  with  an  optimum  guidance 
system,  terrain-following  accuracy  is  limited  by  the  dynamic  response 
characteristics  of  the  airframe.  Hence,  comparing  the  terrain-following 
accuracy  of  joined  wing  and  conventional  airframes,  each  having  indi¬ 
vidually  optimized  guidance  systems,  provides  one  facet  of  an  overall 
comparison  of  airframes.  To  perform  such  a  comparison,  criteria  for 
optimality  must  first  be  defined.  This  can  be  done  as  follows. 

Current  and  projected  cruise  missiles  employ  downward -loo king  terrain¬ 
following  systems  which  measure  the  missile's  clearance  above  the  terrain, 
i.e.,  hc  in  Fig.  47.  (Downward,  rather  than  forward-looking  guidance  is 
employed  because  it  is  less  easily  detected  by  the  enemy.)  A  bias  set¬ 
ting  is  incorporated  in  the  guidance  logic,  corresponding  to  a  desired 
constant  clearance  hg.  The  difference  between  hg  and  hj.,  at  any  instant, 
is  the  height  error,  he.  If  he  has  a  symmetrical  statistical  distribu¬ 
tion,  the  mean  value  of  he,  E(he),  may  be  assumed  to  be  zero.  To 
achieve  stealth  hg  should  be  low  and  upward  deviations  from  hg  should  be 
small.  With  a  low  hg,  downward  deviations  from  hg  must  also  be  small 


FIGURE  47.  CRUISE  MISSILE  FLYING  ABOVE  TERRAIN 
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to  reduce  the  probability  of  hitting  the  terrain  (usually  called  "proba¬ 
bility  of  clobber,"  Pc) •  Therefore  a  reasonable  optimality  criterion 
for  the  guidance  system  is  that  it  should  minimize  the  mean  square  height 
error,  E(h|).  This  criterion  has  the  advantage  of  being  mathematically 
convenient,  and  it  is  widely  employed  in  Optimal  Control  Theory.  This 
theory  shows  how  to  compute  the  guidance  system  that  (for  a  given  air- 

frame)  yields  minimum  E(h|) .  The  theory  also  yields  short-cut  methods 

2 

for  computing  the  minimum  E(he)  and  also  the  associated  mean  square  con¬ 
trol  deflection  E(6^) . 

2 

While  there  are  advantages  to  employing  minimum  E(hg)  as  a  criterion 

o 

for  terrain  following,  it  should  be  realized  that  the  minimum  E(h|)  system 
is  only  an  approximation  to  the  minimum  pc  system.  This  is  because  pc 
depends  on  the  frequency  characteristics  of  he  as  well  as  its  statistical 
amplitude  distribution.  For  example,  consider  two  Gaussian  he  time  his¬ 
tories,  having  identical  E(h|)  but  with  the  first  having  a  higher  band¬ 
width  than  the  second.  For  equal  flight  durations,  the  high-bandwidth  he 
will  have  the  greater  number  of  exceedances  of  any  specified  boundary 
value.  (For  further  discussions  of  the  relative  utility  of  E(h|)  and  pc 
as  optimization  criteria  see  Ref.  36.  References  37  and  38  are  also  of 
interest:  these  explain  some  of  the  difficulties  and  high  computer  time 

requirements  involved  in  computing  pc>)  Although  E(he)  is  a  less  precise 
measure  of  system  merit  than  pc,  it  provides  a  valid  and  cost-effective 
basis  for  comparison  of  alternative  preliminary  configuration  designs  and 
will  be  used  here  for  that  purpose. 

In  the  remainder  of  this  section  we  present  power  spectra  for  terrain 
to  be  followed  (actually  ocean  waves)  and  also  for  atmospheric  turbulence 
appropriate  to  low  altitudes.  Transfer  functions  are  then  given  for  the 
alternative  joined  wing,  hybrid  (i.e.,  joined  wing  +  canard),  and  conven¬ 
tional  airframes  shown  in  Figs.  31  through  34.  The  best  achievable 
terra in -folio wing  performance  of  each  configuration  is  then  calculated 
(i.e.,  the  minimum  achievable  E(hg) .  It  is  shown  that  the  hybrid  and 
the  joined  wing  configurations  are  potentially  capable  of  closer -terra in 
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following  than  the  conventional  aft  tail  configurations.  This  advantage 
is  due  to  the  fact  that  the  aft  tail  airframes  initially  move  in  the 
wrong  direction  following  elevator  commands  intended  to  change  altitude. 
By  contrast,  the  joined  wing  configurations  have  responses  which  are 
always  in  the  desired  direction.  Finally,  it  is  shown  that,  for  the 
wind  and  wave  conditions  specified,  practical  realization  of  this  advan¬ 
tage  requires  height  to  be  controlled  by  two  control  surfaces.  Thus  the 
hybrid  configuration  requires  simultaneous  control  of  front  wing  flaps 
and  canard  incidence  for  closest  terrain-following.  The  basic  joined 
wing  configuration  requires  the  lift  generated  on  each  wing  to  be  con¬ 
trolled,  e.g.,  by  flaps  on  each  wing. 

4.2  SYSTEM  DESCRIPTION 


Under  this  heading  we  describe  the  mathematical  models  employed  to 
represent  the  airframe,  its  atmospheric  turbulence  environment,  and  the 
"terrain"  (sea  surface).  A  more  complete  description  of  each  of  these 
is  given  in  Appendix  A:  only  the  bare  essentials  are  presented  here. 

Airframe  Equations  of  Motion:  Short-period  longitudinal  equations 
are  employed;  it  is  assumed  that  airspeed  variations  are  suppressed  per¬ 
fectly  by  an  outer  throttle  loop  not  shown  on  Fig.  48.  The  equations  of 
motion  are  transformed  from  the  standard  form  given  in  Ref.  39:  they 
are  written  in  terms  of  height  perturbations  (h)  and  pitch  angle  pertur¬ 
bations  (9)  as: 


S  -  Zw  -Zq 

-sh 

V 

6  + 

+  ZqS/UQ 

_-M^s-Mw  s2  -  (Mq+MQ)s -Ma_ 

0 

m6 

.-Mw  +  MqS/Uo-MwS_ 

6  denotes  a  single  control  surface  deflection,  subsequently  multiple 
controls  are  denoted  as  6^,  62.  Wg  is  the  vertical  component  of  turbu¬ 
lence.  The  terms  Zqs/U0,  Mqs/U0  multiplying  Wg  are  "gust-gradient" 
terms  (Ref.  39,  p.  251).  Following  customary  practice  the  Zq  term  is 
neglected,  since  Zq  «  U0. 
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FIGURE  48.  TERRAIN-FOLLOWING  SYSTEM  BLOCK  DIAGRAM 


The  principal  transfer  functions  of  interest  are  those  describing 
the  height  response  to  control  and  also  to  vertical  gust  components. 


-Zgs2  +  Zg(Mq  +  Mq)s  +  Z6M6  -  MgZa 
s2[s2-  (Mq  +  Mq  +  Zw)s  -  +  Mqzw3 


(11) 


h  _  _ 2w(s  -  2Mq) _ 

wg  s2[s2-  (Mq  +  +  Zw)s  -  ^  +  MqZw] 


(12) 


For  brevity,  the  denominator  of  these  transfer  functions  will  be  written 
s2(s2  +  Bs  +  C) .  Initially  "single-controller"  systems  are  considered, 
so  6  is  a  scalar  quantity.  "Multi-controller"  systems  involving  con¬ 
trols  6j,  62,  ....  will  be  discussed  in  Section  4.4. 

Sea  Spectrum:  As  explained  in  Appendix  A,  a  modified  Bretschneider 
spectrum  is  employed,  appropriate  to  a  short -crested  sea,  of  Sea  State  5, 
and  a  20  knot  windspeed.  For  a  missile  traversing  this  sea  at  M  =  0.7, 
flying  into  the  inertial  wind  direction,  the  wave  height  encountered 
frequency  characteristics  are  described  by  the  following  power  spectrum: 


^hw^w 


818.9  u)2 
(0)2  +  6.32)2 


(13) 


where 

»  Wave  height  above  mean  sea  level,  ft 
o)e  =  Frequency  of  encounter,  rad/ sec 

Figure  49  graphs  this  spectrum.  Note  that  it  is  of  relatively  narrow 
bandwidth,  with  most  of  the  power  concentrated  in  the  region  2.0 
20.0  rad/sec.  This  is  different  from  typical  land  terrain  spectra 
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FIGURE  49.  COMMAND  AND  DISTURBANCE  SPECTRA 
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(Refs.  37  and  38)  which  do  not  attenuate  at  low  frequencies.  The  rms 
wave  height  is  9.18  ft. 


Turbulence  Spectrum:  The  well-known  Dryden  spectrum  is  employed. 
As  shown  in  Appendix  A,  for  the  assumed  flight  condition  this  reduces 
to : 


4vgwg 


(14. 0438) 2  (cd2  +  4.5143)2 
(a)2  +  7.8192)2 


(14) 


The  rms  value  of  wg  is  2.9  fps.  This  spectrum  is  also  shown  on  Fig.  49. 

Stability  Derivatives:  Stability  derivatives  were  calculated  for 
each  of  the  configurations  3GA,  6GA,  9GA,  3F,  6F,  9F,  and  6H  shown  in 
Figs.  31  through  34.  In  all  cases  a  positive  static  margin  was  assumed 
as  listed  on  Table  1,  which  also  gives  weights  and  inertias.  The  vortex- 
lattice  computer  program  of  Ref.  19  was  employed  to  calculate  the  lift¬ 
ing  surface  contribution  to  the  joined  wing  derivatives,  with  fuselage 
corrections  from  standard  texts  (e.g..  Ref.  23).  For  the  conventional 
configuration  3GA,  full-scale  M  *  0.7  A.E.D.C.  Tullahoma  wind  tunnel 
data  on  the  closely  similar  Vought  YGBM-109  Tomahawk  cruise  missile  were 
employed  (Ref.  45):  these  were  scaled  to  Configurations  6GA  and  9GA  by 
standard  methods  (Ref.  23).  The  hybrid  configuration  planform  involves 
three  lifting  surfaces  and  therefore  cannot  be  modeled  by  the  vortex- 
lattice  program  of  Ref.  19  which  is  limited  to  two  planforms.  However, 
a  later  version  of  this  program,  originally  developed  by  J.  S.  Luckring 
for  oblique-wing  planforms,  avoids  this  restriction  and  was  utilized  to 
compute  derivatives  for  the  hybrid  configuration.  Table  2  lists  the 
derivatives  and  the  h/6  and  h/wg  transfer  functions  in  the  forms: 

JL  „  Ah5s  +  Bhgs  +  Ch^ 

<5  s2(s2  +  Bs  +  C) 


_h_ 

wg 


Av,  s  +  Bh., 
“wg 

s2(s2  +  Bs  +  C) 
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TABLE  2.  DERIVATIVES  AND  TRANSFER  FUNCTIONS 


H  WC  -1.00172160  -3.31878960  -6.15518400  -0.37921080  -1.94661954  -5.08876050  -1.54397600 


Table  1  indicates  that  although  the  static  margins  of  all  the  con¬ 
figurations  studied  are  positive  they  are  not  equal.  This  variation 
was  expected  to  have  little  effect  on  terrain-following  performance  of 
single-controller  systems:  this  expectation  was  verified  by  performing 
the  calculations  of  terrain-following  accuracy  twice,  first  with  the 
static  margin  indicated  and  secondly  with  zero  static  margin.  Only 
slight  differences  in  terrain-following  accuracy  were  found.  The  reason 
for  this  insensitivity  is  that  varying  static  margin  has  little  effect 
on  the  numerator  coefficients  of  the  h/wg  and  h/6  transfer  functions, 
and  affects  the  denominators  of  these  transfer  functions  equally,  by 
varying  the  C  coefficient.  [Varying  static  margin  can  have  significant 
effects,  however,  if  a  multicontroller  (e.g.,  f laps-plus-elevator)  height 
control  system  is  employed.  This  is  discussed  in  Section  4.4.] 

4.3  LIMITS  ON  OPTIMUM  TERRAIN-FOLLOWING  PERFORMANCE 

Tor  a  given  airframe,  the  guidance  system  can  be  designed  to  mini¬ 
mize  E(h|)  through  the  use  of  Optimal  Control  Theory.  This  yields  the 
guidance  system  that  gives  the  minimum  achievable  E(he)  for  a  specified 
level  of  mean  square  elevator  deflection  E(6^) .  (This  constraint  may  be 
generalized  for  multicontroller  systems,  e.g.,  elevator-plus-flap  con¬ 
trol,  and  additional  constraints,  e.g.,  g-limits  may  be  introduced.) 

2 

It  might  be  thought  that  any  desired  level  of  E(hg) ,  even  zero, 
could  be  obtained  by  allowing  E(6^)  to  increase.  This  is  not  true  for 
conventional  (aft  tail)  configurations  in  which  elevator  is  employed  to 
control  height.  Such  configurations  have  height -to-elevator  transfer 
functions  which  incorporate  one  or  more  right -half-plane  zeros.  (Trans¬ 
fer  functions  of  this  kind  are  called  "nonminimum  phase.")  It  has  been 
shown  (Ref.  40)  that  for  a  system  having  a  nonminimum  phase  transfer 
function  the  error  cannot  be  reduced  to  zero  even  if  infinitely  large 
mean  square  control  deflections  are  permitted.  For  such  a  system  there 
is  a  finite  minimum  achievable  mean  square  error. 

From  the  Laplace  transform  initial  and  final  value  theorems  it  can 
shown  that  the  presence  of  a  right-half-plane  transfer  function  zero  of 


the  height/elevator  transfer  function  implies  that  the  initial  direction 
of  the  airframe's  height -to-elevator  response  is  in  the  opposite  direc¬ 
tion  to  the  ultimate  response.  This  occurs  because,  in  order  to  climb, 
a  download  must  be  applied  to  the  elevator.  This  download  causes  the 
initial  response  of  the  airframe  c.g.  to  be  downwards. 

By  contrast,  for  a  canard  or  joined  wing  configuration  this  reversal 
of  response  can  be  avoided,  by  placing  the  elevator  control  on  the  front 
wing,  ahead  of  the  c.g.  Thus,  potentially,  such  "minimum  phase"  configu¬ 
rations  can  follow  given  terrain  more  closely  than  conventional  (aft 
tail)  configurations. 

The  foregoing  discussion  has  omitted  many  details  such  as  the  effects 
of  digital  processing  and  actuator  lags,  possibilities  of  multicontroller 
systems  (e.g.,  elevators-plus-f laps) ,  forward-looking  guidance,  airframe 
flexibility,  et.  These  will  be  discussed  in  Section  4.4.  First,  how¬ 
ever,  we  present  the  results  of  some  optimal  control  calculations  for 
downward-looking  single-controller  guidance  systems,  without  lags,  applied 
to  a  variety  of  joined  wing  and  conventional  cruise  missile  configurations. 

Optimum  Single -Controller  System:  Appendix  A  derives  the  equations 
for  the  optimal  single-controller  system,  i.e.,  the  system  that  minimizes 
the  performance  index,  J. 

J  *  E(h2)  +  k2E(62)  (17) 

Our  principal  interest  here  lies  not  in  the  details  of  the  optimal  com¬ 
pensation  in  Fig.  48,  but  rather  in  the  question:  "How  good  is  the 
optimum  system,"  i.e.,  how  closely  does  it  follow  the  terrain.  Thus 
for  a  given  kz  we  wish  to  know  the  value  of  E(he)  obtained,  and  the  asso- 
elated  mean  square  control  deflection  E(60 .  As  shown  in  Appendix  A, 
these  values  are  given  by  expressions  of  the  type  : 
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In  each  expression  the  first  term  gives  the  contribution  to  J  that  stems 
from  the  portion  of  he  or  6  that  is  correlated  with  the  terrain.  The 
second  term  gives  the  contribution  due  to  the  portion  of  he  or  <5  that 
is  correlated  with  the  atmospheric  turbulence.  Only  two  terms  appear 
in  each  expression:  there  are  no  cross-correlation  terms  involving 
^Rd>  ^dR  because  it  has  been  assumed  that  the  terrain  (waves)  and  tur¬ 
bulence  are  uncorrelated, 

E(he)Qpt  and  the  associated  E(<5  )  were  calculated  by  the  computer 
program  described  in  Appendices  A  and  C,  for  each  of  the  configurations 
in  Figs.  31  through  34.  The  results  are  presented  below. 

Conventional  (Aft  Tail)  Configurations:  Figure  50  shows  the  mean 
square  height  error,  E(h|) ,  for  aft  tail  configurations  3GA,  6GA,  and 

2 

9GA  graphed  versus  the  corresponding  mean  square  control  deflection  E(<$  ) 
required  by  the  optimal  control  system.  Values  of  the  weighting  parame¬ 
ter  k  are  indicated  on  the  graphs.  (The  control  deflection,  6,  indicates 
deflection  of  the  all-moving  tail  surfaces  in  units  of  degrees.)  At  high 
values  of  k  (i.e.,  "expensive  control")  corresponding  to  the  left  side 
of  the  graphs  the  control  deflections  are  small  and  have  little  effect 
on  the  missile,  which  flies  more  or  less  in  a  straight  line  with  the 
preponderant  part  of  the  height  error  being  due  to  the  waves.  Thus  the 
asymptote  of  E(h|)  for  k  -*•  00  is  essentially  equal  to  the  mean  square 

n 

wave  height,  81.2  ft  .  As  k  is  decreased  the  system  accuracy  increases, 
but  for  reasonable  values  of  E(6^),  e.g.,  of  the  order  of  10  deg^,  the 
gain  in  accuracy  is  slight.  For  example,  Configuration  3GA  reduces  its 
mean  square  error  only  to  74/81.2  ■  0.91  of  its  open-loop  value.  If 
much  larger  E(6^)  could  be  produced  or  if  the  control  effectiveness 
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parameters  Cl^,  could  be  increased  the  accuracy  could  be  improved, 

but  even  with  infinite  mean’ square  control  deflections  the  best  accuracy 
that  can  be  achieved  with  airframe  3GA  is  only  E(hg)  =  22  ft^,  i.e., 

0.27  of  the  open-loop  mean  square  error.* 

This  situation  is  only  slightly  improved  by  selecting  a  larger  wing 
area,  as  shown  by  the  graphs  for  Configurations  6GA  and  9GA  in  Fig.  50. 
Nor  is  any  significant  improvement  gained  by  reducing  static  margin. 

This  expedient  was  studied  by  recalculating  the  optimum  system  perfor¬ 
mance  with  Mq  set  equal  to  zero,  and  the  remaining  derivatives  unchanged. 

9  9 

For  any  given  k  the  zero  static  margin  values  of  E(hg)  and  E(6  )  were, 
when  graphed  on  Fig.  50,  virtually  coincident  with  the  points  graphed  for 
the  normal  static  margins  listed  in  Table  1. 

It  is  clear  from  the  above  results  that  even  with  ideal  guidance 
system  dynamics  (e.g.,  no  actuator  lags,  no  digital  processing  delays) 
the  gain  in  accuracy  resulting  from  terrain  following  is  meager  for  the 
aft  tail  configurations.  The  root  of  the  problem  can  be  traced  to  the 
nonminimum  phase  characteristics  inherent  in  any  airframe  controlled  by 
an  aft  tail.  To  demonstrate  this,  the  magnitude  of  the  nonminimum  phase 
zero  in  the  h/ 6  transfer  function  was  varied  parametrically,  and  the 
asymptotic  value  of  E(h|)  was  recalculated.  That  is,  the  h/6  transfer 
function  (Eq.  15)  was  expressed  as: 


Ah^Cs  +  a) (s  -  b) 
s^(s^  +  Bs  +  C) 


b  >  0 


and  b  was  varied  parametrically.  The  resulting  E(h^)  was  calculated  from 
Eqs.  B34  and  B35  of  Appendix  B.  Figure  51  illustrates  the  result  of  para¬ 
metrically  varying  b.  It  shows  that  for  1.0  >  b  >  100.0  the  asymptotic 


*The  asymptotes  shown  in  Fig.  50  have  been  verified  by  an  independent 
method  given  in  Appendix  B. 


FIGURE  51.  MINIMUM  ACHIEVABLE  MEAN  SQUARE  ERROR  FOR  AFT  TAIL  CONFIGURATIONS 
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error  is  very  small.  Unfortunately,  each  one  of  the  configurations  3GA, 
6GA,  9GA  has  values  of  b  which  fall  near  the  peak  of  the  graph  in  Fig.  51. 

There  is  no  simple  remedy  for  this  limitation  of  aft  tail  configu¬ 
rations.  Varying  static  margin  has  only  a  slight  effect  on  b;  even 
reducing  the  static  margin  to  zero  induces  only  a  slight  change  in  the 
value  of  b,  which  is  primarily  determined  by  the  tail  moment  arm.  To 
effectively  vary  b  requires  feedback  to  a  second  controller,  i.e.,  the 
guidance  system  should  simultaneously  control  tail  incidence  plus  one 
other  control,  e.g.,  flaps.  This  is  discussed  in  Section  A.  A  where  it 
is  shown  to  be  of  only  limited  effectiveness. 

Effect  of  Varying  Terrain;  Since  the  above  results  have  been  derived 
for  a  particular  sea  spectrum  it  is  necessary  to  investigate  the  effects 
of  deviating  from  this  spectrum.  For  example,  a  different  command  spec¬ 
trum  would  be  appropriate  over  land.  The  effects  of  varying  the  command 
spectrum  can  be  found  as  explained  in  Appendix  B,  Eqs.  B3A  through  B38 . 
These  equations  can  be  applied  to  determine  the  asymptotic  mean  square 
error  of  a  nonminimum  phase  missile  flying  over  land  terrain,  as  follows. 

Various  forms  of  land  terrain  spectra  are  in  use  (see  Refs.  37,  38, 
and  Al)  but  the  simplest  and  most  convenient  form  for  our  purpose  is  a 
first-order  Butterworth  form,  thus: 


/2a[E(R2)] 
s  +  a 


(21) 


In  this  equation  E(R  )  is  the  mean  square  terrain  height  variation,  and 
a  is  the  encounter  break  frequency,  which  depends  on  the  terrain  rough¬ 
ness  and  the  missile's  speed.  From  Eq.  B3A,  the  component  of  mean  square 
height  error  due  to  terrain-following  has  a  minimum  value  for  zero  con¬ 
trol  weighting  (k  0)  which  is  given  by 
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where  Q/ Q  =  (b  +  s)/(b  -  s)  ,  with  b  being  the  nonminimum  phase  factor. 
Since 
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^atECR2)]  •  — : •  — — 
b  +  a  b  -  s 


(23) 


E(he)R  =  4ab  ^  4(b/a) 

E(R2)  (b  +  a)2  [1  +  (b/a)]2 


(24) 


Equation  24  is  of  value  because  it  enables  us  to  state  limits  on  the 
performance  of  single-controller  systems  in  terms  related  to  terrain 
roughness.  For  example,  consider  two  types  of  terrain: 

(1)  Rolling  terrain,  with  a  spatial  correlation  distance 
ag  =  2,000  ft.  (as  is  the  average  distance  between 
two  samples  differing  in  magnitude  by  a  factor  of 
1/e  =  0.368.) 

(2)  Rough  terrain,  with  as  =  500  ft. 

For  a  vehicle  traversing  the  terrain  at  M  =  0.7  (781.9  fps)  the  corre¬ 
sponding  values  of  a  in  Eq.  24  are  a  *  781.9/2,000  **  0.39  sec  1  and 
a  *  781.9/500  ■  1.563  sec  for  rolling  and  rough  terrain  respectively. 
For  Configuration  3GA,  which  is  approximately  representative  of  current 
Tomahawk  class  missiles,  factoring  the  transfer  functions  given  in 
Table  2  yields  b  *  16.9799  rad/sec.  Hence  from  Eq.  24,  for  rolling 
terrain: 
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4  x  0.39  x  16.9799 


Minimum  Achievable 
Mean  Square  Height  Error 
Mean  Square  Terrain 
Height  Variation 


(0.39  +  16. 97 99) 2 


0.008  (25) 


For  rough  terrain: 


Minimum  Achievable 
Mean  Square  Height  Error 
Mean  Square  Terrain 
Height  Variation 


4— ■1---56-3-  16-979^  =  0.308  (26) 

(1. 563  +  16. 9799)2 


A  "worst  case"  type  of  terrain  would  be  one  yielding  a  =  b,  for  which 
Eq.  24  attains  its  maximum  value  of  1.0.  In  this  situation  the  optimal 
feedback  is  zero,  i.e.,  the  best  strategy  is  to  fly  the  missile  controls- 
fixed  in  a  straight-line  flight  path,  making  no  attempt  to  follow  the 
undulations  of  the  terrain.  For  a  Tomahawk  class  missile  flying  at  M  =  0.7 
the  spatial  correlation  distance  of  such  terrain  is  as  *  UD/b  *  46.0  ft. 
This  corresponds  to  exceedingly  jagged  terrain  which  would  rarely  be  en¬ 
countered.  The  "rough"  terrain  quoted  above  is  more  probable:  for  that 
terrain  rms  height  error  exceeds  50%  of  the  rms  terrain  altitude  variation. 

From  the  above  results  and  those  presented  in  Fig.  51,  it  is  concluded 
that,  for  aft  tail  single  controller  missiles  flying  over  land  or  sea  the 
maximum  achievable  closeness  with  which  the  terrain  can  be  followed  is 
severely  limited,  due  to  the  vehicle's  initial  wrong-way  height  response. 

Atmospheric  Turbulence  Response:  The  Dryden  spectrum  employed  to 
describe  low  altitude  turbulence  (Eq.  14)  is  widely  employed  in  aircraft 
guidance  and  control  studies.  The  mean  square  height  error  due  to  turbu¬ 
lence  was  much  smaller  than  that  due  to  terrain  for  all  the  configurations 
studied.  Therefore  the  effect  of  varying  the  Dryden  spectrum  appears  to 
be  secondary,  hence  no  parametric  variations  were  performed.  However 
three  points  of  detail  should  be  noted  for  future  studies. 

(1)  The  Dryden  spectrum  has  a  finite  d.c.  value,  which  may  sometimes 
cause  convergence  difficulties  for  optimization  procedures.  Since  the 
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ultra-low  frequency  behavior  is  of  little  interest  it  is  suggested  that 
the  spectrum  should  be  attenuated  at  low  frequencies  (e.g.,  below 
0.01  rad/sec)  to  remove  the  d.c.  component. 

(2)  The  use  of  gust-gradient  terms  (see  Eq.  10)  in  connection  with 
the  standard  Dryden  spectrum,  or  any  other  turbulence  spectrum  that  has 
a  finite  d.c.  value,  requires  consideration.  These  terms  provide  only  a 
first  approximation  to  the  complex  interactions  of  the  airframe  with 
time-varying  and  spatially-varying  atmospheric  velocity  components  (see 
Ref.  42  for  more  exact  but  more  complicated  representations).  As  a 
check,  we  repeated  most  of  our  optimization  calculations  with  the  gust 
gradient  terms  removed.  Only  minor  differences  in  the  results  were 
observed.  It  is  recommended  that  such  checks  should  always  be  performed. 

(3)  Our  calculations  show  that  the  height  error  induced  by  turbu¬ 
lence  is  much  smaller  than  that  due  to  terrain.  It  does  not  necessarily 
follow  that  the  turbulence- induced  normal  accelerations  are  much  less 
than  those  due  to  terrain-following.  If  subsequent,  more  detailed  studies 
are  performed  the  "g"  response  of  the  optimally  controlled  system  should 
be  assessed.  The  performance  index  can  be  extended,  if  desired,  to 
include  terms  relating  to  normal  acceleration. 

Terrain-Following  Performance  of  Single-Controller  Joined  Wing  Con¬ 
figuration:  The  joined  wing  configurations  3F,  6F,  9F  obtain  pitch  control 
via  wing-warping,  as  described  in  Section  3.  With  a  parabolic  variation 
of  twist,  corresponding  to  uniform  torsional  modulus  GJ  along  both  wings, 
the  center  of  pressure  of  the  lift  induced  by  twist  is  ahead  of  the  c.g., 
and  hence  the  height/twist  transfer  function  is  minimum  phase.  However, 
as  is  apparent  from  Fig.  32,  the  pitching  moment  arm  is  small  for  Con¬ 
figuration  3F.  For  the  hybrid  configuration,  6H  (Fig.  34) ,  a  long  moment 
arm  is  obtained  by  employing  an  all-moving  canard  as  the  primary  pitch 
control. 

Figure  51  shows  the  terrain-following  performance  of  these  configu¬ 
rations  when  optimally  controlled.  The  gust  and  wave  spectra  are  identi¬ 
cal  to  those  employed  previously.  Comparing  Fig.  52  with  Fig.  50  it  is 
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JOINED  WING  CONFIGURATIONS 


seen  that,  as  predicted,  the  joined  wing  and  hybrid  configurations  yield 
a  mean  square  height  error  which  tends  to  zero  as  k  -*•  0.  Unfortunately, 
the  associated  magnitudes  of  control  deflection  E(6^)  are  impractically 
large.  Two  remedies  are  possible:  (1)  increase  control  effectiveness; 

(2)  employ  multiple  controls  (e.g.,  flaps  plus  elevator). 

Increased  control  effectiveness  implies  increasing  both  and 

to  n  times  their  original  levels.*  This  yields  a  mean  square  height  error 

2 

equal  to  that  provided  by  (1/n)  times  the  original  control  deflection  of 
the  original  optimized  system,  for  a  given  value  of  k.  Thus  for  the  aft 
tail  configuration  the  limit  on  terrain  following  accuracy  remains  un¬ 
changed;  however,  it  is  approached  more  closely.  The  chances  of  achiev¬ 
ing  substantially  increased  control  effectiveness  are  slim,  although  some 
gains  could  be  made  by  the  following  expedients: 

(1)  (For  the  aft  tail  and  hybrid  configurations)  gearing 
an  anti-servo  tab  to  the  all-moving  tail  (or  canard) 
is  done  on  many  light  planes  (e.g.,  the  Piper  Cherokee 
series) . 

(2)  (For  the  joined  wing  configurations)  augmenting  or  re¬ 
placing  the  wing  warping  control  with  deflections  of 

a  flap  mounted  on  the  front  wing. 

On  balance  it  appears  to  be  more  desirable  to  employ  a  multicontroller 
system,  as  described  below.  This  provides  large  improvements  in  terrain¬ 
following  ability  for  the  hybrid  and  joined  wing  configurations,  but  only 
smaller  improvements  for  the  aft  tail  configurations. 


*The  effect  of  reducing  static  margin  to  zero  was  also  studied.  This 
yielded  negligible  benefits  because,  as  shown  in  Fig.  53,  it  does  not 
change  the  h/6  transfer  function  much  in  the  frequency  region  around 
6.3  rad/ sec  where  the  wave  spectrum  has  its  peak.  Reducing  static  mar¬ 
gin  would  be  effective  if  the  command  (terrain)  spectrum  is  of  the  type 
that  has  large  low-frequency  components,  e.g.,  land  terrain. 
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FIGURE  53.  HEIGHT/CANARD  DEFLECTION  TRANSFER  FUNCTION  FOR  HYBRID  CONFIGURATIONS 


4 . 4  MULTI  CONTROLLER  SYSTEMS:  HYBRID  CONFIGURATION 


We  use  the  term  "multicontroller’’  to  denote  systems  employing  two 
separate  controls,  e.g.,  flaps  plus  elevator.  The  benefits  of  multicon¬ 
troller  systems  are  most  easily  understood  by  first  ascertaining  the 
reason  for  the  large  control  deflections  demanded  by  the  single-controller 
systems.  Consider  Fig.  53  which  shows  the  h/6  transfer  function  for  Con¬ 
figuration  6H,  where  6  =  deflection  of  the  all-moving  canard,  measured  in 
degrees.  The  equation  for  h/5,  from  Table  2,  is 

JL  _  1.318s  4-1. 287 687s +  500.441103  =  1 ,318(s  +  0.4885  t  19.4797j) 

6  s2(s2 + 2. 083s+  6.019988)  s2 (s2 + 2 . 083s +  6.019988) 

The  Bode  diagram  of  Eq.  25,  given  on  Fig.  53,  exhibits  a  pronounced  dip 
around  19.5  rad/sec,  due  to  the  light  damping  of  the  numerator  factors. 
This  dip  unfortunately  attenuates  the  h/6  response  in  the  region  where  it 
is  most  needed,  i.e.,  in  the  region  around  the  peak  of  the  wave  spectrum. 
If  the  dip  is  eliminated  (or  smoothed)  the  response  improves  dramatically. 
This  can  be  demonstrated  by  putting 

h  100 

■  '  ■■  s  - 


Below  frequencies  of  2  rad/sec  this  has  an  amplitude  characteristic 
similar  to  the  basic  h/6  transfer  function  (Eq.  27).  However,  its  ter¬ 
rain-following  performance  is  much  superior,  as  shown  in  Fig.  54,  which 
graphs  the  results  of  optimal  control  calculations  similar  to  those 

illustrated  in  Figs.  50  and  52.  Note  that  close  terrain  following  can 

2  2 

be  achieved  with  moderate  control  deflections,  e.g.,  E(hg)  =  45.1  ft 
with  E(62)  *  6.8  deg2.  For  zero  static  margin  the  results  are  even 
better:  removing  the  dip  from  the  zero  static  margin  line  on  Fig.  53 
yields : 
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FIGURE  54.  TERRAIN-FOLLOWING  PERFORMANCE  OF  MULTI CONTROLLER  SYSTEMS 
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Figure  54  illustrates  that  this  achieves  E(he)  *  21.9  ft^  with  only 
1.00  deg2  of  mean  square  canard  deflection. 

The  desirable  forms  of  transfer  functions,  presented  above,  can  be 
realized  by  employing  two  controls.  This  will  be  demonstrated  for  the 
hybrid  configuration  6H  controlled  via  canard  incidence  and  front  wing 
flaps.  It  will  also  be  shown  that,  for  aft  tail  configurations,  desir¬ 
able  forms  of  transfer  functions  cannot  be  physically  realized  through 
the  combined  use  of  aft  tail  and  wing  flap  deflections. 

Consider  Fig.  55.  This  shows  a  multicontroller  system  suitable  for 
Configuration  6H,  in  which  two  controls  are  actuated  by  the  optimal  con¬ 
trol  system.  These  controls  are  canard  incidence  (<$j  degrees)  and  deflec¬ 
tion  of  front  wing  flaps  (62  degrees).  The  flap  control  derivatives  are 
M52 »  zfi2’  and  the  crossfeed  in  Fig.  55  has  transfer  function  C(s) .  It 
is  desired  that  the  overall  transfer  function  of  the  blocks  enclosed  by 
dashed  lines  should  be  100/s2.  Then,  from  Eqs.  11  and  15: 


100 


cz<52  Is2  -  %  +  Md)  s  -  [Ma  -  (M62/Z62) ZJ  }  +  Ah(5ls2  +  B^s  +  Ch(5l 


s2 ( s2  +  Bs  +  C) 


(31) 


Substituting  numerical  values  for  Configuration  6H  from  Table  2,  and 
manipulating  Eq.  31  yields  the  required  crossfeed  transfer  function  as: 


_ 98.682s2  +  207 . 0123s  +  101 .5649 _ 

-Z62|s2  +  0.977s+  [5.248  -  864.77(M<52/Z<52)  ] } 


(32) 


C  must  not  be  unstable,  if  it  is  to  be  constructed  physically.  This 
condition  implies  that 


5.248  -  864.77(M,52/Z62)  >  0 


MULTI CONTROLLER 


\ 


A  controller  located  forward  of  the  c.g.,  such  as  inboard  flaps  on  the 
front  wing,  has  (M5 2/ ^62 ^  <  ®  an^  satisfies  this  condition  automatically. 
For  trailing  edge  flaps  occupying  the  entire  span  of  the  front  wing 
Mg ^  ~  0,  and  2,5  2  =  -4.31,  whence: 


98 . 682 ( s  4-  0.7821)  (s  +  1.3157) 
4 . 31 (s  +  0.4885  ±  2.2382j) 


(33) 


Figure  56  graphs  the  Bode  diagram  of  this  transfer  function.  In  the 
frequency  range  of  most  interest  (1  to  10  rad/sec)  the  amplitude  ratio 
is  typically  30  dB,  i.e.,  a  linear  ratio  of  31.6.  From  Fig.  54  this 
implies  that  to  obtain  a  mean  square  height  error  of,  for  example, 

74.87  ft^  the  canard  rms  deflection  would  be  /0.0369  *  0.192  deg  and  the 
rms  flap  deflection  would  be  of  the  order  of  31.6  *  0.192  =  6.07  deg. 

This  performance  is  considerably  better  than  that  of  the  single  control¬ 
ler  system.  Still  further  gains  could  be  achieved  by  tailoring  the  form 
of  C  and/or  reducing  static  margin.  These  gains  would  involve  per¬ 
mitting  higher  rms  canard  deflections,  e.g.,  5  degrees,  while  maintain¬ 
ing  rms  flap  deflection  at  approximately  6  degrees.  This  could  be  done 
either  by  applying  optimal  control  theory  for  multicontroller  systems, 
or  by  iterative  classical  methods  of  control  systems  design. 

In  summary,  the  use  of  combined  flap  and  canard  controls  enables 
the  hybrid  configuration  to  achieve  substantially  better  terrain-following 
performance  than  can  be  attained  using  canard  control  alone.  Similar 
benefits  cannot  be  obtained  for  the  aft  tail  configuration  through  the 
combined  use  of  flaps  and  elevator,  as  will  be  shown  below. 

Multicontroller  Aft  Tail  Configuration:  The  inherent  limitation 
of  aft  tail  configurations  for  close  terrain  following  with  combined 
control  flap  and  elevator  (or  all-moving  tail)  becomes  apparent  upon 
attempting  to  repeat  the  procedure  demonstrated  above  for  the  hybrid 
configuration,  whereby  an  effective  transfer  function  of  the  desired 
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FIGURE  56.  BODE  DIAGRAM  OF  CROSSFEED  TRANSFER  FUNCTION 


5 


form  (Eq.  28)  was  realized.  Configuration  6GA  is  selected  to  provide  a 
numerical  example. 


Substituting  the  appropriate  numbers  for  Configuration  6GA  (from 
Table  2)  in  Eq.  31  yields  an  equation  for  the  crossfeed  transfer  func¬ 
tion,  C,  which  relates  flap  deflection  to  aft  tail  deflection.  Equa¬ 
tion  31  thus  becomes,  for  Configuration  6GA, 


-CZ5;3Js2+  0.9468s +  3  6.14  -  1 ,790.  9(Mg3/Z(s3)[ 

100  =  _ -  3.81859s2  -  3.615s  +  1,883.92 

s2  s2(s  +  3.2372s+  37.7993948) 


(34) 


where  63  denotes  flap  deflection. 

Manipulating  Eq.  34  similarly  to  Eq.  31  yields  the  following  neces¬ 
sary  condition  for  C  to  be  stable 


i.e. , 


M63 

36.14  -  1,790.0  ~  >  0 

4<$3 


% 

z63 


< 


0.02017 


(35) 


(36) 


Let  the  center  of  pressure  of  the  lift  due  to  flap  deflection  be  located 
a  distance  &F  aft  of  the  c.g.  By  definition: 


M6.  W  2F 

Z53  g  Iy 


2,396.8  % 

32.2  '  1,815.6 


0.4099  eF 


(37) 


Hence  the  maximum  value  of  £p  is  given  by 

Up.)  -  0.02017/0.04099  -  0.49  feet 

*  max 


(38) 


This  distance  is  0.317  times  the  mean  wing  chord  of  Configuration  6GA. 
Thus  if  the  vehicle  c.g.  is  at  0.25  c,  for  C  to  be  stable  the  incre¬ 
mental  lift  due  to  flap  deflection  must  have  its  center  of  pressure 
forward  of  0.567  c.  As  shown  in  Ref.  43,  p.  244,  most  trailing  edge 
flap  configurations  yield  incremental  center  of  pressure  locations 
which  are  further  aft  than  this  value.  Offsetting  this,  the  tail  load 
induced  by  flap  downwash  moves  the  c.p.  forward,  but  only  after  a  lag 
of  approximately  0.01  seconds  (the  time  required  for  the  downwash  to 
reach  the  tail).  Hence,  even  if  the  aft-mounted  control  attains 

<  0.49  feet,  it  still  has  a  slower  response  than  a  forward-mounted 
control.  Note  that  should  be  substantially  less  than  0.49  feet, 
because  if  the  inequality  condition  (Eq.  36)  is  only  just  satisfied  C 
will  be  large  in  the  region  to  s  [36.14  -  1 , 790.  OCM^/Z^)  ]  and  the 
system  will  be  impracticable  because  it  will  demand  excessively  large 
flap  deflections. 

It  is  interesting  to  examine  the  effect  of  reducing  static  margin 
on  the  feasibility  of  aft  tail  multicontroller  systems.  If  we  let  the 
static  margin  become  zero  then  in  Eq.  34  the  Mq  term  (=  -36.14)  vanishes 
As  a  result  the  stability  condition  becomes  M53/Z53  <  0,  which  condition 
cannot  be  satisfied  by  any  trailing  edge  flap  system.  Thus  reducing 
static  margin  only  worsens  the  situation  for  the  aft  tail  airframe. 

Multicontroller  Non-Hybrid  Joined  Wing  Configurations:  The  hybrid 
is  more  attractive  than  the  pure  joined  wing  configuration  for  multicon¬ 
troller  terrain  following,  because  of  the  long  moment  arm  of  the  canard. 
If  staging  is  employed,  as  shown  in  Figs.  34,  42,  and  43,  the  canard  is 
no  longer  available  to  the  post-staged  configuration.  To  obtain  close 
terrain  following  after  staging  has  occurred  it  may  be  desirable  to 
mount  the  canard  aft  of  the  nose  tank,  so  that  the  vehicle  remains  a 
hybrid  configuration  throughout  its  flight. 


If  Che  canard  cannot  be  used,  consideration  should  be  given  to: 

(1)  Combined  deflections  of  rear  wing  flaps  and  front 
wing  flaps,  with  the  larger  deflections  on  the 
front  wing  flaps,  to  keep  the  c.p.  of  the  added 
airload  forward  of  the  c.g. 

(2)  Movable  leading  edge  devices  (e.g.,  Kruger  flaps, 
controllable  L.E.  droop)  mounted  on  the  inboard 
part  of  the  front  wing. 


Canard  Configurations.  It  may  be  asked  whether  a  canard  configura¬ 
tion  would  follow  terrain  as  closely  as  the  hybrid  configuration  6H. 
Penalties  of  structure  weight,  induced  drag,  and  increased  wing  volume 
would  ensue  if  the  joined  wing  of  6H  were  replaced  by  a  monoplane  wing. 
Even  if  these  penalties  are  disregarded,  the  control  moment  arm  available 
from  such  a  wing  would  be  less  than  that  provided  by  a  comparable  joined 
wing.  Furthermore  flaps  mounted  on  such  a  wing  would  be  aft  of  the  c.g., 
thus  producing  an  initially  wrong-way  response  which,  as  has  been  shown, 
is  detrimental  to  terrain-following.  It  therefore  appears  that  no  advan¬ 
tage  would  be  gained  by  adopting  a  canard  configuration  employing  canti¬ 
lever  lifting  surfaces. 

4.5  CONCLUDING  REMARKS 


This  section  has  shown  how  configuration  size  and  geometry  affect 
terrain-following  accuracy.  Seven  airframe  configurations  were  studied, 
but  only  a  single  type  of  terrain  was  considered:  a  short-crested  sea. 

It  would  be  of  interest  to  apply  the  techniques  described  here  to  terrain¬ 
following  over  land.  One  would  expect  greater  improvements  to  be  possible 
over  land  because  of  the  larger  amplitude  and  smaller  high-frequency  con¬ 
tent  of  land  terrain  spectra. 

The  system  models  presented  here  are  simplified  in  that  they  omit 
such  things  as  actuator  lags,  digital  processing  lags,  radar  altimeter 
noise,  and  airframe  structural  dynamics.  Because  of  these  simplifica¬ 
tions  the  accuracy  of  terrain-following  calculated  here  tends  to  be 
optimistic.  Further,  more  detailed,  studies  should  be  performed  to 
study  these  effects. 


L 


118 


5.0  PERFORMANCE  COMPARISON 


5.1  INTRODUCTION 

This  section  presents  comparisons  of  calculated  range  and  maneuvera¬ 
bility  for  the  joined  wing,  conventional  and  hybrid  configurations  3GA, 
6GA,  9GA,  3F,  6F,  9F,  and  6H  shown  in  Figs.  31  through  34.  The  compari¬ 
sons  are  based  on: 

(1)  Range  at  sea  level  at  M  *  0.7. 

(2)  Range  at  40,000  ft  altitude  at  M  *  0.7. 

(3)  Sustained  maneuvering  "g"  capability  at  sea 
level,  M  *  0.7. 

All  configurations  are  assumed  to  have  an  initial  cruise  weight  of 
2864.0  lb.  This  is  the  configuration's  weight  when  the  booster  separates, 
hence  all  configurations  employ  the  same  size  of  booster.  The  assumed 
powerplant  characteristics  are  generally  representative  of  projected 
small  turbofans  of  the  class  used  in  current  Tomahawk  cruise  missiles, 
but  do  not  correspond  to  any  specific  existing  engine. 

The  geometry  of  the  airframe  configurations  has  been  described  in 
Section  3.  A  point  that  requires  additional  comment  is  the  increase  in 
the  ratio  of  horizontal  projected  tail  area  S?  to  wing  area  for  Configu¬ 
rations  6GA  and  9GA,  compared  to  Configuration  3GA.  This  is  required  to 
maintain  adequate  tail  volume  ratio  v  -  J^tS^/Sc.  The  mean  chord,  c,  is 
proportional  to  because  a  constant  aspect  ratio  is  retained.  The 

tail  moment  arm  is  approximately  constant,  since  fuselage  length  is 
fixed.  For  constant  v,  S?  would  be  proportional  to  S  x  i.e.,  S^"^. 

However  since  part  of  the  function  of  the  tail  is  to  overcome  fuselage 
destabilizing  moments,  and  since  the  fuselage  size  is  fixed,  a  less 
extreme  variation  of  tail  size  was  adopted:  *  S**^.  This  gives  the 

tail  areas  tabulated  overleaf  and  illustrated  on  Fig.  31.  For  Configu¬ 
rations  6GA  and  9GA  the  curved  folding  tail  of  3GA  (Fig.  33)  is  not 
used:  instead  a  Y-tail  arrangement  having  straight  surfaces  spaced  at 

120  degrees  is  employed,  with  multiple  hinges  for  folding. 
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3GA 

6GA 

9GA 

Gross  horiz.  tail  pro j .  area,  ft^ 

4.83 

10.868 

17.51 

Gross  wing  area,  ft^ 

11.5 

12.0 

34.5 

Tail  area/wing  area 

0.42 

0.4725 

0.5075 

£-p/ mean  geometric  chord 

7.57 

4.569 

3.92 

Tail -volume  ratio 

3.18 

2.16 

1.99 

The  vertical  tail  area  was  varied  similarly  to  the  horizontal  tail  area. 
(An  increase  in  wing  span,  even  with  no  fuselage  change,  requires  an 
increase  in  vertical  tail  area  for  proper  turn  coordination,  as  shown 
in  Ref.  46.) 

For  drag  estimation  purposes  it  was  assumed  that  all  configurations 
had  identical  flush  inlets.  As  noted  previously  the  joined  wing  and 
hybrid  concepts  can  be  employed  with  either  flush  or  protruding  inlets, 
as  desired. 

5 . 2  DRAG  ESTIMATION 


The  procedure  employed  for  estimating  drag  is  generally  similar  to 
the  methods  of  Linden  and  O'Brimski  (Ref.  44)  and  Torenbeek  (Ref.  24) 
in  that  (1)  wetted  area  and  form  factors  are  employed  as  a  basis  for 
calculating  minimum  drag,  and  (2)  an  offset  drag  polar  is  used  with 
minimum  drag  at  Cl  “  Ci^  i  0. 


Cu  *  Cl 


m^n 


(cL  -  ci^)2 


(39) 


The  accuracy  of  the  method  was  checked  by  computing  the  drag  of  the 
Vought  YGBM-110  Tomahawk  missile.  Full-scale  wind  tunnel  measurements 
taken  on  an  actual  YGBM-110  at  H  *  0.7  S.L.  conditions  are  given  in 
Ref.  45.  Figure  57  compares  the  measured  drag  polar  with  the  results 
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of  our  calculations.*  This  figure  shows  that  although  the  predicted 
variation  of  drag  with  lift  is  reasonably  accurate,  the  calculated  CDmin 
is  approximately  1.1  times  the  measured  value.  This  can  be  ascribed  to 
our  assumption  of  100%  turbulent  flow,  i.e.,  transition  at  all  leading 
edges.  It  would  have  been  simple  to  modify  the  method  by  assuming  tran¬ 
sition  at  10%  chord  on  all  surfaces,  which  assumption  would  yield  results 
closely  matching  the  wind  tunnel  data.  However  it  was  decided  to  retain 
the  assumption  of  fully  turbulent  flow  (1)  for  conservatism  and  (2)  to 
allow  for  the  "drag  rectification"  effect  caused  by  positive  and  negative 
control  deflections  which  cause  a  small  drag  increase  regardless  of  sign. 
A  further  consideration  is  that  the  surface  finish  of  an  operational 
vehicle  may  be  inferior  to  that  of  a  factory-fresh  vehicle  selected  for 
wind  tunnel  testing. 

Table  3  shows  the  calculated  minimum  drag  at  Sea  Level  Standard 
Conditions  for  the  joined  wing  and  conventional  configurations:  the 
hybrid  configuration  results  are  given  on  Table  4.  For  the  40,000  ft 
condition  the  minimum  drag  values  were  calculated  by  varying  the  weighted 
Ryenolds  numbers  shown  on  Tables  3  and  4  proportionately  to  the  Re/ft  of 
the  high-altitude  flight  condition. 

Drag  Due  to  Lift:  For  all  configurations  drag  due  to  lift  is 
expressed  as: 

(cL -cl*)  2 

CD  -  Cd  min  *  - -  for  CL  <  CLg  (40) 


(Cl  -  CLv  +  cl-cls)2 

Cd  -  CD  min  *  - ^77 - ~  for  CL  CLg  (41) 


*The  YGBM-110  has  a  fuselage  diameter  of  20.875  in.  Configuration  3GA 
is  generally  similar,  but  has  a  smaller  fuselage  diameter  (20.375  in.). 
This  diameter  was  selected  to  equal  that  of  the  General  Dynamics  YGBM-109, 
which  employs  a  protective  steel  casing  for  underwater  launch,  and  there¬ 
fore  is  appreciably  smaller  in  diameter  than  the  21  inch  torpedo  tube 
from  which  it  is  launched. 
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TABLE  3.  MINIMUM  DRAG  AT  M  =  0.7,  SEA  LEVEL  STANDARD  CONDITIONS 


Conf ig. 

Component 

^wet 

(ft2) 

Re  x  10~6 

cf 

Turbulent 
Flat  Plate 

Form 

Factor 

Cfe 

3F 

Front  wing 

10.059 

2.85 

0.0037 

1.267 

0.0047 

Rear  wing 

9.680 

2.56 

0.0038 

1.267 

0.0048 

Fin 

8.800 

6.2 

0.0032 

1.267 

0.0041 

Fuselage 

87.769 

89.3 

0.00216 

1.324 

0.0028 

TOTAL 

116.308 

68.32* 

0.003274* 

6F 

Front  wing 

18.870 

4.05 

0.00349 

1.267 

0.0044 

Rear  wing 

13.360 

3.63 

0.00355 

1.267 

0.0045 

Fin 

8.914 

6.22 

0.00324 

1.267 

0.0041 

Fuselage 

87.769 

89.3 

0.00216 

1.324 

0.0028 

TOTAL 

128.91 

62.2* 

0.00332* 

9F 

Front  wing 

29.148 

4.985 

0.003366 

1.267 

0.00427 

Rear  wing 

29.040 

4.44 

0.003432 

1.267 

0.00435 

Fin 

13.37 

7.59 

0.00314 

1.267 

0.003 98 

Fuselage 

87.769 

89.3 

0.00216 

1.324 

0.0028 

TOTAL 

159.327 

51.6* 

0.00345* 

3GA 

Wing 

20.6 

5.32 

0.00333 

1.267 

0.00422 

Tail 

10.272 

5.30 

0.00333 

1.267 

0.00422 

Fuselage 

87.769 

89.3 

0.00216 

1.324 

0.0028 

TOTAL 

118.641 

67.44* 

0.00317* 

6GA 

Wing 

41.2 

7.45 

0.00315 

1.267 

0,00399 

Tail 

24.43 

8.86 

0.00306 

1.267 

0.00388 

Fuselage 

87.769 

89.3 

0.00216 

1.324 

0.0028 

TOTAL 

153.399 

54.5* 

0.00329* 

9GA 

Wing 

61.8 

9.197 

0.00304 

1.267 

0.00386 

Tail 

40.56 

10.77 

0.00297 

1.267 

0.00376 

Fuselage 

87.769 

89.3 

0.00216 

1.324 

0.0028 

TOTAL 

190.129 

46.51* 

0.00335* 

*"Total"  indicates  weighted  Re  =  E  Re  Swet/E  Swet  and 
similarly  weighted  Cfe- 


Dmin 

(lb) 


34.3 
33.6 

26.3 
182.4 

276.6 


60.6 

43.7 

26.6 

182.4 

313.3 


90.4 

91.8 

38.7 

182.4 

403.3 


63.2 

31.5 

182.4 

277.1 


119.5 

68.9 

182.4 

370.8 


173.3 
110.9 

182.4 

466.6 


TABLE  4 
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HYBRID  CONFIGURATION  MINIMUM  DRAG  AT  M  =  0.7,  SEA  LEVEL  STANDARD  CONDITIONS 


Conf ig. 

Component 

swet 

(ft2) 

Re  x 10“6 

Cf 

Turbulent 
Flat  Plate 

Form 

Factor 

Cfe 

Dmin 

(lb) 

Initial 

Fuselage 

87.769 

89.3 

0.00216 

1.35 

0.0029 

186.0 

Cruise 

Canard  ^ 

4.112 

4.17 

0.00347 

1.267 

0.0044 

13.1 

Front  wing 

19.190 

3.58 

0.00356 

1.267 

0.0045 

62.9 

Rear  wing 

13.226 

2.48 

0.00379 

1.267 

0.0048 

46.2 

Fin 

3.006 

6.29 

0.00324 

1.267 

0.0041 

9.0 

TOTAL 

127.303 

62.65* 

0.00342* 

317.2 

Before 

Fuselage 

87.769 

89.3 

0.00216 

1.35 

0.0029 

186.0 

Staging 

Canard  § 

2.93 

5.65 

0.0033 

1.267 

0.00418 

8.9 

Front  wing 

19.190 

3.58 

0.00356 

1.267 

0.0045 

62.9 

Rear  wing 

13.226 

2.48 

0.00379 

1.267 

0.0048 

46.2 

Fin 

3.006 

6.29 

0.00324 

1.267 

0.0041 

9.0 

TOTAL 

126.12 

63.14* 

0.00340* 

313.0 

After 

Staging 

Fuselage 
Front  wing 
Rear  wing 
Fin 

67.851 

19.190 

13.226 

3.006 

70.6 

3.58 

2.48 

6.29 

0.00224 

0.00356 

0.00379 

0.00324 

1.324 

1.267 

1.267 

1.267 

0.00296 

0.0045 

0.0048 

0.0041 

150.0 

62.5 

42.7 

8.96 

TOTAL 

103.273 

47.55* 

0.00352* 

264.2 

*"Total"  indicates  weighted  Re  «  E  Re  Swet/E  Swet  and 
similarly  weighted  Cfe. 

A  «  35°,  at  1/4-chord  line. 

®  A  “  62°,  at  1/4-chord  line. 
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where  Cls  is  the  lift  coefficient  at  which  massive  separation  occurs. 

Cls  is  approximately  equal  to  the  buffet  lift  coefficient:  it  was  found 
that  almost  all  the  range  and  maneuver  results  of  interest  here  did  not 
involve  flight  at  Cl  >  Cls,  so  Eq.  AO  was  adequate.  Table  5  lists  the 
key  parameters  describing  drag  due  to  lift.  For  the  joined  wing  configu¬ 
rations  ex  and  Clx  were  estimated  using  wind  tunnel  test  data  as  described 
in  Section  1.  Note  that,  before  staging,  the  hybrid  configuration  has  an 
ex  only  slightly  better  than  that  of  the  conventional  configurations: 
this  is  because  the  lift  carried  by  the  canard  is  distributed  over  a  much 
shorter  span  than  of  the  joined  wing.  The  canard  lift  decreases  as  the 
nose  tank  is  emptied,  but  its  span  also  decreases  because  its  sweep  must 
be  increased  throughout  cruise  to  maintain  a  positive  static  margin.  After 
staging,  the  hybrid  regains  the  full  induced  drag  benefits  of  the  joined 
wing. 


TABLE  5.  DRAG-DUE -TO-LIFT  PARAMETERS 


Configu¬ 

ration: 

3GA 

3F 

6GA 

6F 

* 

6Hi 

* 

6H2 

9GA 

9F 

CLX 

0.066 

0.126 

0.066 

0.126 

0.126 

0.126 

0.066 

0.126 

ex 

0.528 

0.755 

0.528 

0.755 

0.528 

0.766 

0.528 

0.755 

Span,  b, 
ft 

10.5 

9. 1395 

1A.8A 

12.925 

12.925 

12.925 

18.18 

15.83 

Reference 
area ,  S , 
ft2 

11.5 

10.0 

23.0 

20.0 

15.9992 

15.9992 

34.5 

30.0 

A 

9.6 

8.35 

9.6 

8.35 

10.44 

10.44 

9.6 

8.35 

CLS 

0.9 

1.1 

0.9 

1.1 

1.1 

1.1 

0.9 

1.1 

*Subscripts  1  and  2  denote  pre-staging  and  post-staging  configurations, 
respectively. 
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Figure  58  presents  calculated  drag  polars  and  (L/D)max  values  for 
all  the  above  configurations. 

5.3  ENGINE  PERFORMANCE 


As  noted  previously,  the  assumed  specific  fuel  consumption  charac¬ 
teristics  are  representative  of  turbofan  engines  projected  for  advanced 
cruise  missiles,  but  do  not  correspond  to  the  S.F.C.  characteristics  of 
any  particular  engine.  For  flight  at  M  =  0.7  S.L.  the  following  equa¬ 
tion  for  S.F.C.  was  assumed: 


(lb  Fuel /lb  Thrust /hr)  ‘  1-57597295  -  0.30258  (-jL)  +  0.0369  (j^) 

(42) 

when  T  *  installed  thrust  in  lb.  This  gives  a  S.F.C.  of  1.31  at  T  =  100  lb, 
reaching  a  minimum  of  0.955  at  T  =  410  lb.  For  flight' at  M  *  0.7,  40,000  ft 
the  S.F.C.  was  taken  as  0.842  times  the  value  calculated  from  the  above 
equation.  The  maximum  installed  thrust  at  M  =  0.7,  S.L.  was  assumed  to  be 
500  lb. 

5.4  ESTIMATED  WEIGHTS 


Table  6  gives  weight  breakdowns  for  each  configuration.  Standard 
weight  estimation  procedures  (Ref.  24)  plus  published  data  were  employed 
to  derive  component  weights  for  the  conventional  configurations.  The 
estimates  given  are  necessarily  approximate  —  this  is  particularly  true 
for  Configuration  9GA  where  the  wing  folding  system  is  bound  to  be  com¬ 
plex  in  view  of  the  large  dimensions  of  the  wing  relative  to  the  fuse¬ 
lage.  The  weight  penalty  for  such  a  folding  system  can  only  be  determined 
accurately  after  it  has  been  designed  in  detail.  For  current  cruise  mis¬ 
siles  wing  weight  is  determined  by  stiffness  rather  than  static  strength 
requirements.  In  view  of  the  high  stiffness  of  the  joined  wing  (see 
Fig.  8)  this  should  permit  considerable  weight  saving  compared  to  the 
conventional  wing -plus-tail .  It  has  been  assumed  here  that  the  joined 
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Conflg.  (L/DJmax  Cl  at  (L/D) 


AND  6H  (BEFORE  AND  AFTER  STAGING) 


TABLE  6.  WEIGHT  SUMMARY 


"Remainder"  includes:  propulsion  (226),  power  generation  (AA)  ,  orientation  (A3), 
guidance  (113),  environmental  system  (8),  armament  (3),  and  miscellaneous  (30.5), 
totalling  A67.5  lb. 


wing  weight  is  equal  to  that  of  the  conventional  wing.  From  the  analysis 
of  Ref.  3  this  assumption  is  quite  conservative. 

The  fuel  weights  listed  in  Table  5  yield  initial  cruise  weights  of 
2864.0  lb  for  each  configuration.  On  this  basis  the  fuel  weights  for 
the  large-wing  configurations  are  obviously  less  than  those  for  con¬ 
figurations  employing  smaller  and  lighter  wings.  Because  of  this  fuel 
weight  reduction,  adequate  fuel  tank  volume  is  available  even  for  the 
large-winged  configurations  where  much  fuselage  volume  is  occupied  by  the 
folded  wings.  If  the  initial  cruise  weight  were  allowed  to  increase, 
volumetric  fuel  capacity  would  become  a  limiting  factor,  particularly 
for  the  conventional  configurations.  This  is  discussed  in  Section  5.5. 

5.5  RANGE  AND  MANEUVERABILITY  COMPARISONS 

Since  payload  and  Mach  number  have  been  fixed  for  all  configurations, 
the  major  performance  comparisons  relate  to  range  and  sustained  g's. 

These  have  been  calculated  for  each  configuration,  and  the  results  are 
compactly  presented  on  Figs.  59  and  60.  In  these  figures  "sustained  g's" 
denotes  the  normal  acceleration  at  which  the  drag  equals  the  maximum 
installed  thrust  (500  lb) .  The  weight  assumed  for  calculating  sustained 
g's  equals  the  weight  of  the  missile,  plus  payload,  oil,  and  10%  of  the 
initial  cruise  fuel.  In  both  Figs.  59  and  60  the  sustained  g's  shown 
are  appropriate  to  M  =  0.7  at  sea  level,  however  in  Fig.  60  the  range 
shown  applies  for  cruise  at  M  =  0.7,  40,000  ft.  The  ranges  given  do  not 
include  allowances  for  headwinds,  loiter,  climb  and  descent,  and  unusable 
fuel  (which  is  typically  1%  of  total  fuel  tankage) . 

Figure  59  shows  that  for  the  conventional  configuration  increasing 
wing  area  does  not  significantly  increase  sustained  g  capability  (although 
it  raises  non-sustained  g  capability).  This  is  because  the  parasite  drag 
of  the  increased  wing  area  approximately  cancels  the  reduction  of  induced 
drag  obtained  through  the  reduced  span-loading.  This  result  is  of  course 
dependent  on  altitude  and  also  on  the  maximum  engine  thrust  assumed,  but 
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SUSTAINED  ‘g's  AT  SEA  LEVEL 


FIGURE  59.  MANEUVERABILITY  VS.  RANGE  AT  SEA  LEVEL 
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FIGURE  60.  MANEUVERABILITY  VS.  RANGE  AT  40,000  FT  ALTITUDE 
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for  sea  level  maneuvering,  changes  of  ±10%  in  maximum  thrust  will  not 
greatly  affect  this  conclusion.  To  improve  maneuverability  it  is  neces¬ 
sary  to  employ  either  a  higher  aspect  ratio  monoplane  configuration,  or 
a  joined  wing  or  hybrid  configuration,  in  order  to  get  lower  induced 
drag.  For  1800  n.m.  range  the  gain  in  maneuverability  offered  by  the 
joined  wing  and  hybrid  configurations  is  considerable,  as  shown  below: 


Configuration 

Conventional 

Joined  Wing 

Hybrid 

Sustained  g 

A  7 

capability 

J  •  J 

J  •  J 

Ratio 

1.0 

1.34 

1.57 

Figure  60  illustrates  the  high  altitude  range  performance  of  the 
configurations  studied.  It  is  necessary  to  employ  wing  areas  of  at 
least  20  ft^  in  order  to  stay  below  the  buffet  limit  (Cls)  during  the 
initial  cruise.  The  low  induced  drag  of  the  joined  wing  configuration 
6F  is  the  main  factor  providing  1.38  times  the  range  of  its  conventional 
counterpart  (6GA) .  The  hybrid  configuration  6H  has  a  maneuverability 
advantage,  but  no  range  advantage  over  6F.  This  is  because  only  a  single 
hybrid  configuration  was  studied.  The  maneuverability  advantage  can  be 
traded  for  a  range  advantage  by  varying  the  hybrid  wing  area.  Parametric 
variations  of  canard/ front /rear  wing  geometry  would  also  yield  improved 
range  for  the  hybrid  configuration,  since  Configuration  6H  was  sized  for 
stability  and  control,  and  convenient  folding  rather  than  for  maximum 
range  or  maneuverability. 

The  range  advantage  of  the  joined  wing  configuration  can  be  summar¬ 
ized  by  comparing  the  6F  and  6GA  ranges  in  nautical  miles  graphed  on 
Figs.  59  and  60,  thus: 
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Range  of  Joined  Wing  Config. 
Range  of  Conventional  Config. 


1636 

1206 


3020 

2166 


1.36  at  sea  level 


1.39  at  40,000  ft 


The  advantage  of  the  joined  wing  decreases  as  wing  area  is  reduced. 
However,  it  should  be  noted  that  the  wing  areas  of  the  6F  and  6GA  con¬ 
figurations  approximately  equal  the  minimum  wing  areas  for  which  flight 
at  40,000  ft  is  possible. 

In  summary,  the  joined  wing  configuration  has  a  marked  advantage 
in  range  for  missiles  that  are  required  to  cruise  at  high  altitude  and/or 
to  maneuver  above  the  3  g  level  at  low  altitude. 

Volume  Constraints:  As  shown  in  Section  3,  joined  wings  have  less 
volume  than  equivalent  conventional  wings.  This  offers  a  potential  range 
advantage  for  a  cruise  missile  since  the  total  vehicle  volume  is  limited 
by  the  dimensions  of  the  torpedo  tube  from  which  it  is  launched.  The 
range  advantage  could  be  significant  for  a  tactical  type  of  cruise  mis¬ 
sile,  equipped  with  a  large  (non-nuclear)  warhead  which  leaves  relatively 
little  volume  available  for  fuel  tanks.  For  the  strategic  type  of  missile 
considered  here  the  benefits  are  modest,  as  shown  below. 

Assuming  a  30%  ullage  factor,  the  volume  occupied  by  the  folded  wings 
is  subtracted  from  the  total  volume  available  for  wings-plus-fual  tanks 
(23.013  ft^) .  The  remainder  is  the  available  fuel  volume: 


Other  things  being  equal  (e.g.,  L/D,  initial  cruise  weight)  the  range  of 
the  joined  wing  configuration  would  be  greater  than  that  of  the  corre¬ 
sponding  conventional  configuration  in  (approximately)  the  ratio  indi- 
dated.  No  attempt  was  made  to  exploit  this  advantage  in  the  performance 
calculations  presented  above  because  it  would  have  necessitated  raising 
the  initial  cruise  weight  above  the  uniform  value  of  2864  lb  selected 
for  all  configurations.  This  would  introduce  new  factors  (such  as  dif¬ 
ferent  sizes  of  boosters)  which  would  complicate  the  performance  compari¬ 
sons.  However  the  low  volume  of  the  joined  wing  facilitates  the  task  of 
accommodating  payload  and  systems  within  the  vehicle  and  would  be  more 
significant  for  a  short-range,  densely-packed  missile  in  which  the  avail¬ 
able  fuel  volume  is  critical.  For  all  configurations  discussed  here  the 
fuel  volume  was  adequate  to  contain  the  fuel  appropriate  to  the  assumed 
initial  cruise  weight  of  2864  lb. 
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6.0  CONCLUSIONS  AND  RECOMMENDATIONS 


6.1  CONCLUSIONS 

1.  The  joined  wing  concept  has  been  described,  and  its  application 
to  cruise  missile  design  has  been  studied.  For  cruise  missiles  the  joined 
wing  yields  significant  gains  in  range  and  maneuverability.  It  is  also 
simple  to  fold  and  occupies  much  less  volume  inside  the  missile  body  than 
a  comparable  cantilever  wing. 

2.  A  study  has  been  made  of  terrain-following,  comparing  the  accu¬ 
racy  of  joined  wing  and  conventional  cruise  missiles,  each  equipped 
with  optimal  terrain-following  guidance  systems.  The  joined  wing  yields 
closer  terrain-following  provided  that  two  forward-mounted  controls  are 
employed.  This  can  be  done  by  combining  a  joined  wing  with  a  canard  to 
form  a  hybrid  configuration,  having  control  surfaces  on  the  canard  and 
on  the  front  wing. 

3.  If  desired,  the  hybrid  configuration  can  be  staged  by  employing  a 
j ettisonnable  nose  tank  to  which  the  canard  is  attached.  Staging  improves 
range  and  maneuverability  and  also  reduces  the  size  of  the  missile  during 
its  terminal  phase  so  that  it  is  less  easily  detected. 

6.2  RECOMMENDATIONS 


The  joined  wing  is  a  highly  integrated  concept  involving  structures 
and  aerodynamics.  Further  research  in  both  these  areas  is  required,  as 
indicated  below. 

(1)  Wind  tunnel  tests  (subsonic,  transonic  and  supersonic) . 

(2)  Extension  of  monoplane  airfoil  design  methods  to  joined 
wings. 

(3)  Development  of  procedures  for  calculating  induced  drag 

and  trimmed  of  joined  wing  and  hybrid  configura¬ 

tions. 
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(4)  General  structural  analyses  to  determine  the  optimum 
joined  wing  geometry  for  maximum  stiffness  and  minimum 
weight.  In  particular,  studies  should  be  made  of  the 
effect  of  taper  ratio  and  spanwise  location  of  the  wing 
joint  on  stiffness  and  weight. 

(5)  More  detailed  analyses  of  terrain-following,  including 
the  effects  of  control  system  lags  and  alternative 
terrain  spectra. 

(6)  Design  studies  of  joined  wing  alternatives  to  current 
and  projected  aircraft  and  missiles. 


136 


REFERENCES 


1.  Wolkovitch,  J.,  Joined  Wing  Aircraft,  U.S.  Patent  3,942,747,  1976. 

2.  Wolkovitch,  J.,  An  Introduction  to  the  Joined  Wing  Airplane,  EEC 
Rept.  No.  771,  Engel  Engineering  Co.,  28603  Trailriders  Drive, 

Rancho  Palos  Verdes,  CA  90274,  1977. 

3.  Fairchild,  M.  P.,  Structural  Weight  Comparison  of  a  Joined  Wing  and 
a  Conventional  Wing,  M.S.  Thesis  in  Engineering  Mechanics,  Univer¬ 
sity  of  Texas  at  Arlington,  Arlington,  TX,  1980.* 

4.  White,  E.  R. ,  Preliminary  Force  Test  Data  for  a  Joined-Wing  Aircraft 
Conf iguration,  Internal  Memo  File  V-19100/OLTR-094,  NASA  Langley 
Research  Center,  Hampton,  VA,  8  April  1980. 

5.  Robertson,  F.  H.,  et  al . ,  Report  on  Comparison  Between  the  Joined 
Wing  and  Orthodox  Layouts  When  Applied  to  a  Long  Range  V.I.P. 

Transport  Airplane  Design.  Vought  Rept.  No.  2-31200/8R-51546,  Vought 
Corporation,  Dallas,  TX,  1978. 

6.  DeYoung,  J. ,  Fuselage  Structural  Efficiency  of  N-Plane  Compared  with 
Single  Wing  Transport  Aircraft.  Vought  (Hampton  Tech.  Center)  Memo 
V-l 9000/ 5NASA-3 68 ,  Dec.  1975. 

7.  Bottomley,  J.  W. ,  "The  Tandem-Wing  Concept  Applied  to  Modern  Trans¬ 
ports,"  Aeronautical  J.,  Nov.  1974,  pp.  523-524. 

8.  Reid,  E.  G.,  Applied  Wing  Theory,  McGraw-Hill,  1932. 

9.  Igoe,  W.  B.,  R.  J.  Re,  and  M.  D.  Cassetti,  Transonic  Aerodynamic  Char¬ 
acteristics  of  a  Wing-Body  Combination  Having  a  52.5  deg  Sweptback 
Wing  of  Aspect  Ratio  3  with  Conical  Camber  and  Designed  for  a  Mach 
Number  of  Vl ,  NASA  TN  D-817,  May  1961. 

10.  Durand,  W.  F.,  Aerodynamic  Theory,  1936;  reprinted  by  Dover  Press, 

New  York,  1963  (section  by  L.  V.  Kerber) . 

11.  Glauert,  H.,  The  Elements  of  Aerofoil  and  Airscrew  Theory,  Cambridge 
University  Press,  England,  1926. 

12.  Laitone,  E.  V.,  "Positive  Tail  Loads  for  Minimum  Induced  Drag  of 
Subsonic  Aircraft,"  J.  Aircraft.  Vol.  15,  Dec.  1978,  pp.  837-842. 

13.  Miranda,  L.  R. ,  and  G.  L.  Dougherty,  Transonic  Wind  Tunnel  Testing 
of  a  Low  Induced  Drag  Lifting  System,  Vols.  I  and  II.  Contract 
N00019-73-C-0499,  Lockheed  Aircraft  Corp.,  Burbank,  CA,  Feb.  1974. 


*Published  in  a  shortened  form  as  AIAA  Paper  81-0366,  with  the  same 
title. 


137 


I 


14.  Wolkovitch,  J.,  "Subsonic  VSTOL  Aircraft  Configurations  with  Tandem 
Wings,"  J.  Aircraft,  Sept.  1970,  pp.  605-611. 

15.  Hoerner,  S.  F. ,  Fluid -Dynamic  Drag,  Hoerner  Fluid  Dynamics,  Brick 
Town,  NJ,  1965. 

16.  Henderson,  W.  P.,  and  J.  K.  Huffman,  Aerodynamic  Characteristics  of 
a  Tandem  Wing  Configuration  at  a  Mach  Number  of  0.30,  NASA  TM 
X-72779,  Oct.  1975. 

17.  Letcher,  J.  S.,  "V-Wings  and  Diamond-Ring  Wings  of  Minimum  Induced 
Drag,"  J.  Aircraft,  Aug.  1972,  pp.  605-607. 

18.  Kuhlman,  J.  M. ,  and  T.  J.  Ku,  Numerical  Optimization  Techniques 
for  Bound  Circulation  Distribution  for  Minimum  Induced  Drag  of  Non- 
planar  Wings,  NASA  Tech.  Progress  Rept.,  Grant  NSG1357,  Nov.  1979. 

19.  Margason,  R.  J.,  and  J.  E.  Lamar,  Vortex -Lattice  Fortran  Program 
for  Estimating  Subsonic  Aerodynamic  Characteristics  of  Complex 
Planforms,  NASA  TN  D-6142,  1971. 

20.  Lamar,  J.  E.,  A  Vortex -Lattice  Method  for  the  Mean  Camber  Shapes  of 
Trimmed  Noncoplanar  Planforms  with  Minimum  Vortex  Drag,  NASA  TN  D-8090, 
1976. 

21.  McCormick,  B.  W.,  Aerodynamics,  Aeronautics,  and  Flight  Mechanics, 
Wiley,  New  York,  1979. 

22.  Ashley,  H.,  Engineering  Analysis  of  Flight  Vehicles,  Add ison -Wesley , 
Reading,  MA,  1974. 

23.  Hoerner,  S.  F. ,  and  H.  V.  Borst,  Fluid -Dynamic  Lift,  Hoerner  Fluid 
Dynamics,  Brick  Town,  NY,  1975. 

24.  Torenbeek,  E.,  Synthesis  of  Subsonic  Airplane  Design,  Delft  Univer¬ 
sity  Press,  Holland,  1976. 

25.  Lacey,  D.  W.,  Aerodynamic  Characteristics  of  the  Close-Coupled  Canard 
as  Applied  to  Low-to -Moderate  Swept  Wings.  Vol.  2:  Subsonic  Speed 
Reg ime ,  David  W.  Taylor  Naval  Ship  Research  and  Development  Center, 
Bethesda,  MD,  Rept.  DTNSRDC-79/002,  1979. 

26.  Addoms,  R.  B.,  and  F.  W.  Spaid,  "Aerodynamic  Design  of  High  Performance 
Biplane  Wings,"  J.  Aircraft.  Vol.  12,  No.  8,  Aug.  1975,  pp.  629-630. 

27.  Smith,  A.  M.  0.,  "High-Lift  Aerodynamics,"  J.  Aircraft.  June  1975, 
pp.  501-530. 

28.  Lange,  R.  H.,  et  al . ,  Feasibility  Study  of  the  Transonic  Biplane  Con¬ 
cept  for  Transport  Aircraft  Application.  NASA  CR-132462,  June  1974. 


138 


29.  Zimmer ,  H.,  Das  Hochauf tr iebsverhalten  Beim  Rautenf lugelkonzept , 
Dornier  GMBH,  Friedrichshafen,  Germany,  DGLR  Rept.  No.  78-114,  1978. 

30.  Cahill,  J.  F.,  and  D.  H.  Stead,  Preliminary  Investigation  at  Subsonic 
and  Transonic  Speeds  of  Aerodynamic  Characteristics  of  a  Biplane 
Composed  of  a  Swept-Back  and  Swept-Forward  Wing  Joined  at  the  Tip, 
NACA  RM  L53L34B,  Mar.  1954. 

31.  Henderson,  J.  S.,  Wind  Tunnel  Tests  of  a  Model  of  a  Warren-Young  Wing 
Aircraft ,  unpublished  report,  Department  of  Engineering,  University 
of  Bath,  England,  1976. 

32.  Nicolai,  L.  M. ,  ”A  Perspective  on  the  Requirements  and  Design  for 
Advanced  Cruise  Missiles,"  AIAA  Paper  79-1817,  1979. 

33.  Cronvich,  L.  L.,  and  H.  P.  Liepman,  Advanced  Missile  Technology: 

A  Review  of  Technology  Improvement  Areas  for  Cruise  Missiles ,  NASA 
CR-3187 ,  1979. 

34.  Washington,  W.  D. ,  R.  E.  Wittmeyer,  and  W.  H.  Appich,  Jr.,  "Body  Slot 
Effects  on  Wing-Body  and  Wing-Tail  Interference  of  a  Typical  Cannon- 
Launched  Projectile,"  AIAA  Paper  80-0260,  1980. 

35.  Levinsky,  E.  S.,  et  al. ,  Semispan  Wind  Tunnel  Test  Evaluation  of  a 
Computer-Controlled  Variable  Geometry  Wing  (Self-Optimizing  Flexible 
Technology  Wing  Program) ,  AEDC-TR-78-51 ,  1978. 

36.  Magdaleno,  R.  E.,  and  J.  Wolkovitch,  Performance  Criteria  for  Linear 
Constant-Coefficient  Systems  with  Random  Inputs,  USAF  ASD  TR  62-470, 
ASTIA  297805,  1962. 

37.  Cunningham,  E.  P.,  "Probability  of  Crashing  for  a  Terrain-Following 
Missile,"  J.  Spacecraft  and  Rockets,  Vol.  11,  No.  4,  April  1974, 
pp.  257-260. 

38.  Cunningham,  E.  P.,  "Probability  of  Crashing  from  Monte  Carlo  Simula¬ 
tion,  "  J_I_S£acecraft_and_Rockets ,  Vol.  16,  No.  5,  Sept. -Oct.  1979, 
pp.  348-350. 

39.  McRuer,  D.,  I.  Ashkenas,  and  D.  Graham,  Aircraft  Dynamics  and  Auto¬ 
matic  Control,  Princeton  University  Press,  Princeton,  N.J.,  1973. 

40.  Kwakernaak,  H.,  and  R.  Sivan,  "The  Maximally  Achievable  Accuracy  of 
Linear  Optimal  Regulators  and  Linear  Optimal  Filters,"  IEEE  Trans. 

A.C . ,  Vol.  AC-17,  Feb.  1972,  pp .  79-86. 

41.  Cunningham,  E.  P.,  "Single-Parameter  Control  for  Terrain  Following," 
To  be  published  in  J.  Aircraft.  Dec.  1980. 

42.  Etkin,  B.,  "The  Turbulent  Wind  and  Its  Effect  on  Flight,"  AIAA 
Paper  80-1836,  1980. 


139 


43.  Abbott,  I.  H.,  and  A.  E.  von  Doehnhoff,  Theory  of  Wing  Sections, 

Dover,  N.  Y.,  1960. 

44.  Linden,  J.  E.,  and  F.  J.  O'Brimski,  "Some  Procedures  for  Use  in  Per¬ 
formance  Prediction  of  Proposed  Aircraft  Designs,"  SAE  Paper  650800, 
1965. 

45.  Sanders,  M.  E.,  and  D.  E.  Wieduwilt,  Performance,  Stability  and 
Control  Characteristics  of  the  Full-Scale  Navy/Vought  Tomahawk  Sea- 
Launched  Cruise  Missile  (YGBM-TIO)  at  Subsonic  Mach  Numbers,  AEDC- 
TR-77-18,  1977  (UNCLASSIFIED). 

46.  Larrabee,  E.  E.,  "The  Aerodynamic  Design  of  Sailplane  Tail  Assemblies, 
Technical  Soaring,  Vol.  V,  No.  1,  1977. 


140 


APPENDIX  A 


CALCULATION  OF  OPTIMUM  TERRAIN-FOLLOWING  PERFORMANCE 

Al.  GENERAL  OBJECTIVES 

Figure  Al  shows  a  cruise  missile  flying  above  terrain.  As  discussed 
in  Che  main  text  the  missile  is  assumed  to  have  a  downward -looking  terrain¬ 
following  system,  schematically  illustrated  in  Fig.  A2 .  The  purpose  of 
this  system  is  to  maintain  the  missile  close  to  a  desired  clearance  above 
the  terrain.  This  desired  clearance  (or  height  bias)  is  denoted  as  hg. 

We  are  interested  in  the  question  of  what  is  the  minimum  value  of  hg 
that  can  be  employed.  This  depends  upon  the  statistical  variation  of  the 
height  error  he,  which  at  any  instant  is  given  by: 

he  =  h  -  (hw  +  hB) 

whence 

hw  *  Wave  (or  terrain)  inertial  height, 

h  «  Inertial  height  of  some  specified 
point  on  the  airframe 

If  we  assume  that  the  atmospheric  turbulence  and  the  wave  or  terrain  height 
have  Gaussian  distributions,  then  to  minimize  hg  for  a  given  probability 
of  clobber,  pc>  the  control  system  should  be  designed  to  minimize  the  aver¬ 
age  value  of  h|.  Applying  the  ergodic  hypothesis  reduces  the  problem  to  a 

standard  type  of  optimization  problem,  i.e.,  to  find  the  control  system 

2 

that  minimizes  E(he),  where  E  denotes  expected  value. 

Chang  (Ref.  Al)  and  numerous  other  references  describe  how  such  a 
control  system  can  be  found  by  application  of  optimal  control  theory.  Two 
alternate  types  of  solution  techniques  are  available: 

1.  Time-domain  methods  (popularized  by  Kalman, 
et  al . ) 

2.  Frequency-domain  methods  (popularized  by 
Weiner,  et^  al . ) 


1 


FIGURE  Al.  CRUISE  MISSILE  FLYING  ABOVE  TERRAIN 


FIGURE  A2.  TERRAIN-FOLLOWING  SYSTEM  BLOCK  DIAGRAM 


A2 


For  our  purposes  the  frequency-domain  approach  is  preferable,  since 
it  clearly  identifies  the  airframe  characteristics  which  (for  some  con¬ 
figurations)  limit  the  closeness  with  which  the  airframe  can  follow  the 
terrain.  The  frequency  domain  method  is  admissible  since  it  may  be  assumed 
that  the  statistical  properties  of  the  signal  (terrain)  and  noise  (atmos¬ 
pheric  turbulence)  are  both  stationary. 

This  appendix  presents  the  details  of  the  optimization  calculation. 

It  should  be  noted  that  our  principal  interest  is  in  the  performance  of  the 

optimum  system  rather  than  how  it  is  mechanized.  Thus  we  are  primarily 

y  9 

concerned  with  E(hg)  and  the  mean  square  control  deflection  E(6  )  rather 

than  the  form  of  the  compensation  transfer  functions  Gc  and  H  in  Fig.  A2. 

In  the  following  sections  we  discuss: 

1.  The  mathematical  model  of  the  airframe. 

2.  The  mathematical  model  of  the  "terrain" 

(sea  surface) . 

3.  The  mathematical  model  of  atmospheric 
turbulence. 

4.  Solution  for  the  optimum  system. 

5.  Calculation  of  the  performance  of  the 
optimum  system. 

This  appendix  presents  the  required  mathematical  models  and  theory.  Program 
listings  for  the  HP-41C  computer  are  given  in  Appendix  C.  Results  of  the 
optimization  calculations  are  presented  in  the  main  text  of  this  report. 

A2.  MATHEMATICAL  MODEL  OF  THE  AIRFRAME 

It  is  assumed  that  the  missile  is  equipped  with  a  speed  control  (via 
the  throttle)  such  that  it  flies  over  the  terrain  at  essentially  constant 
inertial  speed.  Thus  the  missile's  longitudinal  dynamics  can  be  approxi¬ 
mated  by  the  short-period  equations  of  motion  (Ref.  A2) .  These  equations 
are  written  in  terms  of  pitch  attitude  (0)  and  angle  of  attack  (a),  however 
it  is  more  convenient  to  rewrite  the  equations  in  terms  of  6  and  height  of 
the  center  of  gravity  (h) ,  where  h  *  /  U0(6  —  a)  dt.  The  desired  form  of 
the  short-period  equations  then  becomes: 


(Al) 


s  Zw 

2a 

-sh 

"  26 

s2  _  (Mq  +  Hd)s  -  Ma 

e 

m6- 

whence 

JL  _  ~Z(SS2  +  Zg(Mq  +  Mg)s  +  ZgMg  -  M6Za 
5  s2 [ s2  -  (Mq  +  +  Zw)s  ~Ha+  MqZw] 

J3_  _  (ZgM^  +  M5)s  +  Z^  -  MgZy 

5  s [ s 2  -  (Mq  +  +  Zw)s  -  Mq  +  MqZj 


(A2) 


(A3) 


Also  of  interest  are  the  transfer  functions  giving  the  response  of  the  air 
frame  to  vertical  gusts  of  velocity  Wg  fps.  These  transfer  functions  are 
formed  following  the  procedure  of  Ref.  A2,  page  251,  which  accounts  for 
gust  gradient  effects  by  replacing  the  right  side  of  Eq.  27  by 


— Zw  +  Zqs/U0 
-Mw  +  Mqs/U0  -  M*s  J 


Since  Zq  «  U0  it  may  be  neglected  over  the  frequency  range  of  interest 
here,  whence,  after  some  reduction: 


h  _  ^y(s  2Mg ) 

wg  s^sp 

where 

*sp  -  s2  "  <Mq  +  Md  +  Vs  -  Ma  +  MqZw 


(A4) 


A4 


i 


S2(Mq  ~  Md)  +  sC-Ma  ~  MqZw) 
sAc 


-»sp 


(A5) 


A3.  WAVE  SPECTRA 

Fortenbaugh  (Ref.  A3)  discusses  various  wave  power  spectra  for  short- 
crested  (confused)  seas  and  recommends  the  Bretschneider  spectrum,  which 
has  the  equation 


where 


Sw  (w) 


u)  =  Wave  frequency  in  rad/sec 

A  =  483.5  H^/T^ 

B  *  1944. 5/T4 
o 

TQ  =  Modal  wave  period  (sec) 

Hs  *  Significant  wave  height  (ft) 


(A6) 


The  parameters  T0  and  Hg  are  correlated  with  sea  state  and  wind  condition 
in  Ref.  A3,  page  19.  We  shall  assume  a  Sea  State  5  with  a  20  knot  head¬ 
wind  as  a  typical  condition.  This  gives  Hs  =  12  ft,  T0  =  12.07  sec,  whence 


3.28  -0.09162/uf 


“w 


U) 


(A7) 


For  a  moving  observer  the  encountered  wave  frequency,  ue,  is  related 
to  w  by 


ui 


e 


ojV 

g 


cos 


(A8) 
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where 

V  =  Inertial  speed  (fps) 

ys  =  The  observer's  heading  with  respect 
to  the  predominant  wave  direction 

For  a  cruise  missile  flying  at  airspeed  corresponding  to  M  =  0.7 
directly  into  a  20  knot  headwind  the  inertial  speed  is  748.12  fps.  For 
this  condition  the  Bretschneider  spectrum  is  shifted  in  frequency,  as  shown 
in  Fig.  A3.  The  amount  of  the  shift  varies  with  the  frequency,  hence  the 
shape  of  the  spectrum  changes.  This  is  shown  in  Fig.  A3  which  compares 
Bretschneider  spectra  for  a  stationary  observer  and  for  an  observer  moving 
at  748.12  fps  inertial  speed. 

The  formula  for  the  Bretschneider  spectrum  is  inconvenient  for  opti¬ 
mization  calculations;  for  such  calculations  it  is  desirable  to  employ 
spectra  which  can  be  expressed  as  ratios  and  products  of  linear  factors 
of  the  frequency  variable  joj,  i.e.. 


Power  Spectrum 


HpUj  +  ju>e?  ^  ~  jtae) 

nq(t>i  +  jwe)(bi  “  j^e) 


»  q  >  p 


(A9) 


Figure  A3  shows  such  an  approximation 


we 

Sw  -  2,238,721.0  - - - -  (A10) 

(0)2  +  6.32)6 

e 

For  our  optimization  calculation  Eq.  A10  was  modified  to  increase  the  power 
at  low  frequencies.  This  was  done  to  model  the  long  waves  and  swells  char¬ 
acteristic  of  fully  developed  seas  (as  opposed  to  the  confused  seas  modeled 
by  the  Bretschneider  spectrum),  and  also  to  model  possible  maneuvering 
commands  (e.g.,  for  evasive  action).  Such  commands  are  expected  to  be 


A6 


i  39  69) 


FtCURF  A3.  COMPARATIVE  HAVE  SPECTRA 


largely  in  the  1—10  rad/sec  frequency  range.  Accordingly  the  input  com¬ 
mand  spectrum  $hwhw  emPl°yed  is: 


Equation  23  is  also  graphed  on  Fig.  A3.  E(h^)  equals  81.2  ft^. 

A4 .  ATMOSPHERIC  TURBULENCE 

Only  the  vertical  component  of  turbulence,  w  ,  is  modeled.  The  well- 

- & 

known  Dryden  spectrum  (Ref.  A4)  gives  this  as: 


^3  *  Spatial  frequency  (rad/ft) 

L3  «  Turbulence  scale  length  (ft) 
°wg  *  rms  turbulence  intensity  (fps) 


(A12) 


*Note  that  the  spectrum  convention  employed  in  Ref.  A4  differs  by 
a  factor  of  tt  from  the  definition  used  here.  The  necessary  corrections 
have  been  made  throughout  the  present  report. 


We  shall  assume  *  100  ft.  For  the  missile  flying  parallel  to  the  iner¬ 
tial  wind  direction  at  airspeed  VA  fps,  the  encounter  frequency  in  rad/sec 
is  given  by: 


toe  =  QVa  (A13 

The  vertical  gust  r.m.s.  amplitude  is  assumed  to  be  2.9  fps,  based  on  data 
presented  in  Ref.  A5,  page  427,  whence: 


^w  w  ) 
wgwg 


1.037(1  +  0 . 221518jaie) 

1  +  0.255787juie  +  0.016357(juJe)2 


14.0438(s  +  4.5143) 
(s  +  7.819)2 


(A14; 


In  essence,  the  purpose  of  the  optimization  is  to  find  the  system  that 
will  achieve  the  best  possible  compromise  between  following  the  commands 
and  not  being  excessively  disturbed  by  atmospheric  turbulence. 

AS.  OPTIMAL  CONTROL  SYSTEM 

For  purposes  of  this  report  the  control  system  is  assumed  to  be  noise- 
free,  e.g.,  the  measured  height  of  the  missile  above  the  local  wave  surface 
is  assumed  to  be  known  precisely.  Furthermore,  internal  system  lags  due  to 
digital  effects  and  actuator  lags  are  neglected.  Such  effects  are  impor¬ 
tant  in  designing  real  hardware,  but  can  be  neglected  for  comparative 
studies  of  joined  wing  versus  conventional  airframe  configurations. 

The  performance  index  selected  is 

J  ■=  E(h2)  +  k2E(62)  (A15! 

e 

where 

k  is  a  constant,  discussed  below 
6  ■  Control  deflection  (from  trim  position) 


Only  a  single  (inner-loop)  control  is  considered,  e.g.,  deflection  of  the 
all-moving  tail  on  the  conventional  configuration.  Combined  elevator- 
plus-flap  control  is  not  considered  (see  discussion  in  main  text) .  How¬ 
ever,  it  is  assumed  that  throttle  is  employed  as  an  outer-loop  control  to 
maintain  constant  airspeed,  thus  suppressing  phugoid  motions. 

The  form  of  performance  index  given  in  Eq.  A15  reflects  our  views  on 
what  aspects  of  system  performance  we  consider  to  be  important.  Obviously 
it  is  desired  to  minimize  E(h^),  but  at  the  same  time  it  is  undesirable 
that  this  should  be  accomplished  at  the  cost  of  excessive  control  deflec- 
tions.  The  parameter  k  represents  this  tradeoff,  and  we  shall  obtain 
solutions  for  a  wide  range  of  k^  to  illustrate  its  effects.  As  shown  in 
Ref.  Al,  additional  parameters  can  be  introduced  into  the  performance  index. 
For  example,  we  could  select  an  index  of  the  form 

J  -  E (h2)  +  kgE(62)  +  k2E(nz  -  l)2  (A16) 

where  k^,  kn  are  constants  and  nz  «  normal  acceleration  measured  in  g's. 

Such  an  index  would  direct  penalize  systems  which  require  the  airframe  to 
operate  near  its  maximum  "g"  capability  for  a  large  portion  of  the  flight. 
However,  increasing  the  complexit-.  jf  the  performance  index  introduces 
unwelcome  complications  into  the  optimization  procedure.  It  should  also 
be  realized  that  the  weighting  E(6  )  also  implicitly  weights  E(nz  —  1)  . 
Hence  Eq.  A16  appears  to  be  adequate  for  the  present  study. 

A6 .  DERIVATION  OF  THE  OPTIMAL  CONTROL  SYSTEM 

The  derivation  given  below  draws  upon  concepts  and  techniques  presented 
by  Whitbeck  in  Ref.  A5,  which  should  be  consulted  for  a  more  detailed  dis¬ 
cussion. 

For  the  system  shown  on  Fig.  A 2  the  control  deflection  can  be  written 
in  vector  form  as 


A10 


s 


6  =  [I  +  GC(H  +  HTn)G1]_1GcR  -  [I  +  GC(H  +  Hm)G^]~1  GC(H  +  Hm)G2d  (A17) 


or  more  compactly 


6  =  WaR  +  Wbd 


(A18) 


Equation  A18  defines  Wa  and  Wb . 

The  corresponding  expression  for  the  height  error  vector  he  is 


r  _  G2d  -  Gi[WaR  +  Wbd] 

(A19) 

[1  -  G]Wa]R  -  (G2  +  G1Wb)d 

(A20) 

WjR  +  W2d 

(A21) 

Equation  A21  defines  Wb  and  W2 . 

The  performance  index  is,  by  definition. 


Put 


J  =  E  .trace  hehg  +  k2  trace  6  6T  ^ 


<5  -  [Wa  I  Wb] 


R 

d 


=  WRX 


(A22) 


(A23) 


[(1  -  G]Wa)  I  -(G2  +  GLWb) ] 


gAr1 


(A24) 


[i,  ;  -G2l 


-  ol[va  :  wb] 


gBR1  ■  Gl1®!  (A25) 


Rl»  gA>  gB»  w  are  defined  by  the  above  equations. 


All 


(  X 

J  =  E  R1[Gg  -  GxW]T  (Gb  -  GXW] +  k2R^WTWR1 J 


(A26) 


where  denotes  transpose 


To  form  the  equivalent  expression  for  stationary  random  inputs  R  and  d, 
having  power  spectra  and  4>dd,  we  employ  the  principle  of  the  transient 
analog  (Ref.  A7)  which  for  the  scalar  he  and  6  considered  here  yields  J 
as  the  following  frequency-domain  integral 


J 


where 


jl  r 

2-n-j  I 

— jao 


'Wl^RRWl  +  W2*ddW2 


k2[W. 


*RRWa 


VddVj  <A27> 


=  Terrain  spectrum 
4>dd  *  Turbulence  spectrum 

U)  =  Frequency  variable,  rad/sec 
j  =  & 


and  the  superposed  bars  (Wj,  W2,  Wa,  W^) 
denote  complex  conjugate  quantities 


In  Eq.  A27  the  first  two  terms  give  the  mean  square  height  error,  while 
the  second  two  terms  give  R  times  the  mean  square  control  deflection.  The 
principal  focus  of  the  discussion  of  this  appendix  relates  to  computing 
the  integrals 


and 


1 

2TTj 


(w1$rrWi  +  W24>ddW2)  djU) 


1 

2rj 


(Wa*RRwa  +  VddV 


dju) 
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for  the  condition  when  Wa  and  WB  are  jointly  selected  to  give  the  minimum 
possible  performance  index  J,  for  a  specified  value  of  the  weighting 
parameter  k^. 

To  derive  the  Wiener-Hopf  equation  for  the  optimal  system  first 
differentiate  Eq.  A26  with  respect  to 


ITT  TTT  2TTi 

EWRjGg  -  R^W  HGgR-L  -  G-jWR^  +  k  R^W  WRj  j 


E  JrigJgbR1  -  R^gJg1WR1  -  R^WTGBRi  + 


+  k2R^WTWR1| 


(A28) 


(A29) 


*  e|-R^(Gb  -  G1W)R1  +  k2WRxR^  I 


(A30) 

(A31) 


3J 


awT 


E]-[Gi(Gb  -  GjW)R  ]Ri  +  k2WRxRi  | 


(A32) 


For  the  variation  to  vanish  we  must  have  the  standard  Wiener-Hopf 
condition: 

in  rr  m  m  j 

<P  (A33) 


E  |[k^  +  G^GjJWR-j^R^  —  G^GjjR^rJ  | 


where  is  a  vector  that  exists  only  for  negative  time,  as  explained  in 
Ref.  A6. 

Substituting  G  ■  [1  !  — G2 ] ,  W  =  [Wa  !  WB]  yields: 


eJir2  +  G^G1][wa  ;  Wb] 


RRt  RdT 
dRT  ddT 


*  Gn[l  J  -G2] 


RRT  RdT 


dR1  dd 


(A34) 
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Applying  the  principle  of  the  transient  analog  yields  a  corresponding 
spectral  matrix  equation  which  reduces  to  two  scalar  equations  since  the 
wave  "command"  R  and  atmospheric  turbulence  "disturbance"  d  are  uncorre¬ 
lated,  i.e.,  =  0.  The  scalar  equations  are: 


[k2  +  -  Gi4>rr  =  ^ 


(A3  5) 


[k2  +  GiGi]Wb4>dd  +  G1G24>dd  =  i|>2 


(A36) 


where  ipo  are  functions  of  s  that  have  only  right  half  plane  poles. 

Both  Eqs.  A35  and  A36  have  the  form 

FFWGG  +  N  =  ip  (A37) 


which  has  the  standard  Wiener-Hopf  solution  (Ref.  A6) : 


W  « 


(F)-1  N(G) 


(A38) 


Equation  A38  is  applied  below  to  write  solutions  for  the  specific  Wa  and 
of  interest  here. 

A7 .  SPECIFIC  OPTIMAL  CONTROL  SYSTEMS 

For  the  system  considered  here  we  have: 


818. 9s' 


(s  +  6.3)4  ’ 


s2D 


<*dd> 


14.0438(s + 4. 5143) 
(s  +  7 . 819) 2 


,  G?  ‘  ~ 
’  1  sD 
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Substituting  in  Eq.  A36  yields  after  reduction: 


s2P  n1sZ  (818.9)  (s  +  6.3)' 

k2A  A(s  +  6.3)4  (818. 9)s2 


(s  +  6.3) 


k2A  A(s  +  6.3)4 


whence 


wa  = 


818.9D  Nls2 
k2A  [a(s  +  6.3)4J  + 


s2D  NiN2(s  +  4.5143)  (s  +  7.819)2 


Ak2[sDA(s  +  7.819)2  J+  (s  +  4.5143) 


Now  let  us  define  two  polynomials  in  s,  La  and  L^,  as  follows: 


A(s  +  6.3)“ J+  A(s  +  6.3)' 


N1N2(s  +  4.5143) 
sDA(s  +  7 . 819) 2 


sDA(s  +  7.819)' 


(A39) 


(A40) 


(A41) 


(A42) 


(A43) 


(A44) 


A15 


1 


The  solutions  for  Wa  and  can  then  be  written  as: 


Wa 


(AA5) 


«b 


-sLb _ 

k2A(s  +  A.51A3) 


(AA6) 


The  method  of  solution  for  La  and  Lj,  will  now  be  described.  It  is  rela¬ 
tively  lengthy;  this  is  a  consequence  of  the  repeated  factors  in  (J)^  and 
4>dd>  which  complicate  the  task  of  forming  partial  fraction  expansions. 

A8.  CALCULATION  OF  RESPONSE  OF  OPTIMUM  SYSTEM 


In  this  subsection  we  outline  the  equations  employed  to  calculate 
the  mean  square  values  of  h|  and  5^  from  Eq.  A27  (which  see).  Appendix  C 
presents  computer  programs  for  solving  these  equations.  The  programs  are 
written  for  the  HP  A1C  calculator  system.  For  ease  of  reference  to  these 
programs  the  notation  used  in  these  sections  indicates  the  registers  in 
which  parameters  are  stored,  e.g.,  the  h/6  transfer  function  (Eq.  A2)  is 
written  as: 


R01s2  +  RQ2s  +  R03 
s2(s2  +  R^s  +  Rq5) 


(AA7) 


This  indicates  that  in  Eq.  A2  the  quantity  -Z$  is  stored  in  Register 
No.  .01,  while  Z^Mq  +  M^)  is  stored  in  Register  No.  02,  etc.,  etc. 
Because  of  the  limited  memory  of  the  HP-A1C  it  has  been  necessary  to 
reassign  some  registers  during  the  program.  Hence,  where  Rq^  -  Rqj  are 
reassigned,  h/<5  is  correspondingly  reassigned  as  follows: 


_h 
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R33S2  +  &32s  +  ^33 
s^(s2  +  R34S  +  **35) 


(AA8) 


A16 


It  will  be  clear  from  the  context  whether  Eq .  A47  or  Eq .  A48  is  being 
employed . 


The  gust  response  transfer  function  (Eq.  A4)  is  written  as: 


_h_ 

w_ 


=  G2 


r42s  +  r43 


s(s2  +  R^s  +  R-^) 


35; 


(A49) 
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A. 8.1  Root  Square  Locus  (Programs  "OCT",  "POLY", 

"Zfl/N",  and  "QR" 

The  first  step  in  the  optimization  is  to  calculate  the  variation  of 
the  poles  of  Y  (Eq.  A50)  as  a  function  of  1/k2. 


Y 


(A50) 


where  and  D-^  are  respectively  the  numerator  and  denominator  of  Eq .  A2. 
Rearranging  Y  give' 


(A51) 


Chang  (Ref.  Al)  describes  a  graphical  procedure  analogous  to  root  locus 
construction  for  finding  the  roots  of  YY  as  functions  of  1/k^.  This  pro¬ 
cedure  is  known  as  the  "root  square  locus"  method.  The  graphical  aspects 
of  the  method  are  not  of  interest  here,  but  the  mathematical  methods  are 
employed  to  find  the  poles  of  Y.  We  put 

M  =  sADD  +  (l/k2)N1N1  (A52) 

*  A^  +  Rq^A  +  RQ2A  +  Rq^A  +  Rqq 

where  A  =  s2. 


A17 
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The  calculation  procedure  is  as  follows: 

1.  AA  is  calculated  (by  Program  "OCT"). 

2.  The  four  X-roots  of  AA  are  found  (by  Program  "POLY"). 

3.  The  eight  s-roots  of  AA  are  found  (by  Program  "Z+l/N") . 

4.  The  four  s-roots  with  L.H.P.  real  parts  are  used  to 
form  factors  of  A  (by  Program  "QR") . 

5.  The  factors  are  multiplied  to  give  A  in  polynomial 
form: 

A  **  s4  +  Rq^s^  +  Rq^s^  +  RQgS  +  Rq9  (A54) 

where  Rgg  -*•  R^g  are  given  as  output  of  Program  "QR" . 


A. 8. 2  Solution  for  Wa  (Program  "AB") 

Equation  A43  can  be  generalized  to: 


Ni  s^ 


A(s  +  6.3r 


Ja 


(s  +  6.3)4  & 


(A55) 


B 


(s  +  6.3)4  (s  +  6.3)- 


(s  +  6.3)" 


(s  +  6.3) 


+  ~  (A56) 


We  require  expressions  for  A,  B,  C,  and  D.  These  expressions  should  pre¬ 
ferably  be  recursive,  i.e.,  B  =  f (A) ,  C  *  f(A,B),  D  =  f(A,B,C)  to  reduce 
the  amount  of  computation  required. 

4 

A  can  be  found  by  multiplying  both  sides  of  Eq.  A56  by  (s  +  6.3)  and 
then  putting  s  *  -6.3,  yielding 


A 


A 


s»-6 . 3 


R14R00 

R10 


R17 


(A37) 


A18 
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The  symbols  Ri4(=  N^)  ,  R00(®  6.3),  R-io(*  A),  Rtt,  indicate  storage 

s  ““  o  •  j 

register  allocations  in  the  program  "AB"  described  later.  The  coeffi¬ 
cients  B,  C,  D  are  found  as  follows. 

Rearranging  Eq.  A56  yields 

Nxs2  -  AA  =  BA(s  +  6.3)  +  CA(s4-6.3)2  4-  DA(s  +  6.3)3  +  Ra(s  +  6.3)4  (A58) 

Differentiating  with  respect  to  s 

2NlS  +  s2n[  -  AA'  =  BA' (s  4-6.3)  4-  BA 

4-  2CA(s  4-6.3)  4-  CA'  (s  4-  6 . 3)3 

4-  3CA(s4-6.3)2  4-  DA' (s  4- 6. 3) 3 

4-  4Ra(s  4-  6 . 3) 3  (A59) 

Putting  s  ■  -6.3  yields  the  desired  result 

B  =  (2Nls  4-  s-N^  -  AA’ )  /  A  (A60) 

Differentiating  Eq.  A59  with  respect  to  s  and  putting  s  *  -6.3  yields: 

2N,  4-  4sNi  4-  s2n;  -  AA"  -  2BA’ 

C  -  — ± -  (A61) 

2A 

Similarly,  differentiating  Eq.  A59  twice  yields,  with  s  =  -6.3 

6Ni  4-  6sNV  -  AA'"  -  3BA"  -  6CA' 

D  -  — - - - - - -  (A62) 

6A 

The  storage  register  allocations  for  Eqs.  A60,  A61,  and  A62  are  as 
indicated  below. 


A19 


B 


(A63) 


t 


=  r18 

C  =  R16 

D  =  Rj5 

where  the  asterisk  indicates  that  register  R]j  is  filled  twice  (first  with 

_ I 

N^,  which  is  subsequently  cleared  and  replaced  with  D) . 

Having  obtained  A,  B,  C,  D,  as  described  above,  these  are  now  com¬ 
bined  to  form  La  in  Eq.  A55. 

La  =  A  +  B(s  +  6.3)  +  C(s  +  6 . 3) 2  +  D(s  +  6.3)3 

The  storage  register  assignments  require  clearing  Rqq  (formerly 
Rqo  =  s  »  -6.3)  and  putting  Roo  *  -s  *  6.3,  giving 

La  -  Ri5  (s3  +  3RoQs2  +  3r£qs  +  R^0)  +  R16(s2  +  2^3  +  R^0) 

+  R1g(S  +  Rqq)  +  R^y  (A66) 

The  coefficients  of  powers  of  s  in  Eq.  A66  are  evaluated  and  stored  as 
follows. 

La  *  RjiqS^  +  R^^s2  Ri2s  ^  Ri3  (A67) 


(2RqqRi4  +  Rqori5  “  r17R11^/R: 


10 


2r14  +  4r0Qr15  +  2R00R01  ~  2r18r11  ~  R17R12 

2R10 


(A64) 


6 (r15  +  2R00R0l)  ~  r17r13  ~  3R18R12  ~  6R16R11 

6R10 


(A65) 


A. 8. 3  Solution  for  (4>RR)+Wa 

The  spectrum  of  that  portion  of  the  control  deflection  that  is 
correlated  with  the  terrain  is  given  by  4>RRwaWa.  Integrating  this  spec¬ 
trum  over  the  complete  frequency  range  gives  the  appropriate  component 


A20 


mt. 


of  mean  square  control  deflection.  The  integration  is  most  easily  per¬ 
formed  by  splitting  the  spectrum  into  left-half  plane  and  right-half  plane 
factors,  thus: 


^RRWaWa  =  Wa  x  ^RR^  Wa 

This  subsection  describes  the  evaluation  of  the  left-half  plane  factor 
^Rr)  Wa' 

From  Section  A6  we  have 

Wa 


(818.9/k2)s2DL£ 
A(s  +  6.3)4 


(A68) 


To  evaluate  the  denominator  of  Eq.  A68  put 

(s  +  6.3)4  -  s4  +  25.2s3  +  238.14s2  +  1,000.188s  +  1,575.2961  (A69) 

=  s4  +  R2Qs3  +  R2]s2  +  R22s  +  R23  (A70) 

Multiplying  the  right-hand  side  of  Eq.  A70  by  the  right-hand  side  of  A54 
and  assigning  appropriate  registers  yields: 

A(s  +  6.3)4  -  s8  +  R24s7  +  R25s6  +  R26s5  +  R2ys4 

+  R2gs3  +  R2gs2  +  R30s  +  R31  (A71) 

To  evaluate  the  numerator  of  Eq .  A68  put 


A21 


Reference  A8  gives  formulas  for  evaluating  E(p2)  in  terms  of  the 
coefficients  of  A(s)  and  B(s).  For  the  8th  order  denominators  of  interest 
here  these  formulas  proved  to  be  ill-conditioned.  Attempts  to  evaluate 
the  Ig  formula  of  Ref.  A8  were  made  using  first  the  HP-41C  and  subse¬ 
quently  the  CDC  6600  computers.  Neither  computer  was  able  to  give 


A22 


I 

i 

I 

I 


reliable  results.  For  example,  negative  values  of  E(p2)  were  obtained. 
The  problem  may  be  associated  with  the  repeated  denominator  roots  of  the 
command  and  disturbance  spectra. 

After  repeated  attempts  to  employ  the  formulas  of  Ref.  A8  proved 
fruitless  it  was  decided  to  approximate  Eq.  A75  by: 


E(p2) 


_1_ 

TTj 


1_ 

■n 


A(.1m)A(-ju)) 

B(joOB(-jU)) 


A(jui)A(-jfaj) 

B(jU))B(-jw) 


djto 


do) 


u>  in  rad/sec 


(A7i 


(A7 


As  will  be  described.  Program  "I8N"  is  employed  to  compute  Wa4>RRWa  as  a 
function  of  u).  The  integration  was  then  performed  by  a  standard  HP-41C 
math  module  routine  "INTG".  An  interval  (Aoj)  of  1  rad/sec  was  employed, 
giving  100-step  integration:  no  significant  gain  in  accuracy  was  obtained 
by  increasing  the  number  of  steps.  Some  checks  were  made  regarding  the 
effect  of  the  finite  upper  limit  in  Eq.  A77.  It  was  found  that  the  value 
of  the  integral  typically  changed  less  than  1 %  if  the  upper  limit  of 
integration  was  raised  to  200  rad/sec. 

The  "I8N"  program  employs  nested  expressions  to  reduce  computer 
time.  Thus  the  RHS  of  Eq.  A73  is  rewritten  as: 


( (  (-R34<j^  +  ^35)^  “  R-38)u^  +  R40)  +  jw(  ( (“R33W2  +  R35)io2  -  R^)!*)2  +  R39) 

(  (  ( (o)^  —  ^25 ^ R2y)w  —  R29  4"  R31)  ^  ju>(((  — R24ii)^  ^  R2g ) R2 g) w  4*  R3Q ) 

a)  is  stored  in  Rjj,  other  register  assignments  are  obvious  from  the 
listing  of  "I8N"  given  later. 


A23 


The  next  subsection  describes  the  method  used  to  calculate  the  mean 
square  value  of  the  second  part  of  the  control  deflection  of  the 
optimally-controlled  system,  i.e.,  that  part  which  is  correlated  with 
the  atmospheric  turbulence.  The  procedure  follows  similar  procedures 
to  those  demonstrated  above  for  computing  the  mean  square  value  of  the 
component  of  control  deflection  correlated  with  the  command  (terrain) 
input.  It  was  found  that  the  latter  was  dominant  for  all  cases  con¬ 
sidered  here  and,  as  a  first  approximation,  the  turbulence-correlated 
component  could  be  neglected.  In  general  it  will  not  be  known  a  priori 
that  one  component  is  negligible,  hence  the  second  component  of  control 
deflection  should  be  calculated  as  described  below. 


A. 8. 5  Solution  for  (ftdd)'*'  WB  (Program  "WBC") 


Equations  A46  and  A44  yield: 

(14.0438/k2) (-sLb) 


(<t>dd)+wB 


A(s+  7.819)2 


where 


Lb  -  sDA(s  +  7.819)  ‘ 


N1N1(s  + 4.5143) 


L  sDA(s  +  7.819)^  J 


(A80) 


(A90) 


To  calculate  the  [  ]+  term,  first  expand  its  kernel  thus: 


* 


N1N2(s  +  4.54143) 
sD(s  +  7.819)2 A 


4  3  2 

L4S  +  L3S  +  L£S  +  L^s  +  Lq 

sD(s  +  7.819)2 


-L^s3  +  B2S2  +  B^s  +  Bg 


(A91) 


The  -L4  term  in  Eq.  A91  arises  because  the  LHS  of  Eq.  A91  comprises  a 
4th  order  numerator  whereas  the  denominator  is  9th  order.  Thus  Eq.  A91 
can  be  written  as: 


NjN2(s  +  4.54143)  *  (L4s4  +  L3S3  +  L2s2  +  LjS  +  Lq)A 

+  (-LAs3  +  B2s2  +  B-jS  +  Bq)sD(s  +  7. 819) 2  (A92) 

with  storage  register  assignments  as  follows: 

N2  *  R42s  +  R43  »  N1  =  **01 s2  “  R02S  +  R03  ’  D  “  s2  +  R04s  +  R05  ’ 

(s  +  7 . 819) 2  =  s2  +  R10s  +  R11 


Program  "WBC"  solves  Eq.  A92  for  L^,  L3,  L2,  L^,  Lq  by  equating 
coefficients  of  equal  powers  of  s.  The  register  assignments  are: 


sD(s  +  7. 819)  2 


(A93) 


4.5143  *  R18 


(A94) 


The  LHS  of  Eq.  A92  is  expanded  as  follows: 

n1n2  *  (fy)l®2  ~  ^2®  +  ^3)  (R42®  +  r43^ 

3  2 

■  ^14s  +  ^15s  +^16s  +  r17 


(A95) 

(A96) 


A25 


(A97) 


(s  +  4 . 5143)  ■  +  R2qS^  +  +  R]_^s  + 


The  RHS  of  Eq.  A92  is  expanded  making  use  of  Eq.  A93  to  yield  a  polynomial 
formed  by  the  sum  of  the  terms  listed  below: 


s® 

S7 

s® 

S5 

s* 

s2 

s 

L4 

l3 

l2 

Ll 

L0 

-L4 

"r06l4 

_r06L3 

-R06L2 

"RObLl 

-Roeki 

(1*0) 

“r12L4 

R07L4 

Rq7l3 

R07L2 

R07Ll 

R07LC 

~r13l4 

~r08l4 

-R08l2 

-RosLi 

~R08l0 

r12b2 

-r10l4 

r09l4 

RQ9L3 

Rq9l2 

R{)9l1 

B1 

r13b2 

-R11L4 

r11b2 

R11B1 

R11B0 

r12b1 

r10b2 

R10B1 

R10B0 

B0 

R13B1 

R13B0 

r12b0 


Evaluating  this  polynomial  to  the  RHS  of  Eq.  A96  yields  immediately 

R09L0  "  r17 


(A98) 


L0  ■  M  ■  1115 


(A99) 


The  remaining  terms  L4  ...  Lj_,  B2  ...  Bq  are  given  by  the  following  matrix 
equation. 


t 


The  first  step  in  the  solution  is  to  reassign  the  RHS  elements  of  Eq.  A100 
thus : 

RHS  *  [0,  0,  0,  Ri4»  J  r20»  r15»  r16^  (A101) 

The  most  straightforward  way  to  solve  Eqs.  A100  and  A101  is  to  invert  the 
matrix  by  use  of  the  HP-41C  math  module  "SIMEQ"  routine.  This  routine 
was  not  available  at  the  time  when  program  "WBC"  was  written,  hence 
Eq.  A101  was  solved  by  a  partitioning  method  which  is  briefly  summarized 
below,  for  the  benefit  of  readers  who  wish  to  check  the  listing  of  "WBC" 
given  later.  Other  readers  should  skip  this,  and  go  directly  to  Eq.  A114. 

1.  The  matrix  is  partitioned  as  indicated,  giving  with  obvious 
notation 


whence 


'*11 

Uill 

V 

-Yr 

m  _  _  __ 

.*21 

!  ?22  J 

B  . 

.Y2. 

*22® 

*21Lb 

“  Y2 

(A102) 


(A103) 


A27 


2.  P22  is  calculated  algebraically 


1/RH 

-Rlo/Rll 

(R10  “  r11r13)/r11 

r17 

r18 

r21 

P22  * 

0 

1/R11 

-r10/rii 

s 

0 

r17 

r18 

0 

0 

1/R11 

0 

0 

R17 

(A104) 


3.  Since 

B  -  P^2|Y2  -  p2lM 


B  can  be  expressed,  after  reassignment  of  many  registers,  as: 


[B] 


b2 

r20 

0  R22 

r23 

r21 

B1 

* 

r15 

- 

0  0 

r22 

r23 

®0 

r16 

0  0 

0 

r23 

L4 

l3 

L2 

L1 


4.  From  Eq.  A102 


(A105) 


(A106) 


p12B  +  PllLb  *  Y1 


(A107) 


Substitute  Eq.  A107  into  Eq.  A106;  reassign  the  first  column  of 
P11  as  fR24>  p25>  r26>  r27^»  and  multiply  out.  After  reassign¬ 
ment  this  yields: 


1 

1 

i 

j 

i 


L. 


r24 

r28 

r32 

r36 

r25 

r29 

r33 

r37 

r26 

r30 

R34 

r38 

r27 

r31 

r35 

r39 

T.4 

r13 

l3 

e 

r12 

l2 

R11 

Ll 

- 

R10 

(A108) 


A28 
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5.  Solve  for  [L4,  L3 ] T  in  terms  of  [L2,  Li]t  by  partitioning  and 
inverting  the  resulting  2><2  submatrices,  yielding  aftei 
reassignment: 


f  1 

L4 

r16 

1 

H4 

ho 

1 

R14 

r  *1 

l2 

L3 

r17 

r13 

R15 

li 

Substituting  Eq.  A109  in  Eq.  A108  yields  after  reduction: 


(A109) 


R34  -  R26r12  -  r30r13 

r38  ”  r26r14  “  r30r15 

l2 

R11 ~  r26r16  +  r38r17 

r35  -  r27r12  “  r31r13 

r39 ~  R27R14  “  R31R15 

L1 

R10 ~  R27R16 ~  R39R17 

(A110) 


After  reassigning  registers  this  becomes 


r20 

r22 

l2 

= 

R11 

r21 

r23 

Ll 

R10 

(Alll) 


Solving  Eq.  Alll 
yields,  from  Eq. 


l3 


and  reassigning  Lq  =  R^g , 
A109 

*  R16  ~  R12R11  ~  R14R10 

*  r17  "  r13r11  ”  r15r10 


*  Rll*  L1  *  Rio 

R13 

(A112) 

r12 

(A113) 

whence  the  desired  result: 


Lb  -  L4S4  +  L3S3  +  L2S2  +  L^s  +  Lq  *  Rus^  +  R-^s3  +  Rlls2  +  R10s  +  R14 

(A114) 


A29 
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The  next  subsection  describes  the  procedure  employed  to  calculate  the 
mean  square  value  of  that  component  of  mean  square  height  error  corre¬ 
lated  with  the  terrain.  For  the  cases  considered  here  this  proved  to  be 
the  dominant  component;  the  turbulence-correlated  component  was  negligible 
in  comparison. 


A. 8. 6  Solution  for  (4>rr)  Wj_  (Program  "DE") 


Referring  back  to  Eq.  A27  it  can  be  seen  that  the  portion  of  the  mean 
square  height  error  that  is  correlated  with  the  terrain  (wave)  input  is 
determined  by  (4>rr)+W]_'  From  Eq.  A20  this  can  be  written  as: 

(818.9)(s*A 

(<*W+  W1  *  - - - —  (A78) 

A(s  +  6.3)4 

The  first  term  in  the  numerator  parentheses  is: 

s^A  *=  s  +  RggS^  +  Rq^s^  +  Rqqs^  +  Rq£js^ 

The  second  term  is: 

(Rios3  +  Rus2  +  R12s  +  r13)(r01®2  +  *02®  +  *03) 

where 


A30 
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The  complete  numerator  of  Eq.  A78,  including  the  818.9  factor  (stored  as 
R34)  is  given  by  Program  "DE"  as  follows.  The  lead  term  R33  is  zero,  as 
explained  below. 

^33s^  ^34s^  +  ^35s^  +  ^36s^  +  ^37s^  +  R3gs^  +  R39®  +  R40 

R33  is  zero:  this  is  a  dummy  term  required  for  compatibility  with  the 
"I8N"  program.  R33  through  R^g  are  overwritten,  since  these  registers 
were  previously  used  for  the  numerator  of  Wa(<t>j^)  +  .  The  denominator  of 
(^RR)  W1  is  identical  with  that  of  (<t>RR)+Wa  so  registers  24  through  31 
need  not  be  overwritten. 

The  calculation  of  is  performed  using  Program  "I8N"  as  for 

^RRWa^a>  anc*  t^le  subsequent  integration  is  performed  by  Program  "INTG", 
as  previously. 

Having  found  L^,  it  can  be  substituted  in  Eq.  A80,  giving 

-(14.0438/k2)(R1,s5  +  R-.-s'4  +  R,,s3  +  Riris2  +  R,,s) 

W  =  - id - U. - ii - — - i±—  (A114) 

A(s  +  7 . 819; 


The  denominator  of  this  expression  is: 

(s  +  R06s3  +  Rq^s3  +  Rggs  +  ^09^  (s  ^15S  +  ^16^ 

=  s^  +  R2QS3  +  R-21s^  ^22s^  R23s^  +  ^24S  +  ^25  (A115) 

The  numerator  of  Eq.  114  is  written  as 

-(R26S3  +  ^27®^  r28s^  R29®^  R30®  R31^ 


A31 


whence 


-(<W  wb 


^26s^  +  +  RoqS^  +  R-»nS  +  R 


27 


^28 


29 


^O6 


31 


6  5  A  3  2 

s  +  ^20s  +  ^21s  +  ^22s  +  ^23s  +  ^24s  +  ^25 


(A116) 


A. 8. 7  Integration  of  Spectral  Factors  (Program  "ISN") 

Section  A. 8. 6  describes  the  ill-conditioning  problems  encountered 
with  integration  of  Wa0RRWa.  To  avoid  a  recurrence  of  these  difficulties 
in  the  evaluation  of  »  a  numerical  integration  procedure  was 

employed,  similar  to  the  "I8N"  program.  The  program  for  evaluating 
^ddWbWb  33  a  functi°n  of  frequency,  to,  is  labeled  "I6N".  It  represents 
Eq.  A116  as  a  nested  expression. 


(  (R22co^  —  R2g)to^  +  R-j^)  +  j  (  ^26^^  —  ^28^^^  R3Q)to 

-Wdd^  wb  *  - 7 - o - 7 - 7  7 

(((-ur  +  R2l)°jZ  ”  R23)toz  +  R25)  +  j((R20w  “  R22)(i)  +  R24)a) 

(All 7) 

This  is  multiplied  by  its  complex  conjugate  and  integrated  between  limits 
of  (o  =  0  and  <0  =  100  rad/sec  using  the  HP-41C  math  module  subroutine  "INTG" 
with  100  intervals  of  1  rad/sec,  similarly  to  the  integration  of  the  "I8N" 
program  output  described  in  Section  A. 8. 4. 

A. 8. 8  Solution  for  (<j>dd^+^2  (Program  "W2C"/ 

From  Eq .  A20 


W2  -  -(G2  +  GjWb) 


(A120) 


r 


w 


2  = 


R42s  +  R43 


~(Rq1s2  +  Rq2s ->■  Rq3)s  Lb 
,2  /-_2 


s  (sz  +  R04S  +  Rq5)  sz (s^  +  R04s  +  RQ5)k^A(s  +  4.5143) 


(A121) 


(1/k3) (Rq^s^+Rq2S+Rq3) (Ri3s^+^12S^+^HS^+^10S+^14^-^^S+^ -5143) (R42S+R43) 


s  (s2  +  Rq4s  +  Rq3)A(s  +  4. 5143) 


(A122) 


The  second  term  in  the  numerator  of  Eq.  122  is  evaluated  by  assigning 
4.5143  to  R15,  subsequently  reassigning  R^3  as  indicated  below. 


(R42S  +  R43)  (s  +  R15)  ~  R-iys^  +  R^gS  +  R15 


(A123) 


Other  storage  registers  are  reassigned  as  indicated  below. 


(R17s2  +  R16s  +  R^g) (s4  +  RnAs3  +  Rn7s2  +  Rnfts  +  Rnq) 


06b  TR07S  TR086TR09J 


Ri7s  +  R^gSJ  +  R-j^s4  +  R26sJ  +  R27s^  +  R16s  +  R15 


(A124) 


The  complete  numerator  of  Eq.  A122  is  computed  as  follows,  using 
(1/k2)  =  R41 

(R0iS2  +  Rq2S  +  Rq3)  (R]_3S^  +  R12S^  Rlls^  R10s  R14^R41 


-  R17s6  -  Rlgs5  -  R^qS4  -  R26s3  “  R27s2  -  R16s  “  R15 


(A125) 


Rj7s°  +  R]_gSJ  +  Ri9sH  +  R26s  J  +  R27s<i  +  R16s  +  R15 


A3  3 


We  now  have 


Rj^s^  +  R-^gs^  +  R^gs^  +  R2gs^  +  R27S^  +  R^gs  +  R]_g 

w2  *  - 

s(s2  +  R04s  +  RQ5)A(s  +  4.5143) 


(A126) 


It  can  be  shown  that  R]_g  should  be  zero.  The  proof  of  this  is  given 
in  Section  A. 8. 9,  to  avoid  a  digression  at  this  point.  Due  to  roundoff 
errors  the  computed  R]_g  will  not  be  exactly  zero.  The  program  "W2C" 
prints  out  R^5  as  an  indication  of  the  magnitude  of  the  roundoff  errors. 

With  R15  dropped,  (<t>dd)+w2  is  computed  using  Eq.  A14  with  d  =  wg 


<*dd>+ 


14. 0438 (s  +  4.5143) 
(s  +  7.819)2 


(A127) 


W2<W 


14.0438(R^7s5  +  R]_gs^  +  R^s^  +  R^gs2  +  R27S  +  Ri6) 
(s2  +  R04s+R05)(A)(s  +  7.819)2 


(A128) 


The  denominator  factors  (A) (s +  7 . 819) 2  are  already  stored  from  Program 
"WBC"  as: 

(A)  (s  +  7 . 819) 2  *  s^  +  R2qs^  +  R2iS^  +  R22s^  "*■  ^23s^  +  ®^24s  +  ^25 


Multiplying  by  (s2 + Rq4s + Rqj)  yields  the  complete  denominator  as: 

s®  +  R24S2  +  R25S®  +  R2gS^  +  R27S^  +  R2gS2  +  R29S2  +  R^qS  ^31 

The  "W2C"  program  prints  out  R24  through  R4q  to  identify  ^(^dd)"^  in 
the  following  form 


A34 


/■ 


s+tt  R33s7  +  R34s6  +  R35s5  +  R36s4  +  R37s3  +  R38s2  +  R39s  +  R40 

W2  “  q  7  ft  ft  7  o  7 

s°  +  R^s  +  ^25s°  +  R26s  J  +  R27“  +  R28s  J  +  R29s^  +  R3QS  +  R31 


(A129) 


Integration.  Since  (4>8cj)  W2  has  an  eighth-order  numerator  it  can  be 
integrated  by  Program  "I8N"  exactly  as  described  in  Section  A. 8. 4. 


A. 8. 9  Proof  That  s  Is  a  Factor  of  the  Numerator  of  W2 


This  proof  is  required  for  the  simplification  of  Eq.  A126.  From 
Eqs.  A120  and  A46, 


W2  =  -G2  -  GlwB  *  — 


-N2  Nx  s  Lb  (1/k2) 
sD  s2D  A(s  +  4.5143) 


(A130) 


-N2/A(s  +  4.5143)  +  Nx/ Lb  (1/k2) 


srD A  (s  +  4.5143) 


(A131) 


But  from  Eqs.  A91  and  A44, 


Lb  •  A  +  B  •  sD(7.819)  =  N^Cs  +  4.5143) 


(A132) 


where  B  =  -L^s3  +  BjS2  +  B-^s  +  Bq  in  Eq.  A91.  If  s  =  0  is  a  factor  of 


the  numerator  of  W2  then 


[-N2A(s  +  4.5143)  +  N1Lb(l/k^)]s=0  =  0 


(A133) 


But  from  Eq.  132,  for  s  =  0,  I^A  ■  N^N2(s  +  4.5143) 


|i  1  .  .I),IIUIU.^J 


I 


4 


-N2A(s  +  4.5143)  + 


N1N1N2(s  +  4.5143) 


(A134) 


-N2(s  +  4.5143)  /  _  i  _ 

— - Z -  AA  -  ±  NlNl 

A  \  V. 


(A135) 


But  from  Eq.  A52,  by  definition. 


AA  =  s4DD  +  -^  N1N1 
k 


whence  the  term  in  the  inner  parentheses  in  Eq.  A135  becomes 


s*DD  +  N1N1 


0  ,  for  s  *  0 


(A136) 
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SYSTEM  PERFORMANCE  WITH  ZERO  COST  OF  CONTROL 

r 

» i 

B . 1  BASIC  RELATIONSHIPS 

i  ' 

t 

The  performance  index  is: 

r  i 

J  =  E(h|)  +  k2E(62)  (Bl)  ! 

For  a  nonminimum  phase  airframe  Jmin  does  not  tend  to  zero  as  k2  +  0.  The  j 

value  of  E(h2)  achieved  for  k2  +  0  is  of  interest  because  it  is  the  mini¬ 
mum  possible  mean  square  height  error  achievable  with  a  single  control  <5. 

With  the  same  symbols  employed  in  Appendix  A,  the  total  mean  square 
height  error  is  given  by: 


where 


E(hg) 


(W1«>RRW1  +  VddV 


denotes 


d  j  to 


(B2) 


The  mean 

control  square 

control  deflection  is  given  by: 

f* 

(B3) 

E(62)  - 

J  (W.  4»rrWb  +  wb*ddwb> 

From  Eq.  A20 

Vl  - 

1  -  GxWa 

(B4) 

W2  - 

-G2  “ 

(B5) 

Bl 


♦ 


With  =  0,  from  Eqs.  A35  and  A36, 


GlGlWa(t)RR  '  g1^RR  =  ^1 


(B6) 


GlGlwb^dd  +  GlG2^dd  *  ^2 


(B7) 


Write  G^  as  =  NQ/D,  where  Q  is  the  product  of  all  RHP  factors,  e.g.. 


G  «  K(s  +  a) (s  -  b) 
s2(s2  +  2^0)ns  +  co£) 


(B8) 


where  N  =  K(s  +  a) ,  D  *  s2(s2  +  2£WnS  +  oj^)  ,  Q  «=  (s  -  b)  .  Then  Eq.  B7  can 
be  written  as: 


AA.Sl. 

5  D 


Wa  ’  (<*W 


*1 


(B9) 


From  Eq.  A3 8 


wa  “  —  • 


-D  1 


QN  (4>rr)+ 


(Jr )(4>y(-M 

\Q  n/  \  D  /Vtjjj)  / 


-D  1 


QN  (<^r)  + 


•=•  <W 


Whence 


W,  -  1  -  GLWa 


i  _  d?  '  1 


*  W  ^rr>+ 


—  CA  1+ 
Q  ( ^ 


J  + 


1--2-  -1 


Q  (Orr) 


Q  (A  & 
L  Q  1111 


(BIO) 


(Bll) 


(B12) 


(B13) 


B2 


I 


From  Eq.  A38, 


W 


b 


-D 

QN 


1 

(4>dd)  + 


1  _  N  Q 

(4>dd>"  D 


(B14) 


whence 


Wv 


—  — 1  ■  Sr  *  G2<4)dd> 

QN  (4>dd)  Q 


Wn  =  ~ G2 


Q  1 


Q  (*dd)+ 


f'G2<W  + 


(B15) 


(B16) 


B.2  EVALUATION  OF  MEAN  SQUARE  QUANTITIES 

The  contribution  to  E(h|)  that  comes  from  that  part  of  he  which  is 
correlated  with  the  command  E(h^) R  is  given  by: 


4>RRw1w1  “ 


^rr  ~ 


■f  ■  w*s)+ 

Q  ,4,  V  + 

-  _  •('J’rr) 

-  ’  ^RR^ 

Q 

Q 

Q 

If  •  (4>rr) 


J  + 


4 


(B17) 


<t>rr  -  M[M]+  -  M[M]_  +  [M]+[M]_ 


(B18) 


where 


63 


(B19) 


r 


^RRWlWl 


f( 


+  [M]+[M]_J 


The  terms  marked  ^  do  not  contribute  to  the  integral  since  they  have 
poles  only  in  either  the  left  or  the  right  half  plane.  Hence: 


E(he) ^ 


J  <t>RRWiWx  -  J  (^rr  -  [M]+[M]_j 


(B20) 


This  can  easily  be  evaluated  by  noting  that  *  MM,  whence: 

E(he)R  =  f  ([M]+  +  +  ®.)  ~  M+M- 


(B21) 


j  ([tJ>j<M]+  +  [M]+[M]_  +  [M]_[M]+ 

+  [yyt£]_  -  [m]+[m]_) 


(B22) 


rif  <w+ 

((*W~ 

/_j°o  Q 

Q 

djw 


(B23) 


The  above  integral  is  simple  to  compute  because  typically  [  (Q/ Q)  ( 4>rr )  ] 
has  only  one  denominator  factor . 


of  he 


2 

The  second  of  the  two  contributions  to  E(he)  comes  from  that  part 

2 

which  is  correlated  with  the  disturbance,  E(hg)^,  and  is  given  by: 


*ddW2W2 


Put 


L 


*ddW2W2 


*>ddG2G2  "  G2  (<W 


q  G2(W 


-  G2  (4>dd)  + 


4-  G2(4)dd) 


G2  (4>dd) 


-2_  7.. 


G2^dd^ 


(B24) 


g2  “q"  (<*>dd> 


J  (lL  -  [L J/fL]+  -  [L]+[L]_  -  [L]+ 
-  +  [L]+[L]_  ) 

j  (lL  -  [L]+  [L]_ j 


[L] 


(B25) 


(B26) 


(B27) 


n  o  0 

[L]+[L]_  +  p<fL]+  +  [L]_[L]+  +  -  [L]+[L]_) 


/  tLMLJ. 


(B28) 

(B29) 


B5 


I 


Thus  the  remaining  contribution  to  the  mean  square  error  is: 


E(h~) 


f*/[’G2Q(*dd>‘ 

G2Q(4>dd)+ 

Ml  «  J 

+ 

Q 

(B30) 


It  is  also  of  interest  to  compute  the  mean  square  control  deflection 
associated  with  k  =  0.  This  is  given  by: 


■/ 


E(^)  =  (»atfa  +  VddV 


(B31) 


r\ 

The  contribution  to  E(<5^)  that  stems  from  the  component  of  6  correlated 
with  the  command  is: 


««2>R  ■  W.  '  /  (sH 


■w  (tW+ 


J+ 


(♦rr)' 


(B32) 


The  remaining  contribution,  which  stems  from  the  component  of  6  correlated 
with  the  disturbance,  is: 


E(62) 


♦ddwb“b 


f/sl 


l  NQNQ 


-2-  G2(0dd)+ 


■q-  G2<W' 


.) 


(B33) 


B6 


J 


t 


B. 3  COMPUTATION  OF  RESULTS:  MEAN  SQUARE  ERROR 


Writing  Eq.  B23  in  full: 


e(Vr 


-l-  p 

27Tj  J-y 


^rr) 


_S_ 


(4>rr) 


djw 


The  command  spectrum  is: 


(4>rr) 


818.9  s2 
(s  +  6.3)4 


Q  -s  -  b  b  +  s 

»  *  -  as  ■  ■ 

s  -  b  b  -  s 

1  .  2b(818 . 9) b2 

b  *  s  (b  +  6.3)4 


A. 


<h B? 


[b+s  .  818.9  s2 


b  "  s  (s  +  6.3)4 


E(he)^ 


^(818.9) 


2b3  \2  _  J_ 

(s  +  6.3)4 /  2b 


For  the  configurations  of  interest,  we  have: 


(B34) 


(B35) 


(B36) 


(B37) 


(B38) 


Conf igurat ion 
b 

E(he)2  ft2 


3GA 

16.9799 

21.944 


6GA 

21.7433 

17.0406 


9GA 

24.3630 

14.73 


B7 


Writing  Eq.  B30  in  full: 


E(he)‘ 


<*dd> 


G2  * 


G2^^dd)+ 


1  f2“  [^Q^dd)” 

+  l 

G2Q($dd)+ 

2*2  J-i-l  Q  J 

Q 

1 4 . 04  3 ( s  +  4.5143) 

(s  +  7.819)2 

R42S  +  R43 

_2_  = 

s(s2  +  R34S  +  R33) 

> 

Q 

R42S  +  R43  b  +  s 

14 . 043 (s - 

s(s2  +  R34S  +r35> 

1 

(s+7. 

dj  to 


b  +  s 
b  -  s 


(B39) 


(B40) 


(B41, 

B42) 


b  -  s 


s=b 


(R42b  +  R43)2X:  14 .043 (b  +  4.5143) 
Xb2  +  R34b  +  R35)(b  +  7.819)2(b-  s) 


(B43) 

(B44) 


E(he) 


28.0876(R42b  +  R43)(b  +  4.5143)|  { 

(b2  +  R34b  +  R35)(b  +  7.819)2  )  2b 


(B45) 


For  all  the  configurations  of  interest  we  have: 


Configuration 

3GA 

6GA 

9GA 

r42 

1.2276 

2.2904 

3.1404 

R43 

1.0017216 

3.31879 

6.155184 

r34 

1.7606 

3.2372 

4.4166 

r35 

14.8089 

37.79935527 

63.83748685 

b 

16.9799 

21.7433 

24.3630 

E(he)2,  ft2 

0.00012212 

0.00013692 

0.00014698 

AD-AQ96  950 
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B.4  COMPUTATION  OF  RESULTS:  MEAN  SQUARE  CONTROL  DEFLECTION 


Writing  Eq.  B32  in  full: 


E(62)' 


rj“A 

27r3  J_jco  NQ 


(*RR)+ 


D 

NQ 


(W 


dju) 


(B46) 


First,  let  us  evaluate  [  (Q/Q)  ] 


"S’  (<()RR)  + 


(<bs)+  ~ 


<'W  + 


(B47) 


b  +  s  818.9  s2  818.9  *  2b3  1 


b_s  (s  +  6.3)4  (b  +  6.3)4  b-s 


(B48) 


_818 

(b 


H-  - L- 7  IsW-— ^(S+6.3>4| 

~ s)  (s  +  6.3)4  (  (b+6.3)4  \ 


(B49) 


818.9 


(b-  s)(s+  6.3)4(b  +  6.34) 


|  (s3  +  bs2)  (b  +  6.34)  -2b3(s+ 6.3)4} 

(B50) 


Put  6.3  *  c.  Then  it  can  be  shown  that  the  (b-s)  factor  cancels,  giving: 


-818.9 j-2b3s3  +  [-b4  -  4cb3  +  6c2b2  +  4c2b  +  c4 ]s2 

_ +  [8c3b2  +  2c4b  Js  +  2c4b2  1 

(s  +  6.3)4(b  +  c)4 


(B51) 


Note  that  this  expression  is  third/fourth  order  in  s. 


( 


To  compute  E(6r) 


2 

R 


the 


following  quantity  is  required. 


_D_ 

NQ 


s(s2  +  2Su)ns  +  c4)(818.92b3s2+  •••) 

(s  +  a)(s  +  b)(s  +  6.3)4(b  +  6.3)4 


(B52) 


But  since  this  expression  has  a  numerator  which  is  sixth  order  in  s  and 
a  denominator  which  is  also  sixth  order  in  s  the  corresponding  integral 
is  unbounded.  This  is  a  consequence  of  the  particular  forms  of 
Gj  and  (4>rr)+  selected  here.  Within  the  limits  of  the  short-period 
airframe  dynamics  approximation  the  form  of  G^  is  fixed,  hence  D/NQ  must 
have  a  third  order  numerator  and  a  second  order  denominator.  To  obtain 
a  finite  E(6)^  it  is  thus  necessary  for  (<1>rr)+  to  have  a  denominator 
which  is  of  (n  +  2)th  order,  where  n  is  the  order  of  the  numerator  of 
^RR^+"  This  could  be  done  by  introducing  a  high-frequency  break  in  the 
<f>RK  spectrum. 

From  Eq.  B33: 


/JOO 

»2*dd»2 


But 


Wo  -  -Go  +  - ” — T" 

2  2  0  <+dd>+ 


-S-  <vw+ 


(B53) 


(B54) 


G2 


r42s  +  R43 
s(s2  +  250^3  +  a)2) 


r42s  +  r43 
s(s2  +  R34S  +  R35) 


(B55) 


<*dd>  + 


14.0438(s  4-  4.5143) 
(s  +  7 .81 9) 2 


(B56) 


BIO 


j 


(B57) 


(4>dd)+w2  *  -(^dd)+G2  +  ^ 


dd'  2 


G2(<W  + 


/  *dd"2»2  •  f  (*. 


ddG2G2  "  (^dd>  G2 


^dd^  G2 


r«.  -2. 
Q 


G2(*dd)+ 


G2(W 


.) 


(B58) 


/' 


^ddG2G2  “ 


^  G2(<t,dd)  + 


[c.c.]+  - 


■q*  G2(W 


[c.c.  ] 


"§"  G2(<W' 


J+ 


JL 


Q  G2((!>dd) 


(B59) 


where  [c.c.]  denotes  the  complex  conjugate  of  the  immediately  preceding 
term  in  [  ] .  Hence: 


/' 


/*(< 


^ddW2W2  "  I  l $ddG2G2  ~ 


"q-  G2(W" 


*  [C.C.].) 


(B60) 


To  evaluate  this  expression  note  that: 


»dd+G2 


14.0438(s +  4. 5143)  (RA2s  +  R43)  2nd  order 


(s  +  7.819)^(s)(s^  +  R3As  +  R33) 


(s) (4th  order) 


(B70) 


B 1 1 


1  (  2b(-RA2b  +  R43)(-b  +  4.51A3)  )  =  constant 

S  +  b  ((-b)(b2-R34b  +  R35)(-b  +  7.819)2  j  lst  order 

(B71 

s  unbounded,  due  to  the  free  s  in  the  denominator  of 

<f’dd'r^2  ‘ 

To  obtain  a  finite  E(tS2),  it  would  be  necessary  to  modify  the  low- 

Q 

frequency  end  of  the  Dryden  atmospheric  turbulence  spectrum,  i.e.,  to 
replace 


-Q-  62(*dd>  ] 
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Hence  E(62),  i 
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(W 


14 . 0438 (s  +  4.5143) 
(s  +  7 .819) 2 


by 


<w 


s  .  14 . 0438 (s  +  4.5143) 
(s  +  h)  (s  +  7.819)2 


(B72' 


(B73; 


where  h  is  a  small  frequency  (e.g.,  0.1  rad/ sec) . 
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